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1 . 1  MBS  Srowth  and  Characterization* 

During  this  period,  a  considerable  effort  was  placed  on  the  preparation 
and  characterization  of  single  and  multiple  period  modulation  doped 
( A1 ,Ga) As/CaAs  heterostructures,  A  theoretical  study  of  these 

he taro structures  was  also  undertaken  and  a  good  agreement  with  experiments  was 
obtained.  Segregation  effects  of  Sn  in  GaAs  grown  by  MBE  was  also 
investigated  in  detail. 

Single  and  multiple  period  modulation  doped  structures  with  mobilities 
several  times  those  reported  by  other  laboratories  ware  prepared.  Transport 

properties  of  these  structures  were  studied  at  low  and  medium  electric  fields 

and  at  room  and  cryogenic  temperatures.  The  maximum  room  temperature 

mobility,  3,400  cm" /Vs,  was  obtained  when  the  AlAs  mole  fraction  was  about 
3'5^.  In  addition,  the  single  period  structures  so  far  appeared  to  exhibit 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 

(U.S.  Army,  U.S.  ?Javy,  u.S.  Air  Force")  under  contract  N0001 4-79-C-0424 ,  and  by 
the  Air  Force  Office  of  Scientific  Research  under  contract  APOSR  30-0034. 
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l^r^er  room  tampsr^tur?  mobilities  then  the  multiple  periol  ones.  At  '<  eni 
!  ■)  however,  oompereble  results  were  obteinei.  77  '<  mobilities  of  100,000 
om- ^Vs  enl  10  '<  mobilities  of  about  200,000  -  210,000  err 'Is  were  obteinei 
both  with  single  end  multiple  period  structures.  If  the  AlAs  mole  fraction  is 
made  larger  than  about  the  direc  t- indirec  t  transition,  the  electron 

mobilities  steadily  decline  up  to  the  AlAs  mole  fraction  of  lOOt. 

Secondary  ion  mass  spectroscopy  (STMS)  has  been  used  to  investigate  the 
accumulation  of  3n  at  the  growing  film  surface  and  the  corresponding  3n 
depletion  near  the  film- substrate  interface  during  IIBE  deposition  of  3n-doped 
SaAs.  3n  surface  concentrations  more  than  two  orders  of  magnitude  larger 

than  the  steady  state  3n  doping  level  have  been  observed.  in  films 

grovm  at  approximately  constant  As^Sa  surface  concentrations  was  fo^ani  to 
increase,  while  remained  constant  with  increasing  growth  temperatures  for 
”^3  i  530°C.  Both  and  decreased  at  higher  Tg  values.  These  results, 
as  well  as  those  reported  by  previous  investigators,  have  been  explained  by  3n 
surface  segregation  luring  deposition  in  which  the  segregation  is  driven  by 
the  relief  of  strain  energy  due  to  the  incorporation  of  over-sized  3n  atoms. 
Galculatei  Sn  profiles  were  shown  to  provide  a  good  fit  to  experimental  data. 

'  .  2  Transport  ? roperties  of  'letero  junction  Layers* 

Transport  properties  of  modulation  doped  structures  uo  to  an  electric 
field  strength  of  2  cm  were  investigated.  The  differential  electron 
mobility  drops  cuite  rapidly  up  to  200  '/''■cm  and  then  stays  about  constant  up 
to  2  <V-'cm.  Higher  electric  fields  cause  negative  differential  resistance 

effects  -/ihich  can  result  in  nonuniform  electric  fields.  Tnce  the  field  is 
nonuniform,  no  reliable  u(T)  data  can  be  obtained.  The  electric  fields  in  1  u 
gate  normally-off  ?3Ts  are  estimated  to  be  about  2  hV/cm;  therefore,  our 
results  should  be  applicable  to  predicting  device  performance.  Structures 
with  an  extremely  high  low-field  mobility  showed  substantial  decrease  with 
electric  field,  particularly  at  low  temperatures.  Generally  the  room 
temperature  mobility.  3,i00  cm- ''Vs.  was  maintained  up  to  2  xV'em.  However,  at 
■<.  mobility  at  2  </''cm  was  about  2c-dT'o  of  its  zero- field  value.  The  '  2  < 

*  This  wor'x  was  sunoorted  by  the  Air  Force  Office  of  Scientific  Research  under 
contract  AFOSR  30-0031,  the  Office  of  Haval  .Research  under  contract 
'I0001  l-''o-0-0306,  and  by  the  Army  Research  Office  under  contract  OAAG 
29-30-0-0011. 
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mobility  at  2  !cY/cm  was  about  15-20o  of  its  zaro-fieli  value-  These  results 
iuiioate  that  a  differential  mobility  of  about  25,000  cm^^Ys  at  2  'icY/cm  with 
an  electron  concentration  of  about  can  be  obtained  at  7"^  X.  This 

represents  an  intrinsic  speed  improvement  of  5  and  a  speed-power  product 
improvement  of  25  if  normally-off  FETs  were  fabricated  and  operated  at 
Room  temperature  operation  can  result  in  an  improvement  of  50^  and  about  200o 
in  speed  and  speed-power  product,  respectively.  3y  "improvement"  we  mean  to 
compare  it  to  a  conventional  normally-off  CaAs  FET  using  a  bulk  channel  layer. 

We  have  done  considerable  work  during  this  period  on  the  real- space 
electron  transfer  effect,  in  which  electrons  in  GaAs-AlGaAs  he tero structures 
transfer  at  high  fields  from  the  CaAs  to  the  AlGaAs  layers.  The  resulting 
negative  resistance  as  been  simulated  by  Monte  Carlo  techniques  and  has  been 
observed  experimentally.  This  effect,  which  has  great  potential  for  device 
applications,  is  discussed  further  in  Unit  4-. 

1  . 3  Theory  of  Surfaces  and  Retero junction  Interfaces* 

The  deep  levels  associated  with  surface  defects  in  compound 

semiconductors  have  been  studied  theoretically;  the  results  indicate  that 
Fermi-level  pinning  and  Schottky  barrier  formation  is  produced  by  surface 
antisite  defects  in  most  cases,  rather  than  by  surface  vacancies  or  subsurface 
defects.  We  also  have  investigated  the  possibility  of  designing  Schottky 

barrier  heights,  and  barrier  heights  at  he tero junctions  and  in  superlattices, 
via  the  deposition  of  selected  Impurities.  Predictions  of  intrinsic  surface 
states  and  surface  core  excitons  in  III-Y  and  II-YI  semiconductors  were  found 
to  be  in  excellent  agreement  with  the  available  experimental  data.  Finally, 

it  was  found  that  impurities  that  are  shallow  donors  or  acceptors  in  the  bulk 

can  become  deep  traps  at  semiconductor  interfaces. 


*  This  work  was  supported  by  the  Joint  Ser>/ice3  Electronics  Program 
(U.3.  Army,  U.S.  Navy,  U.S.  Air  Force)  under  contract  NCD001 4--''9-j-Dl2  4-. 
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2 . 1  Intro iuction 

The  prixary  goal  of  this  research  is  to  evaluate  iae  potential  of 
high-purity  In?  and  InSaAsP  for  optical  and  xicrowave  ievice  applications. 
Besides  investigating  the  various  growth  paraxeters  and  techniques  that  can 
influence  the  purity  of  the  grown  layers,  an  ixportant  part  of  this  work  is 
the  characterization  of  the  purity  and  crystalline  quality  of  the  samples 
groTO.  The  total  Lmpurity  concentration  is  studied  through  Rail  coefficient 
measurements  at  30Q°K  and  at  77°K.  The  origin  and  influence  of  growth 
conditions  on  acceptor  concentrations  are  studied  through  low  temperature 
photol'jminescence  measurements.  Deep  levels,  and  the  influence  of  different 
growth  parameters  on  the  concentrations  and  energy  levels  of  these  centers, 
will  be  studied  through  DLT5  measurements.  The  crystalline  quality  and  degree 
of  lattice  match  are  studied  through  x-ray  diffraction  measurements.  The 
residual  shallow  donor  impurities  are  studied  using  far  infrared  Fourier 
transform  measurements  of  photothermal  ionization  spectroscopy.  This 
technique  is  being  used  to  study  the  residual  donor  impurities  in  high-purity 
laAs  and  InP.  Far  infrared  photoconiuc tivity  measurements  on  InSaAsP  and 
InDaAs  have  been  used  to  determine  the  effective  mass  in  these  matarials,  and 
this  technique  may  be  useful  for  studying  alloy  disorder  in  these  materials. 
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2.  SEMICONDUCTOR  MATERIALS  AND  DEVICES 
2 . 2  Hydride  Vapor  Phase  Crowth  of  High  Purity  InP* 

The  vapor  phase  system,  constructed  vrith  funds  from  ONE,  is  now  fully 
automated  and  under  computer  control  for  added  reproducibility.  Previously 
high  purity  GaAs  (^77  ^  “  38,000  cm^/V-sec,  077  '^  =  1.T  x  10^^  cm"^ )  had  been 
grown  in  this  system  but  comparable  InP  had  not.  Difficulties  encountered  in 
growing  InGaAs  lattice  matched  to  InP  led  to  a  more  detailed  and  systematic 
analysis  of  the  system's  growth  parameters  for  InP. 

Prior  to  growth,  the  substrates  are  etched  in  situ  with  a  0.13“?  HCl  in  H2 
mixture  to  provide  a  clean  surface  upon  which  to  initiate  growth.  In  addition 
the  Arsine/Phosphine  line  has  been  extended  further  down  into  the  gas  stream. 
The  growth  rate  has  increased  to  the  point  where  it  is  possible  to 
substantially  decrease  the  total  gas  flow  (by  nearly  90“?).  These  lower  flows 
result  in  much  less  downtime  for  cleaning  purposes  and  a  substantial  saving  in 
the  amount  of  source  material  used.  The  InP  grown  in  this  system  is  now 
comparable  to  that  grown  by  any  technique,  with  the  best  material  to  date 
having  a  liquid  nitrogen  mobility  P77  ^  =  70, XO  cm-ZV-sec  and  carrier 
concentration  077  k  =  5  x  1  O' cm“^ . 

2 . d  LPB  G rowth  of  InGaAsP** 

High  purity  InP  and  InGaAsP  epitaxial  layers  are  essential  for  wide 
depletion  width  and  low  capacitance  PIN  and  avalanche  photodiodes  and  for 
materials  characterization  in  general.  Also,  for  small  band  gap  InGaAsP 
alloys,  high  purity  material  is  required  to  avoid  Zener  breakdown  rather  than 
avalanche  breakdown  in  p-n  junctions.  Methods  for  obtaining  high  purity  InP 
and  InGaAsP  epitaxial  layers  on  (IOO)  InP  substrates  using  liquid  phase 
epitaxy  have  been  developed  and  this  high  purity  material  has  been 
characterized  through  Hall-effect  measurements  and  far  infrared  photothermal 
ionization  measurements.  Samples  with  net  carrier  concentrations  as  low  as  1 
X  10''’  cm"^  for  InP  and  2  x  1 0^  ^  cm"^  for  In^^-jGa^^^yAs,  with  corresponding 
liquid  nitrogen  temperature  electron  mobilities  of  67,500  and  57,900  cm-/V-s, 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
(U.S.  Army,  U.S.  Navy,  and  U.S.  Air  Force)  under  Contracts  DAAG-2  9-'^3-G-j01  5 
and  NX01  4-79-0-0424  and  by  DARPA  Contract  NX001  4-77-C-000S6 . 

**  This  work  was  supported  by  the  Naval  Research  Laboratory  under  Contracts 
NX1 75-73-G-01 29  and  NX1 75-79-C-01 34,  and  the  Joint  Services  Electronics 
Program  (U.S.  Army,  U.S.  Navy  and  U.S.  Air  Force)  under  Contracts 
DAAG-29-73-C-OOI  6  and  NX01  4-79-C-0424. 
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respectively  have  been  routinely. 

Both  the  grovrth  solution  baking  technique  and  the  growth  temperature  have 
an  influence  on  the  epitaxial  layer  purity.  Baking  times  of  24  to  43  hours 
are  required  to  reduce  the  carrier  concentration  to  these  levels. 
Incorporation  of  small  amounts  of  water  vapor  (<  .1  ppm)  in  the  hydrogen 
atmosphere  in  an  attempt  to  reduce  silicon  transport  from  the  quartz  tube  to 
the  gro'wth  solution  during  baking  has  not  improved  the  purity,  and 
incorporation  of  larger  amounts  (>  1  ppm)  increases  the  carrier  concentration 
and  reduces  the  mobility.  An  increase  in  the  growth  temperature  has  improved 
the  purity,  especially  for  InOa.As,  presumably  due  to  a  decrease  in  the  silicon 
distribution  coefficient  as  the  growth  temperature  is  increased. 


2 . 4  Fourier  Transform  opectroscopy  of  Shallow  Donor  Levels* 

Shallow  donor  energy  levels  have  been  measured  using  photo  thermal 
ionization  spectroscopy  in  high  purity  LP2,  Hydride  VPS,  and  bulk  InP,  grown 
in  our  laboratory  and  else-where.  This  .echnique,  which  involves  the 

measurement  of  photoconductive  response  versus  incident  photon  energy  at  low 

temperature  (4  '<)  and  in  a  high  magnetic  field  (_£  5.5  T)  has  also  been  applied 

to  high  purity  LPS  InSa.As  and  InGaAsP  samples  lattice  matched  to  InP.  The 

minicomputer  based  Fourier  transform  spectrometer  facility  was  set  up  with  the 
support  of  the  Naval  Research  Laboratory. 

In  the  LPS  and  Hydride  ’/PS  InP  grown  in  our  laboratory  and  the  L?S  In? 
grown  by  L.  ?.  Eastman  at  Cornell,  the  sample  purity  was  sufficient  that  peaks 
corresponding  to  different  donor  species  could  be  resolved  in  the  phototherxal 
ionization  spectra.  Though  it  is  not  yet  possible  to  identify  these  donor 
species,  the  two  donors  present  in  the  Hydride  VPS  material  were  both  common 
to  the  Illinois  LPS  material,  and  one  of  these  was  also  the  dominant  donor  in 
the  Cornell  LPS  material,  A  third  donor  species  may  also  be  present  in  the 
Illinois  LPS  material  at  low  concentration.  It  is  reasonable  to  speculate 
that  a  donor  common  to  both  LPS  and  VPS  InP  may  be  silicon,  incorporated  from 
the  hot  quartz  grcwth  reactors,  especially  in  view  of  the  high  distribution 
coefficient  for  silicon  in  InP, 

*  This  work  was  supported  by  the  National  Science  Foundation  under  Contract 
N3F  DHR-'^''-25  239,  by  the  Office  of  Naval  Research  under  Contract 
NOGOt  4-27-C-Co55  and  by  the  Joint  Services  Electronics  Program  ('J.S.  Army, 
■J.3.  Navy,  and  N.S.  Air  Force)  'under  Contract  NOOOl  4-"'9-C-0424. 
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2.  SEMICONDUCTOR  MATERIALS  AND  DEVICES 

The  magnetic  field  dependence  of  the  peak  positions  in  these  spectra 
allows  accurate  determination  of  bottom  of  the  band  effective  mass  (ni*), 
donor  ionization  energy  R  ,  and  low  temperature  static  dielectric  constant 
Values  of  these  important  materials  parameters  have  now  been  obtained  for  InP, 
In^^^Oa^^Y^s  and  For  these  materials,  the  respective 

values  determined  for  these  parameters  are: 

(m*)  =.  .00796  ^  .005,  -0437  ^  .001  ,  and  .0635  ^  .001 

R*  =»  7.4  .17  meV,  4-30  .17  meV,  and  5.23  .17  meV 

eg  =■  12.06  .13,  13.5  ^  .5,  and  13.4  ^  .5. 

The  peaks  in  the  spectra  of  the  ternary  and  quaternary  alloy  samples  are 
much  broader  than  would  be  seen  in  GaAs  or  InP  of  similar  liquid  nitrogen 
temperature  carrier  concentration  and  Hall  mobility  (iijy  =  3  x  1 0' ^  cm~^ , 

=  53, XO  cm^/v-s).  A  mechanism,  not  present  in  the  binary  alloys,  where 
compositional  alloy  disorder  results  in  a  different  local  alloy  environment 
for  different  donors,  could  account  for  the  excess  broadening  seen  in  the 
ternary  and  quaternary  spectra.  The  degree  of  this  broadening  might  then  be 
interpreted  as  a  measure  of  this  compositional  disorder. 
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3. '  Excitation  ani  3 tudy  of  Highly-Excited  Vibrational  States  of  Molecules 

3 . ' . ’  Eirection  and  Significance 

tuantuTi  electronics  hai  its  roots  in  the  xicroi^aTe  3pe:troscoDy  of  lov 
lyh-ng  energy  le'/els  of  xolecular  gases.  Eptically  p'mpad  far  IE  lasers  aare 
eiiployei  the  funiaaiental  ’/-ibrational-rotational  energy  levels  of  -,he 
xclecules.  Mot  banis  ani  excited  states  of  xolecules  io  to  the 
" guasi- continuum"  have  been  little  used  in  quantum  electronics  ani  the 
spectroscopy  ani  state  xinetics  are  relatively  unexplored. 

Ic  is  the  thrust  of  this  vorx  to  study  the  Nineties  of  the  populations 
and  pclarications  of  xolecular  states  up  to  the  ’’ guasi -continuux"  for 
potential  uses  in  freq^uency  tunable  lasers,  nonlinear  optics  and  chexistry. 
To  iapleaient  the  study  requires  the  ievelopment  of  state-specific  excitatiica 
and  diagnostic  techniques,  the  initial  stage  of  the  fo-h. 

*  This  Yor'x  ^as  suapicrted  by  the  Toint  lervices  Electronics  ‘^'•ograt. 
,‘1.3.  irry,  '1.3.  Mavy,’ 'M.  3.  \ir  7orce')  'under  contracts  TIM  "> "^-'’3-3 -TT1  c  and 
MTTTl  and  by  contracts  M0TT1  i-3T-TTT3T,  MMi  vr,-!  'Oa,  'na-  qpg 

3T-TM3''3,  q7T3R  '’a-OMM,'  and  TTE  3aadii  PT-ITT033. 
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3.1.2  Summary  o f  ?  rogreaa 

A  pulsed  multimode  GO^  laser  has  been  used  to  sxoite  a  hot  band 
transition  in  room  temperature  OGS  gas  to  a  state  .fith  three  q^uanta  of 
ribrational  excitation.  A  transient  microwave  gain  was  obser/ed  and  is  a 
measure  of  the  excited  state  population  difference  between  1  doublet  states. 
This  inherently  small  signal  has  been  observed,  without  signal  averaging,  by 
using  a  microtwave  reflection  bridge  -which  is  transiently  unbalanced  by  an 
infrared  pulse  exciting  the  gas  sample  in  part  of  the  bridge.  An  infrared 
pulse  truncation  shutter  and  lower  noise  micro-wave  detection  system  are  being 
added  for  direct  obser/ation  of  the  excitation  relaxation. 

3 . 2  Tnergy  Transfer  ? recesses  in  Electronically  Excited  Atoms  and  Molecules 

3.2.1  LiOw-?ield  Ionization  Coefficients  of  Rare  Eases 

The  elsctron  ionization  and  excitation  coefficients  for  Ar,  Kr,  and  le 
have  been  measured  in  the  low  E/p  region  (1-12  V/cm  torr)  using  a  drift-tube 
apparatus.  Discrepancies  bet-ween  previous  measurements  have  been  resolved  and 
sho-wn  to  be  due  to  photoelectric  amission  attributed  to  the  rare  gas 

metastables.  Using  this  information  a  revised  set  of  self-consistent  near 
threshold  inelastic  cross  sections  for  each  rare  gas  studied  -was  determined. 

3.2.2  '-Cinetics  of  Interhalogen  Diatomics 

The  spectra  of  the  3  A^II^  states  of  17  and  IDl  -were  obtained  and  for 
the  former  sho-wed  clearly  resolvable  rotational  structure.  The  lifetimes  of 
the  E  states  of  17  and  ICl  were  also  estimated  from  pulsed  discharge 

experiments.  The  vibrational  self  quenching  rate  of  the  IDl  state  -was 

estimated  and  also  the  argon  quenching  rate  of  this  same  state  -was  measured. 

3.2.3  Multip'noton  Interactions 

In  an  unrelated  study,  -we  have  been  exploring  some  effects  associated 
with  the  multi-photon  nature  of  near  resonant  light-matter  interactions 
involving  more  than  one  field  and  a  sat  of  coupled  energy  levels.  'Ve  have 
developed  both  algebraic  algorithms  and  a  graph-algebraic  approach  fer  solving 
this  class  of  problems  in  the  adiabatic  limit.  These  techniques  were  then 
applied  to  the  specific  case  of  resonant  second  harmonic  generation  in 
three-level  systems,  an  example  of  -which  occurs  in  magnetic  resonant  cases. 
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T'nia  stuiy  shova  thit  it  ia  possLbla  to  bave  i  vary  bi?h  coaveraioi, 
than.  jOl,  of  p'j.xp  aaar^y  into  tha  aaconi  harnonia  aniar  aalactai,  bat  not  too 
rastriativa.  toniitiona.  Tna  raaalta  of  this  thaoratiaal  vor<  vill  appaar  in 
a  fo '.’the  ox  in-?  pablieation.  •fa  ara  no-rf  applying  tha  finiings  of  thia  ataly  to 
apacifie  ayataxa  in  orhan  to  latarxina  tha  feasibility  of  fatare  axpanixanta. 

5 . 3  3xaitai  Itata  lhaxiatry  in  laaes 

Tha  raaearah  pro.] ant  on  Txaitai  3tata  Taaxistry  in  Oases  has  grovri  iarin? 
tha  past  yaar  frox  the  i.nitial  iisnharge  wr'-c  on  plasxa  etching  to  inoluia 
plaaxa  annealing  ani  gro-,rth  'by  electron  ixpact'l  ani  xultiphoton  iissociation 
of  sexiconi  ac  to  r  bearing  gases  so  as  to  gro-.v  thin  filxs  of  Silicon  or 
larxaniax  'Ll’D'l, 


3.'.'  ^lasxa  Ttching 

Tha  previous  report  xantionai  tha  aviianca  for  tha  role  of  negative  ions 
in  the  etching  of  3i  or  .Sil-j  yith  iischarges  in  'I?y  Va  hal  foani  that  tha 
etch  rate  of  these  xaterials  vas  greatly  anha.ncal  vhen  tha  saxple  ^as  biasel 
to  a  feif  tans  of  volts  positiva  <ith  respect  to  floating  potential.  Va  have 
since  correlatal  this  bahivior  to  Langxuir  probe  charac  t eri.st i 'is  fnieh 
iniicata  a  xaesive  negatively  chargai  carrier.  Tnis  behavior  i.e  only  aviiant 
for  lischarges  in  'ITg;  -iischarges  in  the  xore  conventional  xixturas  of 
lo  not  exhibit  this  behavior.  Tiis  is  consistent  .^i  th  the  t'lct  that 
'17-^  has  a  xuch  larger  attachxent  cross  section  than  the  other  gases  1' 

!•  ^  a 

Thase  results  have  bean  acceptei  fo'^  publication  _3.o|.  ‘Jnfortunately,  tha 
iiantity  of  tha  ion  is  still  un'^noTO  ani  avaits  the  installation  of  tha  xass 
scactroxetar  on  tha  etching  systax.  Tnis  is  scheiulei  for  coxplation  in  the 
forthcoxing  year. 

’.■^.2  ^lasxa  \nnealing 

In  this  phase  of  nhe  ,^or;<,  va  have  sho.^  that  one  can  use  a  aixpl  e 
abnormal  iischarge  in  a  low  pressure  gas  'typically  heliax  at  abr-ut  xtor”' 
to  obtain  a  larga  area  alactron  baax  which,  in  turn,  can  be  usei  to  anneal 
large  areas  of  ion  ixplantal  sexiconluc tor  xatarials  ,1'"'.  It  has  been  '■  ir 
experience  chat  tha  saxplas  annaaiai  in  this  xa.nnar  for  tixas  less  than  ■' 
seconls  have  sheet  resistivities  consistently  less  than  those  anneal ei  in  an 
oven  far  53  xinutes.  -it  the  present  tixe  we  are  attexpting  to  infer  the 
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t9tperatur<?  of  the  satnpla  by  neaaurlng  tha  change  in  optical  reflactivity  as 
tha  sample  is  irraliatad  by  the  electrons.  ihis  is  a  variation  on  the 

technique  lescribei  by  Olson,  Xikorowski ,  Roth  and  Hess  [12]. 

3.3.3  Laser  Thin  ?ilni  Orowth 

The  objective  of  this  research  program  is  to  investigate  the  chemistry 
and  optical  properties  of  excited  atoms  and  molecular  radicals  produced  by 

ultraviolet  photo  dissociation  and  the  interaction  of  these  excited  specias 
vith  a  surface  to  generate  a  thin  film.  Specifically,  an  ultraviolet  excimer 
laser  is  being  used  to  photodissociate  inorganic  molecules  that  contain  a 
semiconductor  atom  (such  as  3iH^,  3eH^  or  SiCl^) .  The  spectroscopy  of  these 
laser-generated  plasmas  as  veil  as  the  electrical  and  chemical  properties  of 
the  resulting  Se  or  3i  films  are  under  investigation. 

To  date  amorphous  polycrystalline  Si  and  Se  thin  films,  around  1  cm-  in 
area,  have  been  grown  by  photolytic  laser  chemical  vapor  deposition  'uZTi)  at 
room  temperature.  Dilute  SiH^,  (CH-^j^^Si,  or  SeH^  mixtures  in  an  ine^'t  carrier 
gas  were  pho todissociated  using  ArP  (193  nm)  or  KrP  (’4-3  nm)  lasers.  Film 
thicknesses  of  up  to  D.5  um  were  obtained  on  SiO-j  substrates  at  deposition 
rates  ranging  from  1  to  ID  A/sec  for  active  gas  partial  pressures  between  3 

and  3D  torr  and  incident  laser  intensities  up  to  25  T/f/cm-.  Absorption 

coefficients  of  as-deposited  De  and  Si  layers  were  2.S  x  1  D^  cm*’  and  3.'^  x 
t  D^  cm”*,  respectively,  and  sheet  resistivities  _>  !  D’’ /  □  have  been  measured. 
Incubation  times  and  initial  growth  rates  have  been  determined  situ  as  a 
function  of  'JV  photon  flux  from  He-Ne  laser  absorption  measurements.  A1 
concentrations  of  up  to  several  atomic  percent  have  been  incorporated  into  De 
films  by  the  simultaneous  photodissociation  of  DeHj^  and  (DH^)^AI. 

Finally,  the  chemical  characteristics  of  these  films  are  being  studied  by 
infrared  absorption  Auger  spectroscopy,  and  I-ray  scattering. 

References 

’rj]  I.  D.  J.  Thynne,  J.  Phys.  Chem.  2^,  1336  (iDoD). 

(R2I  1.  L.  Dlson,  3.  A.  Xokoro-wski,  J.  A.  Roth,  and  L.  D.  Ress,  "Investigation 

of  Laser  Induced  Solid-^hase  Crystal  3ro-<th  by  Time-Resolvei  Dptical 

Reflectivity,"  Taper  ThC-3,  CLSD  '31  ,  Washington,  D.C.,  June  (1931  ). 
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i. 1  Introduction 

This  research  invol/es  stuiy  of  basic  properties  of  ser.icoTiuctors, 
Tiethois  of  ievice  nrooessing,  ani  new  ievice  concepts,  3oth  theoretical  ani 
eKperinental  lethois  are  eTiployel  in  each  of  these  categories.  './e  are 
examining  a  variety  of  hot  electron  phenomena  and  their  effects  on  present  anl 
future  ievice  performance,  especially  in  connection  ■fith  moiulation  ioping. 
We  are  studying  ion  implantation  ani  annealing  of  31  ani  III -7  compounds, 
including  laser  ani  electron  beam  processing.  These  experimental  studies 
include  examination  of  deep-level  impuritias  and  defects  arising  f-'OT. 
implantation  ani  annealing.  Several  aspects  of  this  jrorx  a^e  done  la 
collaboration  with  other  units,  particularly  the  studies  of  naterials  grc vn  by 
molecular  beam  epitaxy  '433,  'Jnit  t). 
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4. •  2  Laaer  and  Electron  Beam  Annealing  of  Semiconducto ra* 

'lie  have  performed  both  theoretical  and  experimental  vrorlc  on  pulaed  laaer 
annealing  of  amorphized  Si  and  awept  line  electron  beam  annealing  of  Si  and 
SaAa.  In  collaboration  -/rith  T.  T.  Verdeyen,  ^e  have  alao  employed  gaa 
liacharge  annealing  to  proceaa  implanted  Si. 

In  atudying  pulaed  laaer  melt  annealing  of  Si,  ve  have  examined  heating 
by  band-to-band  excitation  (Ndrglaaa,  1.06  pm)  and  free  carrier  abaorption 
(-0,,  10.6  ym)  .  The  techniq^ue  of  finite  difference  equationa  vaa  uaed  to 
aolve  the  one-dimenaional  heat  conduction  equationa  for  the  aolii  and  liquii 
phasea,  talcing  into  account  the  temperature  dependence  of  all  the  parameters 
involved,  eapecially  the  aubatrate  parametera,  to  give  a  realistic  model.  In 
the  CO^  case,  ve  have  calculated  the  temperature  variation  of  free  carrier 
abaorption  in  n-3i.  'Ye  have  taken  into  account  acoustic  ieformation  potential 
scattering,  optical  deformation  potential  scattering,  and  ionized  impurity 
scattering.  The  ieformation  potentials  vere  adjusted  to  fit  the 
experimentally  obaer/ed  values  at  300°K.  'lie  have  investigated  the  ’’thermal 
runaway"  behavior  and  calculated  the  threshold  power  density  to  cause  thermal 
runaway  and  surface  melting  of  Si  during  CO,  laser  annealing.  For  the  pulsed 
Nd:glas3  laaer  annealing,  the  melting  model  ia  supported  by  time-reaolved 
reflectivity  measurements  using  a  He-Ne  and  an  Ar  laser.  -Ye  uaed  SE'i 
channeling  patterns  to  study  the  quality  anl  orientation  of  crystalline 
regroTTth.  To  investigate  the  electrical  activation  we  letermined  the  change 
in  sheet  resistance  of  the  implanted  material  before  and  after  a  laser  pulse. 
In  each  case  we  fcund  excellent  reordering  of  the  amorphized  Si. 

Successful  recryatallization  of  amorphous  BF,*  ion- implanted  Si  has  also 
been  demonstrated  by  the  swept  line  electron  beam  ^ SLE3)  technique.  High 
percentages  of  electrically  active  boron  were  observed  with  a  minimum  of 
spatial  redistribution.  Active  boron  was  also  found  in  the  amorphous- 
crystalline  transition  region,  an  effect  not  obser’y-ed  in  low  temperature 
furnace  annealed  material.  We  are  currently  employing  the  SLSB  method  in 
recryatallization  of  deposited  Si  layers  over  SiS,. 

*  This  work  was  supported  by  the  Joint  Ser/ices  Electronics  Program  'H.S. 
Army,  U.3.  Navy,  and  IJ.3.  Air  Force)  under  Contract  N00S1  by  the 

Army  Research  Sffice  under  Contract  DAAC  29-BO-C-001 1 ,  and  by  the  Office  of 
Naval  Research  under  Oonract  NOOOI  4-''6-C-0306. 
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-t.  SSMrGOND’JGTOR  PRYSIGS 

1.3  Ituiles  of  Traasporr  ^ ropertiaa  in  Semi conductor 3* 

We  have  iavestigalai  carrier  transport  in  aaai'onduotor  hateroj 'inc tion 
layers  and  developed  the  theory  of  mobility  in  modulation-doped  structures. 
Ghe  results  of  this  theory  vere  compared  ’ifith  emperiments  and  found  to  ^ive 
excellent  agreement  in  certain  temperature  and  dopin?  ranges.  Ve  have 
continued  to  uork  on  the  concept  of  real  space  transfer  and  have  shoim 
experimentally  the  validity  of  this  concept  'the  spill  out  of  electrons  from 
potential  vails) .  We  believe  that  this  effect  has  high  device  potential  and 
also  many  Interesting  basic  consequences  for  the  transport  in  layered 
structures  and  superlattices.  Ghe  affect  can  also  be  generalised  as  a  real 
spacs  analog  to  a  x-space  effsct  (the  Gunn  effect''  and  this  generalization 
leads  to  several  other  phenomena  such  as  "resonant"  impact  ionisation  because 
of  the  band  edge  discontinuity  in  hetero.iunction  layers. 

We  have  expanded  our  Wonte  Garlo  program  '^including  the  band  structure  as 
calculated  from  the  empirical  pseudopotential  method)  to  include  transport  in 
silicon  and  GaAs.  We  vere  able  to  calculate  ionization  rates  in  these 
materials  vithout  too  many  assumptions.  Ghese  calculations  include  '  fo"  the 
first  time)  the  population  of  more  than  one  conduction  band.  'We  have  shovn 
that  the  hig'ner  conduction  band  is  crucial  in  .silicon  at  high  electric  fields 
as  they  are  currently  reached  in  small  devices. 

We  have  continued  to  develop  numerical  programs  to  model  transnort 
phenomena  in  various  semiconductor  devices  and  are  currently  including 
transient  phenomena  ( overshoic  t  ^uniershcoz  of  velocity''  in  the  operation  of 
charge  couoled  devices.  Part  of  our  theoretical  vorlx  vas  devoted  to  the 
possibilities  of  ballistic  transport  and  to  the  field  dependence  o.f 
net srolayers  at  lov  electric  fields. 


*  This  Yorx  vas  supported  by  the  Joint  Services  Electronics  'Program 
Army,  'J..".  ''lavy.  and"  'J.c.  .Aar  Pores)  wonder  Gontract  'JGGGi  G -HI  l ,  by  the 

Gfflce  of  haval  Psssarch  under  Gontracts  'IT)G'  l-‘'S-G-GhGb  and 
'IGGG' l-d'G-G-G^o3,  and  by  the  .Army  Research  Gffice  under  Gontract  GAAG 
oa-ph-G-GGl  ’  . 


4.  SEMrCOtfDUCTOR  PH'XBICS 
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4 • 4  Materials  Studies  of  InP* 

As  a  part  of  our  continuing  studies  of  ion  implantation  in  compound 
semiconductors,  we  have  studied  InP  carefully  during  this  period.  The 
redistribution  of  the  compensating  dopants,  iron  or  chromium,  in  semi- 
insulating  indium  phosphide  has  been  studied  using  secondary  ion  mass 
spectrometry  (3IMS).  Annealing  vrith  silicon  nitride  encapsulation  results  in 
impurity  accumulation  within  the  first  1000  A  of  the  surface  followed  by 
depletion  extending  to  a  depth  of  ~  i  ym.  Profiles  resulting  from  the 
implantation  of  "neutral"  elements  (He,  3)  exhibit  accumulation  at  the  surface 
and  also  accumulation  at  the  projected  range  peak.  The  profiles  can  be 
explained  in  terms  of  gettering  of  the  compensating  dopant  to  defect- rich 
regions . 

We  have  examined  polycrystalline  and  single  crystal,  liquid  encapsulated 
Czochralski  (L3C)  grown  InP  for  evidence  of  boron  contamination  using  STMS  and 
photol'uminescence .  Precipitates  of  boron  or  a  boron  compound  have  been  found 
in  In?  grown  by  the  LSC  method  with  boric  oxide  (3^0-^)  encapsulation  and 
pyrolytic  boron  nitride  (pBN)  crucibles.  The  density  of  precipitates  appears 
to  be  solely  responsible  for  the  observed  phenomenon.  At  present,  the  source 
of  boron  contamination  most  consistent  with  our  observations  is  the  result  of 
interaction  between  3,3-^  and  pBN.  A  diffusion  coefficient  of  Di'TSiOj)  <  l  x 
13"'^  cm”'  s”'  for  boron  in  InP  has  been  determined  in  this  work. 

We  have  obser/ed  cracking  in  chemical  vapor- deposited  SiO-y  encapsulating 
layers  on  InP  when  these  samples  were  annealed  at  temperatures  above  c53‘’3. 
Using  optical  microscopy,  scanning  electron  microscopy  and  Auger  electron 
spectroscopy  we  have  studied  the  detailed  nature  of  this  form  of  film  failure. 
Thermally  induced  stress  between  the  film  and  the  InP  does  not  fully  acco'unt 
for  the  phenomena  observed.  The  cracking  apparently  results  from  an 
interaction  between  mechanical  stress,  chemical  effects  and  possibly  defects 
in  the  3i3-,  layers. 


*  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (J.S. 
Army,  U.S.  Wavy,  and  U.3.  Air  Force)  under  Uontract  M333!  4-‘'^-C-3424,  and  by 
the  jffice  of  Maval  Research  under  Contract  M0301  i-Io-C-jlOo. 
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For  ion  Lt plantations  of  In?  that  ara  highly  loninatai  by  aleotronio 
stopping,  the  etaniari  lausaian  approximation  proviles  a  poor  fit  to  the 
observei  experimental  profile.  'Jsing  four  moTiants  base!  on  theoretical 
consiierations ,  >re  have  construoteh  profiles  using  a  terminatel  Figevorth  and 
Pearson  I  distribution.  It  was  found  that  neither  the  Fdgeworth  nor  Pearson  I 
provides  the  accuracy  that  would  be  expected  from  a  four  moment  distribution. 
Phe  moments  of  the  experimental  profiles  have  been  empirically  determined. 
These  moments  together  with  the  Pearson  IV  distribution  ^ere  shown  to  fit 
experimental  profiles  over  at  least  two  orders  of  magnitude. 

Vith  the  background  described  above,  we  are  currently  studying  the 
effects  of  3e,  Mg,  and  Si  implants  in  InP.  These  studies  include  impurity 
profiling,  electrical  activation,  and 
impurities . 


interaction  with  defects  and  background 
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5 . 1  Semiconductor  Crystal  C rowth  from  the  7apor  Phase: 
Ion-3urface,  Plasma,  and  Laser  Stimulated  Reactions 


We  are  studying  energetic  particle-surface  interactions  and  stimulated 

gas  phase  reactions  '(rtiich  control  the  nucleation  and  growth  kinetics, 
chemistry,  and  physical  properties  of  compound  and  alloy  semiconductors  gro'^n 
from  the  vapor  phase  by  UHV  ion  beam  sputtering,  plasma-assisted  chemical 
vapor  deposition,  and  laser- induced  chemical  vapor  deposition.  The  common 
feature  of  these  techniques  is  that  crystal  growth  proceeds  under 
non- equilibrium  thermodynamic  conditions  through  the  production  of  highly 
reactive  gas  phase  species:  excited  atoms,  metastables,  radicals,  and  ions. 
Such  species  transfer  energy  to  the  growth  surface  upon  condensation  thereby 
altering  the  surface  reactivity  as  well  as  adsorption  and  adatom  diffusion 
kinetics  allowing  film  growth  at  lower  temperatures,  a  wider  range  in 
controlling  doping  concentrations  and  tailoring  film  properties,  and  the 
growth  of  unique  metastable  materials.  This  work  is  being  pursued  from  both 
an  analytical  and  an  experimental  point  of  view  to  establish  a  detailed 
understanding  of  fundamental  film  growth  mechanisms.  We  have  recently 
published  two  invited  review  papers  ^5.1,  5.5]  in  this  area.  Results  from 
this  research  have  a  wide  range  of  applications  in  addition  to  crystal  growth 
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inoluiin?  the  aetiYe  research  fiells  of  reacti'/e  ion  etohin?,  ion  bean 
lithography,  nicrochenical  analysis,  plasna  chenistry,  anl  1 aser-nacerial a 
interactions . 

0. ' .  1  Irowth  of  Hi'^h  Quality  3ingle  1  rystal  la-^s  ? i Ins 
by  3putter  Deposition 

Epitaxial  3aAs  filns  have  been  gro-;vn  on  (lOO)  3r-iopei  laAs  substrates  by 
snutter  ieposition  '^3.2^  from  an  -jnloped  laAs  target.  'Jnintentionally  iopei 
samples  were  high  resistivity  105  to  1 3^  (l-cm  n-type  -with  room  temperature 
mobilities  as  high  as  5000  cm- /7-sec.  This  represents  one  of  the  highest 
reported  mobilities  in  high  resistivity  3aAs.  Residual  oxygen  and  oarbon 
contamination  was  minimised  through  the  use  of  a  liquid  nitrogen  cooled  shroud 
surrounding  the  discharge  luring  deposition.  Txcess  arsenic  was  provided  to 
the  gro’nng  film  from  an  effusion  cell.  Film  doping  -was  accomplished  using 
either  an  evaporated  source  {'M)  or  the  addition  of  a  gas  phase  impurity  \H-j3'i 
during  soutter  deposition.  ^dn-doped  p-type  films  were  grown  with  room 
temperature  carrier  concentrations  between  1.3  x  1  d' and  2.3  x  1  o' -  cm"'’  and 
corresponding  hole  mobilities  of  24-0  to  55  cm-7/_3oc.  These  values  agree  very 
well  -.rLth  results  for  both  LPF  and  HBF  films  doped  to  similar  concentrations. 
The  3  iopei  films  were  n-type  (see  section  3.1.2').  In  all  oases,  the  dopant 
incorporation  probability  was  controllable  by  varying  the  negative  bias  on  the 
gro'.wing  film. 

3.1.2  Ion  Bombardment  Effects  on  Elemental  3ticking  Probabilities 

He  have  shorwn  previously  that  low  energy  (<  %'J)  ion  bombardment  of  a 

gro'wing  film  can  play  a  large  role  in  controlling  elemental  stic'icing 


probabilities,  allo'wing  an 

extension 

in 

the  growth 

temperature 

r^nge  over 

which  stoichiometric  III-7 

compounds 

and 

alloys  can 

be  formed 

as  veil 

allo’wing  the  gro’wth  of  non- 

equilibrium 

phas 

es  (see  sec 

tion  3 . 2 1  . 

“^is  odours 

through  a  variety  of  mechanisms  including  self- ion  trapping,  enhanced  surface 
reactivity,  sputter  removal  of  adatoms,  and  imoact-stimulated  dissociative 
chemisorption.  For  a  discussion  of  these  mechanisms  anl  related  equations  see 
the  review  article  by  Ireene  '^3.vl, 
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Sulfur  concentrations  of  up  to  1.2  x  1 02''  cm"^ ,  as  measured  by  secondary 
ion  mass  spectroscopy,  were  obtained  in  sputter  deposited  single  crystal  CaAs 
layers  using  an  H23  glow  discharge  source.  The  incorporated  3  concentration 
increased  linearly  with  increasing  H^s  partial  pressure  but  was  not  affected 
by  varying  the  As/Ca  flux  ratio  between  12  and  4-0.  The  S  incorporation 
probability  varied  from  0.02  for  essentially  thermal  species  to  nearly  unity 
for  sulfur  containing  ions  accelerated  by  a  negative  substrate  potential  of 
325  V.  The  enhancement  in  S  incorporation  probability  was  not  due  only  to 
implantation  processes,  but  was  also  related  to  the  strength  of  the  Ca-3  bond. 
'Jery  sharp  doping  profiles  with  no  indication  of  surface  segregation  were 
obtained  for  programmed  doping  steps  (for  comparison,  see  results  for  Sn 
segregation  in  MBS  CaAs,  from  a  study  carried  out  in  collaboration  'with 
prof.  Xorkop,  summarized  in  Unit  I). 

5 . 2  1 rowth  of  Single  C rystal  Metastable  Semiconductors* 

As  reported  last  year  r2],  we  have  carried  out  the  first  detailed 

study  of  the  growth  of  new  single  crystal  metastable  semiconductors.  The  key 
feature  in  stabilizing  the  growth  of  these  materials  ’/fas  the  controlled  use  of 
low  energy  ion  bombardment  during  deposition  to  modify  elemental  sticking 
probabilities  and  to  enhance  adatom  diffusivities .  During  the  past  year,  we 
have  investigated  the  crystal  growth,  thermal  and  temporal  stability, 
electrical  properties,  and  the  thermodynamics  of  a  new  subclass  of  metastable 
semiconductors,  (III-V)^ 

Spitaxial  metastable  ( laAs)  ^  alloys  with  compositions  between  x  = 

D.10  and  D.''’5  have  been  grown  on  (IOO)  laAs  substrates  by  ion  beam  sputtering 
in  an  ultrahigh  vacuum  system.  Electron  channeling,  double  crystal  X-ray 
diffrac tometry ,  and  X-ray  topography  analyses  indi'^r^'e  that  the  films  are  of 
very  high  crystalline  perfection.  layers  with  x-values  up  to  at  least  x  = 
D.47  were  in  the  zinc  blende  structure.  Either  n-type  or  p-type  conduction 
•with  n  and  p  varying  over  several  orders  of  magnitude  could  be  obtained  by 
varying  the  film  composition,  the  growth  temperature,  and  the  As  overpressure 
during  deposition.  The  equilibrium  laAs-Ue  pseudobinary  phase  diagram  has 

*  This  research  was  supported  by  the  Joint  Services  Electronics  Program 
(U.S.  Army,  U.S.  Navy,  ancf  U.S.  Air  Force)  under  Contract  N0001 4-'^9-C -0424-, 
and  by  the  Department  of  Energy. 
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been  ietenined. 

Recent  themolyna'nic  calculations  combinei  •;fith  differential  scanning 
calori:netry  cieasureTients  have  provided  a  uore  detailed  understanding  of  both 
the  growth  and  the  thermal  stability  of  these  unique  naterials.  RXAFS 

cieasurenients ,  carried  out  in  collaboration  -/rith  Prof.  Then  of  the  Hetallurgy 
OepartTient  at  'JI’JT.  conbined  -with  STSH  and  thermal  annealing  results  have 

sho’im  conclusively  that  collisional  nixing  due  to  low  energy  ion  bonbardnent 
can  be  used  to  retard  phase  transformations.  finally,  we  have  beg-un 
collaborations  -.nth  Prof.  Xlein  of  the  Physics  Department  at  'JINT  and 

Prof.  Raccah  of  the  Physics  Department  at  HITT  to  investigate  optical 
properties  oi  metastaole  f  TaAs'' ^  ^nd  'TaSb')  i  _xTe^.  alloys  using  Raman 

soectrosGopy ,  electro- reflectance ,  and  photolaminescence . 

5 . ~  Reactive  Ion  Itching  of  1 a\ s* 

The  first  detailed  study  of  tne  reactive  ion  etching  of  TaAs  ^o.?j  '^as 
carried  out  by  our  group  during  the  past  year.  Tf  primary  interest  was  the 
effect  of  ion-surface  interactions  on  the  mechanisms  and  icinetics  of  surface 
reactions.  Rxperiments  were  carried  out  in  pure  TTl,p,  ^nd  T?j_  discharges,  as 
w^ll  as  in  dilute  TCl^,  TClo?o»  and  T?^  discharges  containing  RD  molel  Ar. 
Anisotropic  etching  -with  removal  rates,  R,  of  up  to  D.3  ym/min  have  been 

obtained  in  reactive  discharges  operated  at  lO  mTorr  and  -d  <V,  whereas  the 
physical  sputtering  rate  in  pure  Ar  discharges  operated  under  the  same 
conditions  was  only  0.  Ti  ym/min.  A  combination  of  optical  emission  and 
absorption  spectroscopies  have  been  used  to  show  that  in  both  pure  and  dilute 
chlorine-containing  reactive  iischarges,  physical  sputtering  of  atomic  Ta  and 
As  is  not  the  primary  etching  mechanism  for  TaAs  although  ion  bombardment  does 
play  an  important  role  in  the  overall  process.  Rather,  etching  occurs  by  the 
formation  and  desorption  of  volatile  Ta  and  As  chlorides.  In  reactive 
fluorine-containing  discharges,  the  rate  limiting  step  is  the  removal  of 
non-volatile  Ta?^. 


*  This  research  was  supported  by  the  Joint  Services  T1  '  tronics  Program 
'■J.3.  Army.  H.3.  Navy,  and  U.3.  Air  Force)  under  Tontrac'  .iITTi  l-"’o-T-0i2i, 
and  by  the  I3H  Researcn  Laboratory,  East  Fishkill,  New  York. 
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5.4  Laser-Surface  Interactions 

In  collaboration  with  Prof.  Eden  of  the  Electrical  Engineering  Department 
at  UIUC  we  have  initiated  a  program  to  investigate  film  growth  by  laser 
induced  chemical  vapor  deposition  (LCVD).  In  initial  experiments  we  have 
studied  the  growth  of  amorphous  hydrogenated  Si  and  3e  thin  films,  1  cm2  in 
area,  by  photolysis.  Dilute  SiH^,  (CH-5)^Si,  or  CeH^  mixtures  in  an  inert 
carrier  gas  were  photodissociated  using  ArP  (193  nm)  or  KrP  (248  nm)  lasers. 
Film  thiclcnesses  of  up  to  0.5  pm  were  obtained  on  3i02  substrates  at 
deposition  rates  ranging  from  4  to  10  A/sec  for  active  gas  partial  pressures 
between  5  and  30  Torr  and  incident  laser  intensities  up  to  25  H'/f/cm-. 
Absorption  coefficients  of  as-deposited  Oe:H  and  Si:H  layers  were  2.3  x  1 03 
cm"'  and  5.9  x  10“^  cm"',  respectively,  and  sheet  resistivities  >_  10^  Ji/o  have 
been  measured.  Incubation  times  and  initial  growth  rates  have  been  determined 
in-situ  as  a  function  of  UV  photon  flux  from  He-Ne  laser  absorption 
measurements.  A1  concentrations  of  up  to  several  atomic  percent  have  been 
incorporated  into  Ge  films  by  the  simultaneous  photodissociation  of  GeH^  and 

(:h^)3A1. 

A  parallel  effort,  carried  out  in  collaboration  with  the  Physics  Division 
of  Eastman  SCodaic  Research  Laboratories,  Rochester,  New  York,  has  been 
established  to  investigate  liquid  phase  regrowth  of  semiconductors  by  7/1  laser 
annealing.  >fe  have  recently  demonstrated  the  growth  of  epitaxial  Ge/GaAs 
heterostructures  by  scanned  CW  argon  laser  annealing  of  440  nm  thick  amorphous 
Ge  films  on  (lOO)  semi- insulating  GaAs  substrates  [5.3].  Depending  on  the 
incident  laser  power  and  scan  rate,  two  modes  of  film  regrowth  ware  obser/ed. 
At  low  powers  (between  -v  1.6  and  4.0  W  for  a  beam  diameter  of  40  pm)  and 
scan  rates  between  1  and  400  cra/s,  polycrystalline  Ge  with  a  (lOO)  preferred 
orientation  was  formed  by  an  "explosive"  crystallization  mechanism.  At  higher 
powers,  and  over  a  scan  rate  range  of  20  to  400  cm/s,  single-crystal  films 
containing  some  dissolved  GaAs  in  solution  were  obtained  by  liquid-phase 
regrowth.  Typical  film  resistivities  p  were  as  follows;  as-deposited,  130 
G-cm;  polycrystalline  films,  3  x  10"-  fl-cm;  single-crystal  films,  9  x  10""^ 
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6. 1  Introduction 

The  objective  of  the  microwave  acoustic  research  is  five  fold: 

1 .  to  use  the  interactions  between  acoustic  waves  and  free  carriers  in 
OaAs  to  make  charge  transport  devices, 

2.  to  use  traveling  wave  electric  fields  similar  to  those  produced  by 
acoustic  waves  to  evaluate  semiconductor  electronic  properties, 

3.  to  integrate  semiconductor  active  devices  and  microwave  acoustic 
devices , 

4.  to  analyze  surface  acoustic  wave  transducers  and  reflectors  from 
basic  material  parameters,  and 

5.  to  investigate  mode  conversion  techniques  of  line  acoustic  wave 
generation. 


S.  Mahon 
H.  Miller 
S.  Vilkus 


Progress  toward  this  objective  is  described  in  the  following  paragraphs. 
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5 . 2  Electron  Transport  in  Pisgoelectric  Semiconductors* 

Acoustic  surface  and  line  '/raves  propagating  in  piezoelectric 

semiconductors  drag  the  carriers  and  produce  an  acoustoelectric  current,  'rfhen 
the  electric  fields  nroduced  by  the  waves  are  greater  than  the  sound  wave 
velocity  divided  by  the  mobility,  the  carriers  are  dragged  at  the  sound 
velocity.  Fiarther  increases  in  acoustic  power  do  not  increase  velocity  or  the 
acoustoelectric  current.  The  carriers  are  then  said  to  propagate 
synchronously  and  the  acoustoelectric  current  is  saturated. 

The  charge  transport  devices  envisioned  in  this  research  depend 
exclusively  on  synchronous  carrier  transport.  Extensive  study  has  been 
carried  out  to  determine  the  nature  of  charge  transport  in  the  presence  of 
very  large  traveling  wave  fields.  The  effect  of  traps  on  synchronous  carrier 
transport  is  of  particular  concern  and  Cr  doped  ^aAs  has  been  studied 

extensively  because  of  its  high  density  of  traps. 

The  saturated  acoustoelectric  current,  caused  by  the  nonlinear 
interaction  of  photo-generated  electrons  'With  large  amplitude  SA'rf  induced 
electric  fields  in  Cr  doped  CaAs,  is  fo'jnd  to  be  too  small  to  be  explained  on 
the  premise  that  all  electrons  move  synchronously  at  the  SAW  velocity.  This 
finding  is  based  on  the  results  of  large  signal  acoustoelec trie  measurements, 
for  SAW  both  internal  and  external  to  the  CaAs  substrate,  -/-ia  a  separated 
media  device  geometry.  The  expected  linear  dependence  of  saturated  current 
•with  SAW  velocity  is  not  observed  when  SAW  with  different  velocities  are  used 
to  induce  the  saturation.  It  is  concluded  that  some  mechanism  other  than 
complete  carrier  bunching  is  responsible  for  acoustoalectric  current 

saturation  in  CrrSaAs  [r1 ]. 

During  the  last  year,  we  have  demonstrated  that  electrons  can  be 

entrained  by  an  acoustic  wavs  propagating  along  the  surface  of  a  Cr-doped  CaAs 
crystal.  This  result  'was  surprising  in  view  of  the  weak  piezoelectricity  of 
SaAs;  it  was  made  possible  by  the  unusual  property  of  high-frequency  acoustic 
surface  -waves  of  providing  a  concentrated  field  over  a  narrow  cross  section. 


Electronics  Program  under 


*  This  -work  was  sunoorted  by  the  Joint  Services 
conract  NOOOU-^9-T-D‘124.. 
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6.  MICROWAVE  ACOUSTICS 

One  can  visualize  a  new  class  of  high-performance  charge  transfer  devices 
based  on  acoustically  entrained  carriers.  Such  devices  would  have  no  gates 
and  no  complex  multi-layer  patterns  to  interconnect  the  gates;  the  region 
along  »taich  the  electrons  travel  would  be  free  for  pickup  electrodes.  Along 
with  the  simpler  structure,  these  devices  would  have  negligible  transfer  loss, 
making  it  possible  to  obtain  a  time- band wid th  product  orders  of  magnitude 
higher  than  that  of  present-day  CCD's. 

5.3^  New  Method  to  Measure  the  Mobility  of  Amorphous  Semiconductor  Films* 

The  measurement  of  mobility  in  high  resistivity  low  mobility  amorphous 
semiconductor  films  can  be  difficult.  We  have  used  the  acoustoelectric 
current  described  in  the  last  section  as  a  measure  of  material  mobility. 

The  acoustoelectric  method  of  determining  carrier  mobility  u  in 
semiconductor  films  has  been  modified  to  permit  its  application  to 
high-resistance,  low-mobility  films.  The  conventional  method  requires 
knowledge  of  the  acoustic  loss  caused  by  the  mobile  carriers;  in 
high-resistance  films  this  loss  becomes  too  small  to  be  measured.  We  show 
that  the  required  information  may  be  derived  from  knowledge  of  the  acoustic 
power  and  the  device  geometry.  Our  samples  were  amorphous  hydrogenated  St  and 
^^0.6^®0.4  fil®3  on  nonpiezoelectric  substrates,  separated  by  a  convenient  air 
gap  (12.5  pm)  from  a  LiNbO-^  slab  carrying  surface  acoustic  waves.  One  sample 
had  1 0®  /□  and  u  =  0.08  cm-/Y  sec;  another  sample  52/n  and  u  =  0.5  cm~/V 

sec  [0.1]. 

In  amorphous  materials  the  drift  mobility  of  carriers  is  usually  measured 
by  observing  the  time-of- flight  of  a  charge  packet  through  a  sample  in  an 
external  d.c.  electric  field.  The  results,  however,  cannot  be  interpreted  in 
terms  of  a  carrier  mobility  in  the  conventional  sense  due  to  the  spreading  of 
the  packet  by  multiple  trapping  and  release,  and  the  concept  of  time-dependent 
mobility  has  been  used  to  describe  such  dispersive  transport.  Recently  the 
mobility  was  also  deduced  from  transient  photoconductivity  and  photo-induced 
absorption  experiments.  Using  the  acoustic-electric  method  to  measure 
mobility  provides  new  information.  A  simple  coherent  view  showing  the 
connection  among  all  these  different  measurement  techniques  has  been  written 

•  This  work  was  supoortei  by  the  Joint  Services  Electronics  Program  under 
contract  N0001  4-79-1-1424. 
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6 . 4  'Wideband  LAW  Transducer* 

The  fundaTiental  advantages  of  LAW  devices  are  veil  documented  but  there 
has  been  no  efficient  viiebandwilth  transducer  vith  which  to  in\plement  them. 
We  have  recently  devised  a  mode  converter  that  converts  surface  acoustic  wave 
SAW  power  to  LAW  power.  This  device  consists  of  a  number  of  unconnected 
electro! S3  which  extend  across  the  SAW  and  LAW  propagation  paths.  The 

electrodes  are  bent  and  spaced  so  that  the  electrode  spacing,  measured  in 

phase  at  center  frequency,  is  the  same  for  the  LAW  and  the  SAW  beams.  The 

first  generation  device  had  a  5^  bandwidth  and  a  loss  of  5  13 

The  second  generation  has  a  bandwidth  increased  to  1  Ot  and  a  loss  reduced 
to  1  13.  With  this  band'width  and  efficiency,  LAW  devices  may  become  practical 
because  efficient  SAW  transducers  can  be  used  to  generate  the  SAW  beam  which 
in  turn  is  converted  to  the  LAW  beam  by  the  Multistrip  Mode  Converter  (MS'dC). 
The  mode  converter  is  reciprocal  and  LAW  waves  can  be  detected  by  conversion 
to  SAW,  and  then  a  SAW  transducer  is  used  to  detect  them,  slfficient  moderate 
bandwidth  LAW  delay  lines  and  convolvers  are  now  possible. 

6.5^  Theoretical  Tetermination  of  Surface  Acoustic  Wave  Velocity 

and  Impedance  Differences  between  Metal  Strips  and  Tree  Surface 
Regions  of  Metallic  S ratings** 

The  surface  acoustic  wave  is  a  single-moied ,  con- iistiersive , 
one- dimensional  wave  that  is  frequently  modeled  using  a  transmission  line 
equivalent  circuit.  Thin  electrodes  and  shallow  grooves  located  in  the 
surface  acoustic  wave  propagation  path  produce  reflections  and  shift  the 
propagation  velocity  of  the  wave.  The  effect  of  the  electrodes  and  grooves 
are  modeled  as  an  offset  in  characteristic  impedance  'DRLT3RT')  and  a  shift  in 
the  propagation  velocity  (DSLTSHV),  This  impedance  offset  and  the  velocity 
shift  have  been  theoretically  detsrminel  from  the  material  constants  and  the 

*  This  work  was  suonorte!  by  the  Rome  Air  Development  Tenter,  Deouty  for 
Electronic  Technology',  under  Air  Force  contract  ?-l  hoZS-31 . 

**  This  work  was  suooortei  by  the  Rome  Air  Develooment  Tsntcr.  Deouty  for 
Electronic  Technology,"  under  Air  Force  contract  F-l  loZ T-IDID. 
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geonietrical  dimensions.  Since  the  materials  are  piezoelectric,  both  the 
elastic  and  electrical  loading  of  the  surface  have  been  considered.  The 
detailed  theoretical  derivations  for  the  impedance  offsets  and  velocity  shifts 
have  been  published  in  a  series  of  papers  during  the  course  of  the  development 
[5.2, 6. 4-6].  A  report  has  been  written  to  tie  these  studies  together  and  to 
provide  a  simplified  analysis  of  the  surface  loading  produced  by  an  electrode 
that  will  be  useful  to  the  device  designer.  The  report  has  been  made  concise, 
with  all  derivations  taken  from  references  or  appendices  in  order  to  make  it  a 
better  design  tool. 
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7 . 1  Dif fuaion  and  Seaoti'/ity  <2^  Metals'* 

7h3  ienompoaition  of  polyatooiio  ^asea  on  aoliia  is  of  interest  in  both 
oheTiical  ani  eleotronica  technology.  Me  have  been  exaiining  the  interactions 
of  nethane  vith  letala  as  a  simple  but  manageable  example  of  each  processes. 
This  effort  has  bean  concentratel  upon  'unierataniing  the  mechanisms  by  ehieh 
lecomposition  occurs;  of  particular  interest  is  -/ay  there  is  an  activation 
energy  to  iecompositioa  at  the  surface  and  how  this  is  oyercome. 

Measurements  reportei  last  year  '”,31  1  indicatei  that  the  activation  energy, 
leterminel  when  the  temperature  of  both  the  surface  and  the  gases  is  changed, 
is  considerably  smaller  than  the  value  fo'and  some  time  ago  32  i  in  beam 
experiments,  in  which  only  the  gas  is  excited.  7.ae’'e  are  two  dietinct 
possibilities  to  acccunt  for  this: 


^Ihis  worx  was  supported  bv  the  'lational  3cience  Moundation  inder,  grant  hM3 
~3-2oll  -3.  as  well  'as  by  tne  loint  3ervices  Mlectronics  ^'•ogram  '‘J.l.  l-my, 
■J.3.  Mavy,  ’1.3.  \ir  7orce)  inder  contract  -niO!  d-'"}-: -Ul  1. 
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1.  The  beam  experiments  were  done  on  single  erystal  tips,  rather  than  on 
evaporated  films,  whioh  served  as  substrate  in  the  more  recent  worx. 
There  may  be  a  significant  dependence  of  decomposition  rate  upon 
crystallography,  which  could  account  for  the  differences. 

2.  Tf  the  reaction  proceeds  through  a  precursor,  that  is  through  a  gas 
molecularly  bound  to  the  surface  prior  to  dissociation,  different 
activation  energies  are  expected  in  the  two  measurements. 

3oth  types  of  effects  are  likely,  and  we  have  carried  out  several  studies 
to  define  mors  closely  the  important  events  at  the  surface.  Measurements  it 
low  temperatures  indicate  a  heat  of  physical  adsorption  of  5  kcal/mole.  Under 
the  conditions  of  the  decomposition  experiments,  a  reser/oir  of  molecularly 
held  methane  therefore  exists  on  the  surface  and  dissociation  should  take 
place  via  these  molecules.  This  suggests  that  in  a  molecular  beam  experiment 
in  which  the  gas  alone  is  heated,  the  activation  energy  should  be  c  kcal/mole 
higher  than  in  stardard  measurements,  and  should  amount  to  roughly  11 
kcal/mole  on  rhodium  film.  To  test  this  hypothesis  a  new  apparatus,  shewn  in 
Pig.  ■i'.  1  ,  has  been  built.  In  this  equipment,  the  gas  can  be  heated  while 
keeping  the  film,  which  is  in  contact  with  a  large  thermal  mass,  at  a  constant 
temperature.  The  results  obtained  in  this  way  are  plotted  in  Pig.  i'.2. 
After  correcting  for  contributions  from  gas  molecules  that  thermally 
accommodate  to  the  surface,  an  activation  energy  of  10.7  kcal/mole  is  found 
for  methane  decomposing  on  rhodium  films. 

This  is  in  excellent  agreement  with  expectations;  more  than  that,  it  ie 
consistent  .with  other  phenomena  observed  in  the  interactions  of  methane  with 
rhodium.  It  has  previously  been  fo'und  [pM.Rd]  that  exciting  the  v  stretching 
mode,  which  has  an  energy  of  3.5  kcal/mole,  is  ineffective  in  accellerating 
decomposition.  Tnis  is  now  luite  understandable .  An  energy  of  11  kcal/mole, 
much  more  than  is  available  in  the  vibrations,  must  be  supplied  to  methane 
in  order  for  it  to  dissociate.  Most  likely  it  is  the  bending  mode  that 
must  be  excited  in  order  to  enhance  the  reacti-cn  on  rhodium.  Detailed 
predictions  of  the  effect  of  isotopic  substitution  have  been  made,  based  on 
the  assumption  that  vibrational  excitation  is  important;  these  predictions 
are  in  good  accord  with  experiments. 


4  Wire  Electrical 
Feed  Through 


rhodium  films  maintained  at  O^C. 
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The  wor'rc  cirriai  oit  here  has  elarifiel  the  role  of  the  ^as  in  the 
ieoo-npeaitioa  reaction.  Still  to  be  explored  now  are  the  contributions  of  the 
eolii,  and  the  extent  to  which  these  are  speoifioally  affected  by  the  atonic 
arran^exent  at  the  surface. 

7 . 2  diffusion  and  Reactivity  on  Silicon* 

Although  the  ixportance  of  surface  effects  in  solid  state  devices  is 
generally  recognised,  relatively  little  is  'xnown  about  the  properties  of 
surface  layers  on  sexiconductors.  '.is  have  been  concentrating  our  efforts  on 
'Obtaining  inforxation  on  the  spatial  stability  of  such  layers.  As  part  of 
these  invest igations  we  have  developsd  a  techniq^ue  for  direct  electron  beax 
writing  on  silicon  In  this  xetho d ,  xolecular  nitrogen  is  physically 

adsorbed  on  silicon,  xaintained  at  low  taxperatures  at  which  a  xonolayer  or 
xore  is  forxad.  Molecular  nitrogen  on  the  surface  is  then  dissociated  by 
electron  beax  Lxpact.  In  warxing  the  saxple,  undissociated  xolecular  nitrogen 
evaporates,  leaving  behind  a  nitrogen  layer  wherever  the  electron  beax  has 
struc’x,  as  3ho'.^n  in  7ig.  "’.M* 

Current  efforts  have  been  focused  upon  charao tericing  the  conditions 
uadar  which  these  nitrogen  layers  are  forxsi.  Cur  studies  indicate  that 
electron  ixpact  is  quite  efficient  in  proxoting  dissociation  of  nitrogen, 
less  than  5  x  I  C'  '  electrons/cx-  suffice  to  create  a  saturation  laye"; 
further  bcxbardxent  does  not  increase  the  axo'unt  of  nitrogen  chexically  bound 
to  the  silicon. 

These  saturation  layers  are  actually  nitrides,  not  just  chexisorbed 
layers  of  nitrogen.  It  is  clear  from  Fig.  h.q  that  the  3i'L7V)  Auger  spectra 
of  the  layers  produced  by  elecron  boxbardxent  are  in  good  accord  with  the 
spectra  recorded  for  silicon  nitride  produced  by  chexicel  vapcr  deposition 
[do].  Me  have  been  able  to  show,  however,  that  during  the  early  stages  of  the 
process,  what  appears  to  be  a  chexisorbed  layer  of  nitrogen  does  forx.  This 
is  illustrated  in  Fig.  ''•5,  where  the  intensities  of  various  Auger  lines  are 
shown  after  increasing  exposures  of  nitrogen  and  concoxitant  boxbardxent  ^ith 
electrons.  The  dd-ay  line  charac teristi c  of  the  nitride  does  not  appear  at 

*  This  worx  was  supported  by  the  Toiat  der^ices  Tleotronics  '^■’Ograx 
'■J.d.  Arxy.  T.d.  'lavy, '  O.  d.  Air  Force)  under  contract  'lOO'  l-’d-d-Clll. 
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Fig.  7.3 


Secondary  electron  Images  of  nitride  structures  created  by  electron 
beam  writing  on  molecular  nitrogen  layers  held  on  Si(lOO). 
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Fig.  7.5  Auger  Intensities  observed  after  increasing  nitrogen  exposures 
and  electron  irradiation.  The  83-eV  nitride  Line  only  develops 
after  an  induction  period,  during  which  a  chemisorbed  layer  forms 
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very  start.  However,  the  330-eV  line  attribute!  to  nitro'?en  is  apparent 
even  at  /ery  low  exposures,  iniieating  that  there  is  a  nitrogen  layer  present, 
though  it  is  not  a  nitrile.  The  nitride  layers  formed  after  extensive 

exposure  are  more  than  one  Layer  thiox.  The  nechanism  of  formation  of  these 
nitriles  is  not  clear,  but  it  is  quite  complex  and  may  involve  a  competition 
between  electron  impact  desorption  and  dissociation. 

Electron  beam  induced  formation  of  silicon  nitride  is  a  particularly 
appealing  process.  It  gives  us  the  capability  of  forming  nitrile  structures 
under  controlled  and  very  clean  conditions.  A  question  immediately  arises  - 
is  this  technique  applicable  to  the  formation  of  other  films  as  well?  Ve  have 
recently  examined  the  feasibility  of  creating  oxide  layers  on  silicon  in  this 
way.  On  exposing  silicon  to  molecular  oxygen,  a  chemisorbed  layer  forms  on 
the  surface.  At  room  temperature  and  below,  oxidation  is  alow;  however, 

electron  bombardment  enhances  the  oxidation  rate.  As  is  apparent  from  Fig. 

■^.o,  it  is  therefore  possible  to  form  oxide  structures  by  electron  beam 
irradiation  of  an  oxygen  covered  surface.  ?he  most  advantageous  conditions 
for  this  process  are  now  being  explored. 

.  5  Atomic  Exploration  of  1  rystal  Layers* 

Over  the  years  considerable  information  has  been  acc'jmulated  in  this 
laboratory  about  the  behavior  of  individual  metal  atoms  on  metal  surfaces 
This  has  been  possible  through  the  use  of  the  field  ion  microscope  and 
its  ability  to  resolve  single  atoms.  'jrnile  work  on  single  atoms  has 
continued,  we  have  also  been  examining  the  behavior  of  overlayers,  with  the 
intent  of  contributing  to  a  better  understand ing  of  crystal  and  film  growth 
phenomena;  in  the  past  these  phenomena  have  not  been  accessible  to 
observation  on  the  atomic  level. 

Initially,  we  have  concentrated  on  a  survey  of  overlayer  formation  on 
V(nO').  Metals  from  the  platinum  family  form  two-dimensional  layers,  provided 
the  adatom  concentration  on  the  substrate  is  high  enough.  As  is  apparent  in 
the  field  ion  micrographs  in  Fig.  7.7,  individual  atoms  in  the  layers  are 
clearly  resolved.  This  is  quite  different  from  -what  is  found  when  the 

*  This  work  was  3upp<crted  dyr  the  National  Science  Foundation  under  grant  OMR 
30-11061  and  by  the  Joint  cervices  Electronics  Program  'M.o.  Army,  'J.3.  Navy, 
■J.3.  Air  Force)  'under  contract  NOOO' I-’O-O-OIZI. 
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Fig.  7.6  Secondary  electron  image  of  oxide  structures  formed  by  electron 


beam  writing  on  Si(lOO). 
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Fig.  7.7  Platinum  metal  layers  on  W(llO).  Shown  in  order  are  palladium, 
iridium,  and  platinum  monolayers,  as  well  as  a  schematic  of  a 
p(2  X  2)  structure. 


Fig.  7.8  Field  ion  image  of  silicon  layer  on  (110)  plane  of  tungsten 
(left),  and  proposed  structure  at  right. 


■’.  SURFACE  STUDIES 


subatfata  ibaelf  ia  imaged.  3n  a  ’/{(l  1  D)  plane  of  the  aame  aize,  individual 
atoma  cannot  be  diacerned.  The  apacing  ia  the  layera  of  platinum,  palladium, 
and  iridi'um  muat  therefore  be  a  multiple  of  the  aubatrate  apacing.  Thia  ia 
not  neceaaarily  true  for  other  metala;  for  rhenium,  individual  adatoma  are 
not  reaolved  in  overlayera  on  V(110),  and  we  therefore  infer  a  much  ienaer 
packing . 

A  particularly  intereating  behavior  haa  been  obaerved  for  iridi^nm 
cluatara.  Tneir  arrangement  and  atability  ia  strongly  dependent  upon  aize. 
At  4-50  ’<  for  example,  a  cluster  of  16  or  more  iridium  atoma  on  W(l  1  3)  ia 
stable  in  a  rhombohedral  config’xration.  A  cluster  of  only  'c  adatoma, 
however,  occurs  in  a  square  arrangement.  Tnia  ia  not  an  artifact;  the 
cluster  maintains  ita  shape  'under  repeated  heating  intervals.  Tnen  a  single 
atom  is  removed  from  this  layer,  the  arrangement  becomes  'unstable  -  the 
cl'uater  breaks  up  into  individ'ual  iridi’um  chains. 

Thia  sort  of  behavior  ia  not  limited  to  iridium.  As  an  important  part  if 
CUT  studies  we  have  been  loo'tciag  at  silicon  layers  on  W(!13).  Moderate  ioaes 
of  ailicon  equilibrated  at  300  K  yield  a  very  open  overlayer,  aho'wn  in  Fi?. 

made  up  of  eilicon  chains  oriented  along  "".''Oh  The  apacing  between 
chains  ia  close  to  9.5  A.  Closer  packed  layera  can  be  formed,  however,  by 
equilibrating  heavier  doaea  at  higher  temperatures.  3y  field  evaporation  the 
concentration  of  silicon  adatoma  can  be  lowered  at  will.  Repeated  field 
evaporation  followed  by  equilibration  converts  these  dense  layers  into  chain 
arrays,  analogous  to  those  in  Fig.  '^.5.  On  further  reducing  the  aiatom 
population,  these  overlayera  iiaorder  similarly  to  iridium  layera. 

The  formation  of  these  looaely  structured  layers  on  W(n3)  has  not  been 
previously  obseir/ed  by  other  techniques  ^Ro].  A  survey  has  been  made  of 
silicon  layers  on  other  planes.  3o  far  nothing  unusual  haa  been  observed  on 
(310)  or  (l30).  It  ia  already  clear,  however,  that  field  ion 
microacopic  examination  is  capable  of  yielding  interesting  informaticn  about 
overlayera,  and  these  atuiiea  as  well  as  other  measurements  on  individual 
atoms  and  small  clusters  are  continuing. 


4 


40 


STUDIES 


Sefsrenoes 


]  SSL  Annual  Report  of  Progress  1979-1  93D,  7.1,2. 

[R2]  S.  U.  Stewart  and  S.  Ehrlioh,  J.  She®.  Fnys.  4572  (l97S). 

[93]  S.  0.  Brass,  D.  A.  Reed,  and  0.  Ehrlich,  J.  She®.  Phys.  21>  5244  (l979). 

[94]  J.  T.  Yates,  Jr.,  I.  J.  Zinck,  3.  She and ,  and  W.  H.  Ueinberg,  J.  She®. 
Phys.  22,  2256  (l  979). 

i.R5]  R.  Bezel  and  U.  Lieske,  J.  Appl.  Phys.  oL*  (1980). 

[r6]  j.  P.  Siberian  and  S.  A.  3o®orjai,  J.  Vac.  3ci.  Technol.  rS.. 

(1  979). 


3. 


ULSCTROMAGNETTC  COMMUNICATION,  RADIATION  AND  SCATTERING 


Faculty 


R.  Mittra 


Graduate  Students 


M.  Desai  R.  latnpe  3.  Ray 

R.  Kastner  J.  Maguire  T.  Trinh 

R.  Pal  eta 

3.1  Millimeter- Vave  Integrated  Circuits 

In  the  past  twelve-month  period,  we  have  made  significant  progress  In  the 
development  of  devices  and  structures  for  use  in  millimeter-wave  integrated 
circuits.  A  number  of  analytical  techniques,  applicable  to  the  dielectric 
waveguides  used  in  these  circuits,  have  also  been  developed. 

A  particular  area  of  interest  has  been  that  of  dielectric  antennas.  We 
have  investigated  a  type  of  leaky-wave  antenna  which  ccnsista  of  a  dielectric 
waveguide  on  which  are  placed  periodic  metallic  perturbations.  This  antenna 
is  frequency  scannable  and  is  readily  integrated  with  a  receiver- transmitter 
system.  Prior  to  cur  investigation,  it  was  believed  that  the  perturbations 
had  to  be  'uniform.  Our  experiments  have  shewn  that  nenuniform  perturbations 
result  in  significantly  better  antennas.  Antennas  -with  26  dB  sunpressen  of 
sidelcbe  and  endfire  radiation  levels  have  been  built.  These  antennas  can 
have  beamwidths  in  the  plane  parallel  to  the  antenna  axis  as  low  as  '^.0'-  . 


*  This  work  was  supnorted  by  the  -Joint  cervices  Electronics  Program 
'■J.3.  Army,  U.3.  Navy,  'J.3.  Air  rcrce''  'under  contract  NOOQi  i-"’0-0-012  i. 


3.  SLECTROMAaUETIC  COMUNICATION,  RAQIATTON  AJTO  SCATTERIJI:! 


For  some  applications,  the  radiation  pattern  of  this  antenna  in  the  plane 
transverse  to  the  guile  axis  may  be  too  broad.  To  reduce  this  beamwilth,  a 
new  structure  was  developed  in  which  the  dielectric  guile  with  periodic 
perturbations  on  top  is  embedded  in  a  metal  trough  with  a  metal  flare  attached 
on  each  side,  'jfith  this  arrangement,  the  antenna  behaves  like  a  linear  array 
in  the  longitudinal  plane,  while  the  radiation  pattern  in  the  transverse  plane 
resembles  that  of  a  horn.  'Vith  this  structure  we  can  achieve  7°  beamwiiths  in 
the  transverse  plane. 

!^any  passive  components  in  a  dielectric-based  millimeter- wave  integrated 
circuit,  such  as  ring  resonators  and  couplers,  are  comprised  of  curved 
waveguide  sections.  Therefore,  it  is  necessary  to  have  a  good  understanding 
of  the  behavior  of  the  fiali  in  the  neighborhood  of  the  bend  in  order  to 
predict  the  coupling  and  radiation  characteristics.  General  analytical 
approaches  to  this  problem  have  been  very  involved  to  data.  We  have  found  a 
simple,  approximate  technique  for  calculating  tha  field  ^n  a  curved  dielectric 
waveguide  of  rectangular  cross-section. 

Our  analysis  is  based  on  a  conformal  mapping  which  is  used  to  transform  a 
homogeneous  curved  waveguide  into  a  straight  waveguide  with  an  inhomogeneous 
refractive  index.  Solution  of  the  resulting  equations  yields  Airy  functions 
which  accurately  describe  the  field  both  inside  and  outside  the  waveguide. 

In  the  area  of  active  device  design  we  have  worked  on  the  problem  of 

I'iPATT  oscillators  which  are  directly  integrable  with  dielectric  waveguides. 
'Jsually  this  results  in  excessive  mismatch  and  radiation  loss.  Through  the 
use  of  partial  shielding  structures  and  better  matching  networks  we  have 
reduced  these  losses  significantly  and  realized  oscillators  which  work  at  34- 

Tdz. 

A  new  method  of  solving  dielectric  waveguide  discontinuity  problems  is 
also  being  developed.  Unlike  a  closed  waveguide,  dielectric  waveguides 

support  radiation  modes  with  a  continuous  eigenvalue  spectrum.  This  makes  the 
application  of  mode-matching  and  related  techniques  very  difficult.  lur 

method  is  based  on  the  spectral  domain  method  and  utilises  the  F7T  algorithm 
rather  than  matrix  inversion.  At  present  small  steps  on  dielectric  waveguides 
have  been  handled  successfully.  We  are  currently  incorporating  a  variational 
expression  to  enhance  convergence  for  larger  discontinuities . 
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3.2  Electromagnetic  Radiation  and  Scattering 

This  area  of  our  research  is  primarily  concerned  with  the  efficient  and 
accurate  calculation  of  electromagnetic  radiation  and  scattering  problems.  Ne 
have  concentrated  our  efforts  in  the  high  frequency  region  where  moment  method 
solutions  become  impractical  and  where,  either  by  body  size  or  complexity, 
asymptotic  solutions  are  not  applicable.  In  the  past  a  spectral  domain 
approach  has  been  formulated  and  applied  to  thin  plates.  Recently  we  have 
been  able  to  accurately  extend  this  technique  to  two-dimensional  bodies  of 
significant  size,  which  may  be  either  perfectly  conducting  or  dielectric  in 
nature.  As  this  technique  is  iterative  and  uses  the  P?T  algorithm,  it 
contains  a  built-in  boundary  condition  check  and  is  very  efficient. 

Ne  have  also  done  experimental  and  theoretical  studies  of  twist 
polarizers.  These  polarizers  can  be  used  in  radar  systems  where  traditional 
gimbaled  systems  are  too  slew,  as  the  twist  polarizer  is  very  light  weight  and 
can  be  moved  rapidly.  The  theoretical  analysis  uses  the  spectral  domain 
technique  and  gives  good  agreement  with  experimental  results. 


MdiuataKiMGiitii 
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9- t  Statiatloal  ? ropertiea  of  ? laama  Turbulence 

The  ob.jectire  of  thie  reeearoh  is  to  study  the  statistical  properties  of 
plasma  turbulence.  The  system  ’onier  investigation  is  the  positive  column  of  a 
Ic-iischar^e  in  helium  -yhich  is  •unstable  ■fith  respect  to  the  spontaneous 
e'ccitation  of  ion  acoustic  vaves.  cllectric  field  fl’uc  tuations  of  the 
turbulent  ,;ave  field  are  picked  up  by  a  Langm-uir  probe  and  their  probability 
distribution  is  measured  by  a  fast  sampling  method.  The  measured  deviations 
from  daussian  behavior  can  be  fitted  to  a  dram-Charlier  expansion  of  the 
probability  distribution  including  terms  up  to  fourth  order.  Tne  non-daissian 
character  of  the  probability  distribution  allo'<rs  conclusions  to  be  ira-ifn  as  t.^ 
the  nature  of  the  ensemble  of  the  xaves  that  compose  the  turbulent  field. 
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10.1  Rarefied  Cas  Dynaaics* 

Rarefied  gas  dynamics  deals  with  non- equilibrium  gas  flow  problems  in 
which  microscopic  treatment  according  to  ’<inetic  theory  is  necessary  to 
determine  the  effect  of  intermolecular  collisions  and  gas  surface  interactions 
on  both  microscopic  and  macroscopic  gas  flow  properties.  Such  rarefied  gas 
flow  problems  occur  not  only  in  aerodynamics,  but  also  in  electronics, 
aeronomy,  environmental  fluid  dynamics,  and  other  related  fields. 

The  aim  of  this  research  program  is  to  develop  numerical  methods  to  solve 
a  -wide  range  of  problems  wonder  conditions  far  from  and  near  thermal 
equilibrium.  A  Monte  Carlo  method  has  been  developed  at  the  Coordinated 
Science  Laboratory  i_R1  ]  to  solve  directly  the  Roltsraann  equation  and  has  been 
used  by  the  Boltzmann  group  to  solve  the  Boltzmann  equation  for  several 
rarefied  gas  flow  problems  under  a  wide  range  of  nonequilibrium  and  bo'indary 
conditions  [R2-R7].  The  solutions  we  have  obtained  yielded  detailed 
microscopic  and  macroscopic  non-equilibrium  properties,  most  of  which  have 
never  been  treated  and  studied  before.  We  have  also  studied  numerical 
solutions  of  other  kinetic  equations  and  other  numerical  methods  to  solv' 
rarefied  gas  flow  problems,  including  the  direct  simulation  techni:iue. 

*  This  work  -was  sunoorted  by  NATO  Research  Grant  1075  and  by  the  University  of 
Illinois. 
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Nonlinear  evaporation-condensation  problems  are  encountered  in  such 
diversified  areas  as  upper  atmosphere  meteorology,  the  cooling  of  nuclear 
reactors,  design  of  space  experiments,  petrochemical  engineering,  vacuam 
technology,  and  the  interaction  of  high  power  laser  radiation  with  metal 
surfaces.  The  treatment  of  these  problems  requires  first  the  consideration  of 
vapor  fcinetics  problems  characterised  by  the  nonequilbrium  vapor  motion  in  a 
•Cnuisen  layer  at  the  interphase  boundary.  Under  a  joint  research  effort,  we 
have  successfully  studied  the  evaporation  problem  [^^l.  Our  Boltzmann 
solutions  establish  the  validity  of  the  Krook  solution  and  Ytrehus'  kinetic 
theory  approach  [rs].  This  approach  can  be  used  to  calculate  the  jamp 
conditions  across  the  Knudsen  layer  and  the  net  mass,  momentum  and  heat  fluxes 
[R10]  and,  together  with  the  conventional  continuum  method,  to  calculate  the 
flow  parameters  at  the  outer  edge  of  the  Knudsen  layer  for  problems  in  which 
evaporation  occurs  at  the  interphase  boundary. 

Studies  have  also  been  made  for  the  condensation  problem  ’'_10.2J.  The 
non- equilibrium  vapor  near  a  condensing  surface  differs  from  that  near  an 
evaporating  surface  because  of  the  difference  in  edge  condition  at  the 
interpase  boundary  and  its  behavior  is  more  complex.  Our  Boltzmann  solutions 
yield  distinct  non-equilibrium  behavior  in  the  Knudsen  layer  of  condensing 
vapor.  For  example,  the  temperature  gradient  becomes  negative  near  the 
interphase  -under  certain  conditions.  The  condensation  rate  as  a  function  of 
pressure  ratio  was  found  to  depend  strongly  on  a  substance  parameter  which  is 
proportional  to  the  ratio  of  the  latent  heat  and  its  liquid  surface 
temperature.  Some  of  our  results  were  found  to  be  in  agreement  with  those 
obtained  by  using  the  simple  apnroaches  of  Oguchi  ’^Rli]  and  Ytreus  '_Rl2l. 

'tie  plan  to  study  further  the  evaporation  and  the  condensation  problems  by 
considering  more  complex  flow  and  boundary  conditions. 

13.2  Naval  Hydrodynamics* 

Free  surface  wave  problems  encountered  in  Naval  hydrodynamics  are 

characterized  by  complexities  in  flow  geometry,  flow  features  and  boundary 

conditions.  The  flow  has  an  unkno-xn  free  surface  and  it  is  propagative  an  i 

*  This  work  was  sunported  by  the  Office  of  Naval  Research  -inder  Oontract 
N0001 d-B0-C-07d0  and  the  National  Science  Foundation  under  Grant  13?  ENG 
77-204S5. 
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transient.  The  boundary  condition  at  the  free  surface  is  of  a  nixed, 
parabolic  type  and  it  contains  highly  nonlinear  terns.  In  the  steady  state, 

there  also  exists  a  radiation  boundary  condition,  since  the  waves,  once 

generated,  propagate  downstream.  These  complexities  have  led  to  several 
computational  difficulties:  accurately  accommodating  the  free  surface 

geometry,  satisfying  the  boundary  condition  uniformly  over  the  free  surface, 
treating  the  radiation  boundary  condition  and  implementing  the  boundary 
condition  at  the  contact  line  of  the  free  surface  with  the  solid  surface  of  a 
partially  submerged  body. 

In  the  initial  phase  of  our  study,  we  have  developed  two  time-dependent 
numerical  schemes  with  Tulerian  grid  systems  for  solving  steady  and  unsteady 

notential  flows  for  nonlinear  free  surface  problems.  In  one  method 

[aid,RU,Rl6],  we  use  the  finite  element  method  to  deal  '-rith  the  geometrical 
complexity  and  the  free  surface  boundary  condition  of  the  nonlinear  free 
surface  problems.  In  this  scheme,  the  finite  element  method  is  used  to  make 
the  field  calculation  and  the  finite  difference  method  is  used  for  the  time 
evolution.  We  have  used  this  method  to  solve  two  problems:  a  pressure 
distribution  moving  with  a  constant  speed  and  a  moving  submerged  elliptical 
cyliner  or  a  symmetrical  hydrofoil. 

We  have  also  developed  an  explicitly  time-dependent  finite  difference 
scheme  |_R17l.  Rxplicit  schemes  may  be  more  favorable  for  solving  large  scale 
problems  on  "supercomouters"  which  have  more  stringent  implementation 
requirements.  We  have  used  this  method  to  solve  two  problems:  a  pressure 
distribution  moving  with  a  constant  speed  and  an  elliptic  cylindrical, 
surface-piercing  strut  accelerating  from  rest.  Our  solutions  serve  to 
demonstrate  the  feasibility  of  using  the  method  to  solve  two-dimensional  as 
well  as  three-dimensional  problems. 

The  results  of  the  initial  phase  of  our  study  were  presented  at  the 
Second  International  Conference  on  Numerical  Ship  Hydrodynamics  LR’5j. 

In  both  numerical  schemes  mentioned  above,  the  computational  domain  is 
expanded  downstream  periodically  during  the  computation  as  the  disturbance  on 
the  free  surface  is  propagated  close  to  the  domstream  boundary-  Therefore, 
the  ’undisturbed  condition  is  applied  on  the  cut-off  downstream  bo’unlary. 
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In  order  to  obtain  the  longer  time  and  steady  state  solutions,  we  have  to 
use  a  fixed  downstream  boundary  set  close  to  the  disturbance  so  that  we  can 
increase  the  computation  efficiency  by  dealing  with  a  small  computational 

domain.  The  present  phase  of  our  study  has  been  focused  on  the  application  of 

our  numerical  schemes  to  the  nonlinear  free  surface  problems  with  a  fixed, 
small  computational  domain.  There  are  three  major  tasks  involved  in  this 
study;  (l  )  investigate  the  implementation  of  the  open  boundary  condition  at 

the  outflow  boundary  to  allow  the  waves  to  pass  through  it,  (2)  improve  our 

numerical  schemes  so  that  accurate  numerical  solutions  can  be  obtained  using 
the  open  boundary  condition  at  the  outflow  boundary,  and  (3)  epply  the  methods 
developed  to  selected  nonlinear  free  surface  problems. 

We  have  studied  in  detail  the  method  to  implement  the  open  boundary 

condition  used  by  Chan  [ris].  It  is  based  on  Orlanski's  scheme  '.R19J  in  which 
the  Summerfeld  radiation  condition  is  applied  at  the  outflow  boundary  and  the 
phase  velocity  is  calculated  numerically.  We  have  conducted  numerical 
experiments  to  study  systematically  the  errors  of  Chan''s  numerical  scheme 
implement  the  advection  equation  and  to  study  ways  to  minimise  these  errors. 

We  have  studied  two  approaches  for  solving  potential  flows  for  nonlinear 

free  surface  wave  problems  using  a  fixed,  small  computational  domain.  In  each 

approach,  we  have  used  Orlanski's  method  to  implement  the  open  bo’cndary 

condition  at  the  outflow  boundary.  The  two  approaches  differ  in  the  treatment 
of  high  frequency  errors  that  appear  in  the  solution  and  that  inhibit  accurate 
treatment  of  the  open  boundary  condition.  In  the  first  approach,  we  used 
filtering  to  minimize  high  frequency  errors,  however,  it  was  found  to  be 
effective  only  at  an  early  time;  therefore,  it  is  suitable  only  for  obtaining 
transient  solutions.  In  the  second  approach,  we  modified  our  time  marching 
schemes  so  that  the  schemes  are  characterized  by  damping.  This  lamping  is 
more  effective  in  controlling  high  frequency  errors.  Transient  and  steady 
state  solutions  have  been  obtained  for  the  pressure  distribution  and  the 
accelerating  strut  problems. 

Pig.  10.1  shows  the  pressure  distribution  problem  with  the  outflow 
boundary  set  close  to  the  disturbance.  The  appearance  of  a  two-grid  inter/al 
wave  in  the  solution  of  this  problem  near  the  outflow  bo'ondary  is  shown  in 
Fig.  10.2.  The  eliminations  of  high  frequency  waves  in  the  solution  are 
shown  respectively  in  Figs.  10.3  snd  10.4  for  two  methods,  using  filtering 


Fig.  10.3  Elimination  of  2Ax  waves  in  the  solution  of  the  pressure 

distribution  problem  (Fr  =  1/ /l^)  at  t  =  1  by  using  filtering 


Fig.  10.4  Solution  of  the  pressure  distribution  problem  by  using  damping 


in  the  time-integration  scheme. 
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and  damping.  Fig.  10.5  shows  the  steady  state  solution  obtained  by  using 
lamping  for  a  fixed  computational  domain.  Fig.  10.6  shows  a  transient 
solution  of  wave  height  for  the  accelerating  strut  problem. 

The  results  of  the  second  phase  of  our  work  were  presented  at  the  Third 
International  Conference  on  Numerical  Ship  Mydrodynamics  TlO.i]. 
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1  ’  .  1 

Rerfomance  of  Instruction  Picalines 

An  instruction  unit 

pipeline  differs  frea  a  functional  unit  cipeline  in 

that 

it  is  subject  to 

delays  frea  system  phenotiena 

rnal  to  itself 

e. 

g.  Tianicry  ani  functional  units'*  .  "Ohus  a  sagmsnt  in 

an  i 

nstructicn  unit 

*  This  Ycr'K  -vaa  aucucrtal  'qy  the  -Toiat  ^er'ricas  Elactrcnics  Prc^rax.  ‘J.3. 
Arajy,  7.3.  7a'/y.  7.3.  Air  Porca'i  'Xiiar  contract  70001  i-’A-I -IIM,  by  tba 
7.3.  Oavy  uniar  :cntract  7'0003 3-30 -0-0335 ,  ani  by  t'a?  Havl att-'=ac;<a ri  Ooapany. 
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pipsline  nay  take  more  than  one  tine  unit  in  the  event  of  cache  nias,  nencry 
collision,  lata  iependency,  etc. 

An  exact  solution  for  pipeline  perforaance  -.ns  obtained  by  consiiering 
the  execution  time  of  any  given  segment  to  be  a  discrete  ranlom  variable  vhich 
is  governed  by  a  multinomial  iistribution.  The  expected  wait  time  of  an 
instruction  in  a  given  segment  can  be  recursively  defined  in  terms  of  the 
expected  wait  time  of  the  instruction  immediately  preceding  it.  If  the 
difference  of  these  wait  times  is  negative,  the  actual  wait  time  will  be  zero. 
Unfortunately,  this  solution  is  not  in  closed  form. 

The  effective  bandwidth  of  the  pipe  can  be  given  as  a  function  of  the 
expected  time  gap  between  the  ith  and  (i+1  dth  instructions.  The  stationary 
value  of  this  expectation  requires  that  i  be  infinite,  but  it  was  found  that 
for  values  of  i  as  small  as  3.  estimators  can  be  obtained  for  the  parameter 
ranges  of  interest  -with  a  high  degree  of  confidence. 

A  closed  form  solution  was  obtained  heuristically  for  a  pipe  consisting 
of  N  2  segments.  The  two  end  segments  are  variable- time  segments,  and  the 
middle  N  segments  are  all  constant- time  segments.  The  purpose  of  this  study 
was  to  determine  the  degree  to  which  a  finite-length  buffer  could  decouple  the 
bad  behavior  of  two  variable-time  segments.  This  solution  was  found  by 
matching  boundary  conditions  (on  U)  against  an  exponential  function  of  U,  and 
comparing  against  a  wide  range  of  simulated  values  (^29  points  in  all').  This 
solution  had  an  error  range  <  1^,  with  an  RMS  error  of  <  2^  over  the  simulated 
data  used  in  the  comparison. 

11.2  Test  Generation  for  Ucmplex  Digital  Systems 

A  new  and  general  fault  model  for  microprocessors  has  been  developed,  and 
teat  techniques  have  been  derived  to  detect  all  the  faults  in  the  fault  model 
(ll.S].  A  methodology  was  found  for  generating  combinational  structures  from 
high-level  descriptions  (using  assignment  statements,  "if"  statements,  and 
single-nested  loops)  of  register-transfer  level  operators  ;l1.l1.  The 
generated  structures  are  cellular,  and  are  interconnected  in  a  tree  structure. 
A  general  algorithm  has  been  developed  to  test  cellular  tree  structures  with  a 
teat  length  which  grows  only  linearly  with  the  size  of  the  tree.  It  can  be 
proved  that  this  teat  length  is  optimal  to  within  a  constant  factor.  'Vays  of 
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•naicirig  the  structures  sel f-checi'cin?  hare  ^Iso  been  fjunl. 


'1.5  Concurrent  Crrcr  Cetecticn  in  Arithmetic  ani  'Lc^ic  'Jnits 

A  new  metho  i  of  concurrent  error  letection  in  Computer  Arithmetic  ani 
Lo^ic  ’units  has  been  le’/isei.  Chis  methoi  exnloits  the  inherent  ph’/sioal 
red’unlancy  of  components  in  the  time  Icmain.  "This  techniq^'ue  is  applicable  to 
all  operations  of  a  typical  ALU.  In  this  sense,  our  method  is  superior  to 

Arithmetic  Codes,  which  are  applicable  to  only  the  arithmetic  oneraticns  of  an 
ALU.  ?’ur the rmo re ,  cur  methoi  also  works  with  many  types  of  carry  look  ahead 
circuits.  In  terms  of  the  error  coverage  capability,  cur  methoi  compares 
favorably  with  a  duplicated  ALU  method  of  error  detection.  In  terms  of  the 
hard'ware  cost,  o'ur  method  is  superior  to  both  the  Arithmetic  Codes  ani  the 
d’uplicated  ALU. 

' '  .  1  High  ? erfo rmance  and  Reliability  in  .Systolic  Structures 

Systolic  arrays  have  been  proposed  for  convenient  and  efficient 
implementation  on  a  i/L3I  chip.  However,  systolic  system  algorithms  orcposed 
do  not  use  all  the  inherent  parallelism  in  the  hari-.vare.  He  have  found 

tachniq,ues  which  increase  by  a  factor  of  two  or  three  the  performance  of 
systolic  systems  by  rearranging  the  se-puence  of  the  data  flow  and  making 

slight  modifications  to  the  original  systolic  system.  In  addition,  a  high 
level  error  detection  technique  has  been  found  for  use  with  svstolic  systems 
performing  matrix  multiplication.  This  method  does  not  use  m'uch  extra 

hardware,  but  will  provide  protection  against  both  transient  and  permanent 
erro  rs . 

'1.5  Heme ry  Addressing  Architecture 

?or  the  reduction  of  processor-memory  address  bani-.width,  we  have 
developed  an  address  nredicticn  mechanism.  .An  address  predictionstack  holds 
the  addresses  of  instructions  and  data  which  have  a  high  orebability  of  being 
referenced  soon.  T'ae  reduction  of  address  bandTldth  is  achieved  "oy  encoding 
positions  in  the  stack  struct'ure  and  sending  an  encoded  ocsition  to  the  memory 
subsystem  as  an  address.  The  memory  system  decodes  this  stack  oosition  and 
aoplies  it  to  anid-intical  stack  structure  maintained  -vLthin  the  memory,  ’-'sing 


the  information  in  its  stack,  the  memory  system 


can 


a  ccmnlete  address 
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A  full  ailrasa  is  transmittad  cn  a  ataclc  aisa. 

An  LRU  raplacament  policy  and  several  line  ipdate  acheaes  ware  evaluated 
by  aiaulation.  The  aiaulation  results  provide  an  information  theoretic  Icwer 
bound  on  the  avarage  number  of  address  bits  needed  per  memory  referanca.  This 
bound  requires  frequancy  encoding  of  the  stack  positions  and  miss  indicator. 
Conventional  fixed  length  binary  encoding,  which  is  more  easily  implementabla , 
requires  only  ona  or  two  extra  bits  par  address.  For  the  frequency  encoding, 
as  the  stack  depth  increased  and  the  stack  width  remained  relatively  small, 
the  average  number  of  bits  needed  per  memory  reference  continued  to  decrease. 
This  ohenomanon  did  not  occur  in  the  case  of  the  conventional  binary  encoding 
because  after  a  certain  noint  a  price,  in  terms  of  extra  bits  needed  for 
encoding,  was  being  paid  for  the  greater  number  of  stack  locations.  Stacks  of 
modest  size  using  conventional  encoding  can  reduce  the  average  number  of 
address  bits  needed  per  memory  reference  to  below  8  bits,  even  when  the  total 
address  space  is  very  large. 

11,6  Shared  C ache  Organizations  for  Multiple-Stream  Ocmputer  Systems 

Organisations  of  parallel-pipelined  multiple  instruction  stream 
processors  with  shared  two-level  memory  hierarchies  have  been  studied.  The 
cache  memory  interference  and  the  shared  cache  hit  ratios  in  such  systems  have 
been  investigated.  Performance  analysis  is  carried  out  by  using  discrete 
Markov  Ohain  and  probability  based  theorems.  Performance  is  evaluated  as  a 
function  of  hit  ratio,  the  number  of  processors,  the  number  of  pipeline 
segments  and  the  cache  organization  characterized  by  the  number  of  lines,  the 
number  of  cache  modules  per  line  and  the  cache  cycle  time.  Ocme  design 
tradeoffs  are  discussed  and  examples  are  given  to  illustrate  a  variety  of 
design  options  ^11  .12"^. 

11."^  Multiprocessors  with  P rivate  Cache  Memo ries 

We  have  developed  an  approximate  analytical  model  for  the  performance  of 
multiprocessors  -with  private  cache  memories.  All  private  cache  memories  are 
connected  to  a  single  shared  main  memory  via  an  interconnection  network.  We 
have  extended  an  earlier  analytical  model  to  include  complex  cache 
organizations  like  write- through,  load-through  and  buffered  write-back. 
Considering  the  complexity  of  such  organisations,  cur  approximate  analysis  is 
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I'iite  ^cc'i  ever  a  -^lie  ran^e  of  aysteTii  par'ixeters .  T-ieae 
cache  -cias  ratio,  blcc>  transfer  tide,  ani  sice 
incerccnnec ticn  networks 

■  '  . 3  Memo ry  cystsd  Architecture  for  ^ arallel  ^ recessing 

Research  is  in  progress  which  exrlores  the  iesign  of  a  xeTcry 
architecture  for  the  support  of  parallel  processing.  Ve  have  lefineh  a 
preliminary  memcry  system  architecture  which  consists  of  a  set  of  memory 
access  requests,  responses,  ani  the  semantic  action  associate!  with  each 
transaction.  This  architecture  will  ser/e  to  define  a  ccmm-onicaticn  protocol 
which  can  be  applied  to  a  variety  of  parallel  memory  structures  and  parallel 
interconnection  technologies.  The  parallel  architecture  allows  the  sharing  of 
data  in  a  very  intimate  fashion  as  in  the  traditional  shared  memory 
mul tiprocesso r ,  while  also  supporting  the  importing  'caching''  of  shared  data 
for  temporary  use  to  exploit  data  reference  locality.  A  paper  describing  this 
shared  memory  architecture  has  been  submitted. 

t 1 . 9  Computer  Architecture  for  High  bevel  Languages 

Research  has  been  in  progress  in  the  area  of  designing  directly 
interpretable  languages  'TIL's'i  for  high  level  language  support.  Techniques 
for  the  formal  derivation  of  DTL's  designed  to  efficiently  support  high  level 
languages  are  described  in  a  paper  which  has  been  submitted  for  publication. 
Results  indicate  that  it  may  in  the  future  be  possible  to  systematically 
design  efficient  DTL's  from  a  formal  specification  cf  the  high  level  language 
tc  be  supported.  We  are  currently  exploring  the  oossibility  cf  systematicallv 
designing  the  DIL  interpreter  from  a  formal  specification  cf  the  syntax  and 
semantits  of  the  DIL.  The  objective  is  to  explore  techniques  for  a  complete 
systematic  design  of  a  high  level  language  sunpert  processor.  Taus,  the 
design  process  would  not  only  yield  efficient  candidate  DIL's,  but  also 
hard-ware  structures  for  the  efficient  interpretation  of  a  selected  DIL. 


oarameters  include 
and  tyre  cf  the 


'i.'j  Logic  Talidaticn 


An  autemati:  design  verification  facility  is  being  developed  for  cur  ’/’LdI 
lAD  system.  We  have  made  progress  in  defining  the  orcblem  cf  ies'Lgn 


ve-ificati.cn. 


ore  olem 


rour  3Deci*ic 


sub-divides 
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(l  )  ietermining  appropriate  ieaign  language  restrictions  for  the  high-level 
f'onctional  iesign  language,  (2)  translating  the  high-level  proceiural 
iescription  of  the  systss  function  to  an  appropriate  lev-level  representation 
which  is  suitable  for  comparison,  (S)  translating  the  hierarchical  description 
of  the  system  structure  to  the  low-level  representation,  and  (1)  performing 
the  comparison  between  the  low-level  descriptions  of  the  system  function  and 
the  system  structure.  Work  on  the  fourth  sub-topic  may  be  aided  by  recent 
worK  in  the  literature  on  an  analytic  verification  algorithm  '/jhich  uses  binary 
decision  diagrams  to  represent  digital  functions.  We  hope  to  be  able  to 

incorporate  such  a  verification  method  as  part  of  our  hierarchical  Iesign 
verification  system,  since  this  approach  could  potentially  eliminate  the  need 
for  simulation  in  the  design  verification  process. 

11.11  Simulation  ackages  for  Complex  System  Design 

In  any  design  work,  simulation  becomes  an  important  part  of  the  design 
cycle.  Since  today's  logic  simulators  are  based  on  TTL  gate- level  structures, 
many  MOS  circuits  cannot  be  faithfully  modeled  using  them.  Sircuit  level 

simulators  require  a  high  computational  and  memory  space  penalty  and  cannot  be 
used  for  large  circuits.  A  logic  level  simulator  called  FETSIM  for  MD3 

integrated  circuits  has  been  developed.  The  simulator  takes  into  account 
slowdown  due  to  capacitive  loading  and  therefore  gives  accurate  tuning 

estimates  in  addition  to  logic  values. 

A  preprocessor  based  simulation  system  has  also  been  developed  that 
allows  the  specification  and  system  level  simulation  of  digital  systems.  The 
simulator  allows  hierarchical  specification  and  logical  connection  of  modules 
and  comprises  about  5000  lines  of  code  (including  direct  support  code).  It 
has  been  used  to  examine  multiprocessor  bus  protocols,  memory  eganizatien  in 
multiprocessor  systems,  and  failure  modes  for  a  multiprocessor  system.  Since 
the  system  can  deal  with  a  multilevel  specification,  it  is  fairly  efficient 
and  a  real  time  to  simulation  time  ratio  of  50,000  for  systems  with  a  minor 
clock  cycle  of  10  ns  has  been  obsen/ed.  Theoretical  models  for  such  systems 
are  'under  development,  and  will  be  calibrated  using  the  simulator. 
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‘2.1  Introduction 


Tie  prixary  objective  of  the  research  reported  in  this  section  is  the 
ievelopaient  of  efficient  computational  techniq^ues  and  the  analysis  of  the 
capabilities  of  various  models  of  computation.  We  are  concerned  with  the 
resources  -  such  as  time,  equipment,  memory,  interconnection  -  either  used  or 
needed  in  the  al-?orit’nmic  s-clution  of  given  problems.  This  dynamic  discipline 
-  concrete  computational  complexity  -  is  not  only  contributing  to  our  basic 
■under  3  tan  ling  of  computing,  but  it  is  extremely  relevant  to  actual  practice, 
both  in  hardware  and  in  software.  The  relevance  is  greater,  the  closer  the 
adopted  computation  models  are  tc  current  or  projected  computing  systems. 
Tnerefore,  it  is  only  natural  that  the  great  technological  innovation 
represented  by  Very-large-Jcale-Integrated  Jircuitry  has  already  had  a 

substantial  impact  on  the  discipline,  opening  new  horisons  and  pcsing 
challenging  problems.  Indeed  a  substantial  part  of  our  research,  organised 


below  in  four  subsections  VT'a-iH-’l  Computation  in 
Jeometr.",  ''Approximation  Algorithms",  and  "Jtorage 
Information  Retrieval"')  draws  its  motivation  from  this 


Tj^I",  "  Jcmputational 
and  Access  Costs  in 
Lmport.ant  technol'^gical 


*  This  fcr'<  .ns  suoicorted  by  the  Joint  Cervices  31ectroni;3  Program 
'■J.o.  Army,  7.3.  Wavy," ''1. 3.  Air  Force)  -under  contract  'IC01  and  by 

the  'lational  3:ience  Foundation  under  grants  '.103  "’3-1  5642,  MC3  31  -C5552,  .'dCS 
30-33354,  I3T  30-1  2240,  and  I3T  30-1  2242. 
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revolution. 

12.2  Parallel  Oomputation  in  VLoI 

Parallel  eonputation  provides  a  fundamental  answer  to  the  need  of 

increasing  computer  processing  po'^er.  Technological  advances  have  maie  it 
possible  to  realise  parallel  computing  systems.  Our  studies  in  this  area  are 
concerned  -rfith  the  development  of  realistic  computation  models,  the  conception 
of  powerful  interconnections  of  modules  (architectures)  realisable  in  the 
chosen  models,  and  the  analysis  and  design  of  algorithms  to  be  executed  on 
those  architectures.  Vithin  the  commonly  accepted  "synchronous  model"  of  '7L3I 
computation,  we  have  studied  the  realisation  of  algorithms  and  have  developed 
area-time  optimal  schemes  for  matrix  operations  [l  2. 1 d, 12. 1 9],  for  the  PPT, 
and  for  the  multiplication  of  integers  (12.201.  We  have  also  developed  a 

criti:iu3  and  an  appraisal  of  current  '/u3I  models  -  based  on  a  physical 

analysis  of  propagation  delays  -  and  have  begun  the  exploration  of 

three-dimensional  TuST  techniguea  (papers  in  preparation!. 

1 2. d  Oomputational  leometry 

Computational  geometry  is  relevant  to  a  number  of  applicational  areas, 
such  as  operations  research,  statistical  analysis,  and  design  automation.  Ve 
have  continued  our  lines  of  investigation,  and  have  developed  time-optimal 
algorithms  for  a  variety  of  problems  concerning  two  families  of  orthogonal 
segments  (12.6,12.?];  this  problem  is  very  relevant  to  VIj3I  design.  We  have 
also  begun  the  exploration  of  a  class  of  problems  of  planar  realizability  in 
the  presence  of  obstacles,  and  have  developed  a  time-space  optimal  algorithm 
for  ths  case  where  the  obstacles  are  modelled  as  parallel  segments  (l2.16]. 


Droximation  Xlgorithms 


Motivated  by  the  intractability  of  many  important  problems,  we  are 
concerned  with  the  development  of  efficient  approximation  algorithms.  Our 
work  with  two-dimensional  bin  pacicing  (l  2. 1 , 12. 4-,  12. 21  1  models  certain 
scheduling  applications  and  also  relates  to  VIilI  placement  problems.  Wj  have 
recently  studied  '/LoT  street  routing  problems  and  have  obtained  both  a  lower 
bound  on  the  number  of  tracks  required  and  have  developed  an  algorithm  which 
is  at  most  iH’l  worse  than  optimal  (papers  in  preparation!. 
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I  2. 1  5  toraige  ^^ni  A  :cv933  3o3t3  for  Information  Retrieval 

In.  any  infomation-ratriaval  ayatam,  affiniant  storage  ani  accessing  of 
iata  are  essential.  7e  are  stuiying  particular  lata  scructurea  ani  the 
3torage-aoce33  traieoffs  inherent  in  performing  associate!  operations.  Tnia 
i.3  leaiing  C(5  the  ievelopment  of  analysis  techniques  vhich  can  he  appliel  to 
various  retrieval  problems,  using  cost  measures  appropriate  to  the  particular 
applications . 

7e  have  also  stuiiei  '",12.2,12.151  a  general  class  of  resource  traieoffs 
that  arise  in  such  problems  as  parallel  sorting  algorithms,  linear  recursion 
schemata,  T^3I  layouts,  ani  searching  algorithms.  These  traieoffs  span  a 
variety  of  resources,  allowing  us  to  traie  time  against  space,  area,  program 
sice,  ani  n’umber  of  processors. 

A  iata  structure  with  local  access  permits  an  access  pointer  into  the 
structure  to  change  its  value  from  one  location  only  to  an  ai.lacent  location 
in  the  structure.  Multiheai  luring  machines  moiel  storage  ani  retriaval 
operations  in  these  iata  structures:  tree  machines  for  trees, 

multiiimension.al  machines  for  arrays.  Vs  have  usei  these  machines  to  ievelop 
new  methois  for  miniiicing  access  pointers  into  trees  ani  arrays.  Ivery 
multiheai  tree  machine  of  time  complexity  t(  nl  can  be  simulate!  on-line  by  a 
tree  machine  with  only  two  access  heals  in  time  2(t'[n''  log  t(n)  'log  log  t '  n'l . 
Uvery  multiheai  e- iimensional  machine  can  be  simulate!  on-line  by  a 
i-iimensional  machine  with  only  two  access  h^ais  in  time 
t'' n)*'*(l  *■  I  /  i  -  I  /ie'!  log  t'n')')  J  2. 1 1  2.  22  .  The  simulation  for  trees  ia 
optimal . 
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1  S .  1  'Natural  Language  'J  nde  re  tending* 

Research  has  continued  in  both  theoretical  and  practical  areas  relate!  to 
co-uputer  understanding  of  natural  language  and  the  obtaining  of  infomation 
fron  lata  bases.  Attention  was  focused  on  exploring  and  fomalizing  the 
relationship  between  language  and  perception.  ''“!vent  simul ation" ,  the  process 
of  constructing  a  sequence  of  sets  of  assertions  to  represent  the  meaning  of 
sentences,  was  proposed  and  investigate!  as  a  method  for  judging  the  relative 
plausibility  of  various  real  worli  interpretations  of  sentences  lascribing 

*  This  wor){  was  suoported  by  the  Office  of  Naval  Research  under  jontract  No. 
NOOOU-‘'5-G-0‘S1  2. 
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3cen03  -ani  actions;  avan^  simulation  vas  also  shown  to  be  useful  for  filling 
in  ^aps  in  tsxt  in  the  abssnce  of  "scripts".  ?rimiti-/’9  event  simulation  terms 
are  being  iesignel  to  provide  a  sound  basis  for  models  of  language 
acquisition,  concept  learning,  language  generation,  ani  "common  sense" 
J3.29].  The  primitives  are  based  in  part  on  our  study  of  Southeast  Asian  ani 
;Iative  .American  languages,  and  the  methods  they  use  to  eicpress  spatial  ani 
event-related  information. 

The  problem  of  eval-uating  natural  language  systems  was  investigated,  ani 
a  number  of  novel  evaluation  techniques  ani  dimensions  for  eval-uation  were 
developed.  Evaluation  methods  were  designed  to  help  in  the  grinding 
improvement  of  the  tested  systems;  methods  were  tested  on  PLANTES  ani  other 
natural  lang’uage  systems,  and  an  extensive  analysis  of  the  test  results  was 
carried  out  and  documented  [13.57]. 

JTL,  the  relational  data  base  query  language  system  used  in  ?LA?I33  was 
completely  rewritten  to  allow  "cooperative  ’•esponses",  probability,  and  the 
planning  of  optimal  data  base  search  strategies  [13.391. 

Vor'x  was  completed  on  3R0WSSH,  a  system  which  automatically  attempts  to 
find  interesting  patterns  and  facts  in  a  data  base  '_'3.3i1.  BROWSER  is  driven 
by  a  model  of  what  is  interesting,  in  terms  of  statistical  patterns,  measures 
of  significance,  and  models  of  the  real  world  actions  and  ob.jects  which  are 
reflected  in  the  data  base. 

We  have  also  been  working  on  a  computer  system  that  can  learn  from  its 
experiences  with  the  world,  ’yhile  designed  primarily  for  a  natural  language 
processing  system,  the  research  is  applicable  to  other  AI  domains  as  well. 
The  technique  is  called  Explanatory  Schama  Acquisition  and  is  a  kind  of  one 
trial  or  insight  learning;  with  this  technique,  a  natural  lang'uage  processing 
system  ‘will  be  able  to  learn  both  the  vocabularies  and  concepts  of  new 
knowledge  domains  directly  from  its  input  text.  Eor  example,  a  natural 
lang'uage  system  that  'xno-.fs  nothing  about,  say,  blac 'email  will  be  able  to  learn 
that  concept  from  a  input  story  descri'oing  a  blac'xmail  event  in  some  detail. 

Theoretical  work  has  been  carried  out  on  the  use  of  novel  "message 
passing"  algorithms  on  a  hexagonal  array  for  computing  high  level  visual 
features  'e.g.  s’/mmetry  ani  shapel  of  objects,  and  for  defining  appropriate 
neighborhoods  for  obiects  in  scenes  '_!3.S3^.  Results  here  'nave  broad 
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aignifioance ,  a.g.  for  VLSI  design  of  high  level  vision  processing  hardware. 

13.2^  Knowledge- based  System  for  Intelligent  Monitoring  and  Oiagnosis* 

The  goal  of  this  research  is  to  develop  a  knowledge-based  expert  system 
capable  of  performing  intelligent  diagnosis  on  engineer-designed  systems.  As 
traditional  expert  systems  have  been  criticised  for  lack  of  understanding 
about  their  working  domains,  our  effort  is  to  join  the  power  of  expert 
heuristics  with  the  completeness  of  general  domain  knowledge. 

In  our  approach,  an  expert  system  acts  as  the  front-end  for  a  diagnostic 
process.  It  efficiently  postulates  plausible  failures  to  explain  an  ob3er>/ed 
abnormality  based  on  the  powerful,  yet  not  perfect,  diagnostic  heuristics 
acq.uired  from  human  experts.  In  the  background,  postulated  fault  hypotheses 
are  scrutinized  by  a  rationalization  process.  Using  a  constraint  model,  the 
global  effects  of  each  hypothesis  are  thoroughly  deduced  and  then  compared  to 
obser/ed  sensory  data  for  consistency. 

By  emphasizing  coordination  between  proposal  and  verification  processes, 
a  basic  expert  system  can  be  enhanced  by  the  deep- thinking  intelligence.  The 
philosophy  of  our  intelligent  expert  approach  is  useful  not  only  for 
diagnosis,  but  for  other  expert  problem- solving  systems  as  well. 

M .  3  Automated  Air  Traffic  Control** 

Air  traffic  control  can  be  subdivided  into  a  set  of  tasks  each  having 
only  one  primary  concern.  3uch  tasks  include  the  avoidance  of  collisions, 
dangerous  weather  and  restricted  airspaces,  and  the  control  of  flow  rate  to 
airports.  Kowever,  a  set  of  techniques  to  perform  each  of  these  tasks 
independently  does  not  produce  a  system  which  performs  like  a  human 
controller,  '/fhat  is  required  is  an  architecture  that  supports  a  constructive 
dialogue  among  the  experts  comprising  the  system.  3uch  a  system  has  behavior 
that  depends  on  the  power  of  the  expert  interaction  language. 

*  This  work  was  supported  by  the  National  Aeronautics  and  Space  Administration 
under  contract  NASA  NAS  1 -?J. 

*•  This  work  was  suonorted  by  the  U.S.  Department  of  Transportation  inder 
Contract  XT  FA72-WA-noO. 
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Dar  resaarch  has  prolucal  tachniq^ues  li’o  those  usei  by  the  experts,  as 
■.fell  as  a  planning  process  involving  the  ase  of  criticism  ani  suggestions  from 
a  boiy  of  experts.  ilxperts  constantly  monitor  traffic  for  problems  in 
parallel.  3houli  an  expert  letect  a  problem,  it  plans  a  solution  and  submits 
it  to  the  other  experts  for  analysis.  This  process  may  cycle  before  a  final 
solution  is  adopted.  The  experts  share  a  common  conceptual  language 
describing  air  traffic  and  each  'xnoxs  its  own  relative  importance  in  the 
control  task.  Thus  the  experts  can  ''argue"  and  decide  priorities.  This 
architecture  emphasises  modularity  while  maintaining  the  cohesive  nature  of 
human  controller  performance. 

Still  the  system  'will  not  improve  -.with  time  -anless  it  has  a  learning 
capability.  Sur  extensive  rsvie'w  of  literature  on  learning  yielded  no  general 
theory,  although  some  simple  learning  has  been  demonstrated,  Vor;<  is 
continuing  to  develop  learning  techniques  to  improve  expert  performance 
through  experience. 

1 1 .  i  Theoretical  Computer  7 ision* 

k  book  on  the  general  subject  of  modeling  spatial  patterns  was  completed 
^'3.1  1.  'lost  of  the  book  is  devoted  to  a  discussion  of  the  properties, 
33mthe3i3,  exten3ion3  and  perspective  display  of  mosaics,  coverage  patterns 
ani  long  created  wave  patterns,  which  are  explored  as  models  to  planar 
te'ct'ures.  Model  based  texture  meas'ures  are  also  discassed. 

Ipproaches  to  planar  decomposition  for  hierarchical  image  representation 
were  studied  Square  ani  triangular  quad  trees  were  found  to  be  the 

only  feasible  methods,  having  the  same  comp-utational  complexity.  For  grid 
images,  the  type  of  grid  (triangular  or  square)  determines  the  appropriate 
:iu3i  tree  type  (triangular  or  square). 

1  ne'w  classification  of  low  le'/el  image  models  waa  proposed  ander  the 
categories  of  pixel  based  and  region  based  models  Traditional  models 

were  s'ur/eyed  -within  this  framework. 


This  work  was  supported  by  NS7  under  grant  TCl-dl  TolOS. 
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13.5  Pilot  Interaction  with  Automated  Airborne  Decision  Making  Systems* 

The  purpose  of  this  program  of  research  is  to  explore  the  issues 
surrounding  problems  associated  with  human-computer  interaction  in  flight 
management.  During  the  last  year,  three  projects  have  been  pursued: 
multi-task  decision  making,  measures  and  models  of  planning  behavior,  and 
onboard  computer-based  information  systems.  The  multi-task  decision  making 
project  has  been  involved  vrith  applying  the  modeling  and  function  allocation 
notions  developed  earlier  to  realistic  full-mission  simulations  [l3.1sl.  On 
the  basis  of  a  previous  review  paper  ’’^21,  the  planning  project  has  evolved  to 
a  aeries  of  experiments  to  study  the  planning  process  of  aircraft  pilots  in 
emergency  and  abnormal  situations  ’^13.25,13.261.  The  computer-based 
information  system  project  started  with  an  initial  fairly  simple  experiment 
ll3.1'^]  and  resulted  in  the  design  of  a  system  for  full-mission  evaluations 
J  3.27,13.23]. 

’  ^ Problem  Solving  Performance  in  Fault  Diagnosis  Tasks** 

The  purpose  of  this  project  is  to  develop  an  'understanding  of  human 
performance  in  fault  diagnosis  tasks  as  well  as  aids  for  use  in  training 
humans  in  such  tasks.  In  the  last  year,  two  large  experimental  studies  of 
mantenance  trainees  were  performed.  The  first  experiment  was  the  third  in  a 
series  of  studies  of  human  performance  using  computer  simulations  of  aircraft 
power  plants.  The  goal  of  these  three  experiments  was  to  study  how  suboptimal 
performance  is  influenced  by  a  variety  of  factors  131,13.331.  The  second 
experiment  performed  this  year  involved  a  second  study  of  transfer  of  training 
from  computer  simulations  to  real  equipment. 

13.7  Modeling  of  Human  Behavior  in  Beeking  and  Generating  Information*** 

This  project  is  studying  online  use  of  computers  by  scientific  personnel 
in  two  areas.  The  first  involves  the  use  of  online  editors  for  preparation  of 

*  This  work  was  sunnorted  by  the  National  Aeronautics  and  Bnace  Administration 
under  Grant  No.  NS5-2119. 

**  Tnis  work  was  supported  by  the  U.  3.  Army  under  contract  MDA  PD3-''9-T-31-2!  . 

•**  This  work  was  sunnorted  by  the  National  Science  Foundation  under  grant 
I3T-791  oi-^l  . 
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programs  aai  iocuments  ^13. 34^.  The  second  involves  the  online  searching  of 
fnll-text  bibliographic  lata  bases  '^13.361.  The  goals  of  both  of  these 
subprojects  include  developing  mathematical  models  of  human  behavior  that  are 
suitable  for  use  as  design  ails. 
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11.1  Introiuction* 

During  tha  1930-1931  time  perioi  the  Information  Retrieval  Research 
Laboratory  (IRRL)  conducted  a  namber  of  research  and  development  projects  and 
directed  the  operation  of  the  'Jniversity  of  Illinois'  online  search  ser’/ice. 
lajor  activities  included  analysis  of  database  data;  a  test  of  the 
feasibility  of  automatically  determining  the  overlap  between  bibliographic 
databases;  development  of  a  computer-readable  database  directory;  and 
management  and  direction  of  the  University  of  Illinois'  online  search  service. 

11.9  Analysis  o f  Database  Data* 

The  IRRL  maintains  a  body  of  information  about  commercially  available 
databases  called  the  Database  of  Databases.  He  have  continued  to  study 
various  statistical  characteristics  of  this  population  based  on  the  material 
in  our  database.  3uch  statistics  and  analyses  are  useful  both  to  researchers 
in  the  area  and  to  the  users  and  producers  of  databases. 

The  particular  information  developed  in  the  studies  performed  in  the  last 
year  considered  the  age,  size,  type  (scientific,  medical,  etc.'',  and  source 
(government,  private,  etc.')  of  the  databases  covered.  \t  present,  393 

*  This  wor!c  was  supported  by  the  'Jniversity  of  Illinois. 
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iatabases  ara  coverel  by  our  iababaaa  of  iatabaaaa,  but  any  given  atuiy  night 
i.nolule  only  a  specifiei  subset  of  those  lepending  on  the  needs  of  the  study. 
The  nost  recent  analyses  are  for  the  year  1  980. 

The  software  used  in  this  research  consisted  of  a  program  framework 
containing  slots  into  which  logic  could  be  inserted  to  perform  specific  data 
selection,  cleanup,  and  analysis.  The  data  was  first  extracted  from  the 
database  by  our  data  access  program  'JSSR  and  written  into  a  file  of  (tag, 
value)  pairs.  The  analysis  programs  read  that  file  and  selected  the  tags  of 
interest.  An  effort  was  mads  to  overlap  multiple  analyses  by  re-using  data 
selected  and  cleaned  for  prior  use. 

dross  tabulations  were  made  based  on  defined  groupings  or  actual  values, 
depending  on  the  data  involved.  Thus,  cross  tabulations  were  done  involving 
both  raw  year  of  origin  versus  grouped  sise  and  grouped  year  of  origin  versus 
grouped  size.  dtatistical  summaries  included  means,  ranges,  and  standard 
deviations  of  raw  data  and  correlation  coefficients  between  raw  data  and 
grouped  data.  Detection  of  bad  data  included  identification  of  records 
containing  illegal  or  missing  values. 

The  work  involved  in  this  project  also  developed  a  framework  for  further 
studies  of  the  databases  available  to  researchers  today.  'Ye  are  currently 
updating  the  database  and  improving  the  coverage  and  consistency  of  data 
fields  that  are  likely  to  be  important  in  future  studies. 

11.3  Integrated  'fan/Mac 'nine  Interface  to  Facilitate  Network  Resource 

'J  till  cat  ion* 

For  the  first  time  in  history  computerized  information  retrieval  is 
widespread  and  economically  viable.  Rvidence  of  this  is  sho-wn  in  the  volume 
of  machine- readable  records,  databases,  online  software  pac'<age3,  vendors, 
search  services,  and  searches.  The  major  portion  of  the  currently  published 
scientific  and  technical  literature  can  be  identified  through  computer 
searches,  because  the  references  are  in  computer- readable  form.  The  majority 
of  the  world's  currently  published  abstracting  and  indexing  literature  is  In 

*  T'nls  work  was  sunnorted  by  the  "lational  Dcience  Foundation  under  Irant 
No.  'HF  D3I 
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computar- readable  form.  There  are  more  thaa  SCO  publicly  available  databases 
containing  approximately  more  than  ISO  million  records  (citations).  More  than 
‘’5  million  of  these  records  are  now  available  through  the  various  online 
system  vendors  in  the  U.S.,  Canada,  and  Europe.  Some  of  the  major  venlors 
include:  the  Rational  Library  of  Medicine  (RLM)  ,  Lockheed  Information  Ser'/’ice 

(LIS),  System  Development  Corpicration  (SDG),  Bibliographic  Retrieval  Service 
[BRS),  the  Canada  Institute  for  Scientific  Information  (CISTI'',  and  the 
European  Space  Agency  (ESA).  These  online  service  organizations  use  a  variety 
of  sophisticated  online  software  packages  such  as  ELHILL,  DIALOG,  ORBIT. 
STAIRS,  CAR-OLE,  SRI?,  and  RECOR.  Seir/ices  from  the  databases,  through  the 
online  vendors,  are  provided  by  dozens  of  centers,  libraries  and  brokers  and 
they  conduct  nearly  6  million  retrospective  searches  per  year. 

■Jnfortunately,  there  are  no  standards  governing  the  format  of  databases, 
the  elements  and  vocabularies  included  in  them,  or  the  online  systems  for 
searching  them.  The  lack  of  standards  forces  searchers  to  become  familiar 
with  each  existing  database's  coverage  and  vocabulary  and  with  each  system's 
access  protocol,  system  features,  and  command  language  and  to  keep  up  to  date 
with  changes  made  in  all  of  them. 

Many  of  the  differences  and  variations  can  be  made  less  apparent  by 
developing  translators  or  transducers  for  converting  the  procedures, 
conventions,  and  terminology  of  one  system  into  equivalent  procedures, 
conventions,  and  terminology  of  other  systems.  These  converters  would  make 
systems  appear  alike  to  the  searcher  and  make  databases  function  as  if  they 
used  the  same  vocabulary.  Such  aids  would  make  the  differences  transparent  to 
the  users;  thus  they  can  be  called  Transparency  Aids.  An  integrated 
collection  of  transparency  aids  can  be  said  to  comprise  a  Transparent  System. 
The  aids  then  are  the  Transparency  Subsystems. 

Mith  the  reduction  of  differences  the  databases  and  systems  become  easier 
to  use;  the  searching  environment  becomes  more  user-oriented.  In  a  more 
user-oriented  environment  the  requirement  for  intermediary  searchers  decreases 
and  we  begin  to  approach  the  day  when  most  searching  will  be  done  by  end  users 
(those  that  pose  the  search  questions). 
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r'nia  r^3ear::h  prograx  atuiiei  xsthoia  to  inorsase  the  traaapareaoy  of 
retrieval  steps  by  xa-ciag  the  iiscrecsaess  of  separate  tas'o  less  obvious  aai 
by  xa'<iag  the  variety  of  systex  and  database  differences  less  obvious.  Within 
chis  project  we:  )  designed  a  generalized  xodel  for  a  Transparent  Systex, 

1 'i  exaxined  alternative  factors  affecting  the  Transparency  Interface,  ,1) 
built  ccxpcnencs  of  a  user-oriented  terminal,  and  (i)  integrated  example 
Transparency  Subsystems. 


The  Transnarent  Systex  is  comprised  of:  user  interface  programs; 

languages  used  by  the  Transparent  System;  transparency  ails 
-ous  f'unctions,  such  as  conversion  of  command  languages  or 
databases;  data  files  containing  descriptive  and  statistical 
actions  to  databases  and  networx  resources, 
t  Systex  lays  out  the  relationships  among  these 
iy  explores  advantages  and  disadvantages  of 
s.  Tevelopment  of  a  prototype  subsystem 
a  Transparent  System  by  showing  how  subsystems 
can  interact  in  performing  their  functions. 


command 

languages  used 

performing 

various  f'unctions 

selection 

of  databases; 

data  about 

databases;  the  C' 

The  model 

of  the  Transpa 

elements . 

The  alternatives 

various 

model  configurat 

iemonstrat 

S3  the  advantages 

The  generalized  xodel  includes  a  comprehensive  list 


existing , 


proposed,  and  potential  transparency  aids  along  -/rith  estimates  of  their 
intercommunication  requirements.  The  alternatives  review  then  considers  the 
band  width  of  communication  paths  required  for  various  configurations  and  the 
effects  on  response  time  and  performance  of  changes  in  configuration .  The 
Prototype  for  the  Searcher's  Workbench  'PTSW)  implements  a  structured  dialogue 
between  the  user  and  a  model  search  system.  The  dialogue  is  then  translated 
into  the  form  required  by  real-world  systems  to  be  used.  The  structured 
dialogue  integrates  access  to  two  existing  transparency  aids,  a  database 
selection  adviser  and  a  vocabulary  transformer. 


'1.1  directory  o f  Tom put er-Reai able  Bibliographic  Databases* 

The  Tom puter-Read able  Sib lie graphic  Databases  -  Directory  and  Data 
Sourceboo’x,  compiled  and  edited  by  Professor  ^dartha  3.  Williams  and  Sandra 
Rouse,  was  published  in  1 by  Knowledge  Industry  Publi cations ,  Inc.  "Ipdates 
to  the  Directory  were  issued  in  April  '  T’'’  and  Aoril  '  ^"’S.  A  new  ii-ect 

*  This  worx  is  partially  supoorted  by  Knowledge  Industry  PubI i ca tions ,  Inc 
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including  iatn  that  vas  relevant  aa  of  Decaubsr  1 973,  ^aa  published  in  autumn, 
!  97b.  Tlhe  computer  generation  of  the  printed  product  ia  aupportei  by  the 
American  Society  for  Information  Science.  ’Jpdating  of  the  underlying  database 
in  preparation  for  a  new  edition  to  appear  late  in  1931  is  now  underway. 

Data  Collection 

1 .  Database  processors  in  Europe  and  North  America  were  contacted  by 
mail  for  information  on  the  databases  they  provided  access  to. 

2.  The  13  page  questionnaire  designed  for  the  first  directory  was 
revised  to  achieve  a  more  consistent  format  and  expanded  to  include 
non-'J.S.  data.  .A  turnaround  questionnaire  is  being 
computer-generated ,  containing  the  information  in  cur  files  for  the 
databases  covered  by  the  first  directory.  Existing  data  entries  need 
only  to  be  verified  or  changed  and  missing  information  supplied  by 
the  database  producer. 

3.  The  questionnaires  were  mailed  to  the  database  producers  and  D.S.  and 
Canadian  producers  were  subsequently  contacted  by  telephone  for 
verification  of  the  data  and  assistance  with  questions  they  might 
have . 

i.  Information  on  the  availability  of  new  databases  was  obtained  from 
tne  published  literature  and  contacts  with  database 
processors/producers.  The  producers  of  new  databases  were 

inter’/iewsd  first  by  phone.  The  completed  questionnaire  was  then 
mailed  to  them  for  verification  and  if  necessary,  followed  up  by 
telephone.  Dnly  online  databases  that  were  publicly  available  were 
treated  fully.  Databases  available  for  internal  use  only  or  which 
are  only  accessible  offline  will  be  listed  separately. 

0.  Returned  questionnaires  were  checked  and  data  entries  standardized. 
Data  Entry 

1  .  Data  will  be  entered  online  into  the  database  of  databases  file  using 
the  D3IN  program.  Entries  can  then  be  printed  out,  verified,  and 
corrected  and  edited  as  necessary  using  the  DBEDIT  program. 


RST^rmL  i?33i:A3::K 
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Organization 

The  lirectory  i3  or^aniaei  alphabetically  by  iatabase  naae  ani 
an  intro iuction  anl  four  indexes:  Subject  Index, 

;ia!ne/Acron3rm/ Oynonym  Index,  Producer  Index  and  Processor  Index. 

Information  on  each  iatabase  follows  the  general  format; 

(l  )  Basic  Information 

'2)  ?roduc3r''Di3tributor''lenerator  Information 
(d)  Availability  and  Thar^es  for  Database  Tapes 
( V )  jub.ject  'latter  and  Scope  of  Data  on  Tape 
'3)  Subject  Analysis ''Indexing  Data 
(5)  Data  Tlements  Present  on  Tape 
(■'1  Database  Ser'/ices  Offered 
(s')  'Jser  Ails  Available 

Pro  iuc  tion 

The  bibliography  will  be  computer  typeset  using  IPPL  software  ani  a 
special  printer  on  IPPL's  microcomputer  system.  All  indexes  will  then  be 
generate i  automatically  from  the  iatabase  ani  typeset  similarly. 

'  T .  3  'Jniversity  of  I  lino  is  Onl  ine  Search  Bervice* 

During  the  past  year,  the  University  of  Illinois  Online  Search  Service, 
under  the  general  super-zision  ani  coordination  of  IPHL,  expanded  its  ser/ices 
to  the  students,  faculty  and  staff  of  the  university.  In  cooperation  with  the 
University  Library,  searching  is  now  available  at  11  different  locations 
aro'xnd  the  campus,  with  further  growth  pending.  Locations  where  searching  is 
offered  to  'university  students,  faculty  and  staff  include  IPRL,  Reference  Room 
of  the  ^dain  Library,  Ohemistry  Library,  Biology  Library,  Labor  ani  Ini'ustrial 
Relations  Library,  Education  and  Social  Bciences  Library,  Agriculture  Library, 
Uatural  Uistory  Survey  Library,  Physics  Library,  Veterinary  "dedicine  Library, 
Health  Sciences  Library,  Agricultural  Economics  Reading  Room,  Oeology  Library, 
and  the  Oeological  Survey  Library. 

*  This  wor!<  was  supported  by  the  State  of  Illinois. 
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Caline  systems  available  through  the  Search  Service  incluie  Lockheed's 
DIALCS,  System  Development  Corporation's  DRBIC,  the  Bibliographic  Retrieval 
Ser'/ice  system,  and  the  National  Library  of  Medicine  MEDLARS  system. 
Individual  databases  available  from  these  systems  total  176. 

Jse  of  all  systems  luring  ?Y  SI  is  given  in  the  follo'.rLng  table.  The 
costs  given  represent  the  total  amount  paid  by  the  iJniversity  to  the  online 
vendors  concerned.  Approximately  90"?  of  this  vas  recovered  from  the  users  of 
the  Search  Service  through  the  charging  of  fees. 

■JSS  DP  ALL  SYSTEMS  ?Y  SI 


Connect 

Hours 

Costs 

3R3 

S64.0? 

26,666.  hi 

LRS 

1  26 .  S6 

11  ,426.91 

3  DC 

1!  .73 

!  ,610.00 

NLM 

6.  60 

166.00 

TOTAL 

1007.16 

11  ,460.36 

The  total  connect  hours  used  represent  approximately  24DD  individual 
searches,  a  search  being  defined  as  a  single  session  at  the  terminal, 
regardless  of  the  number  of  databases  or  systems  accessed.  Df  these,  veil 
over  SDo  vere  conducted  and  paid  for,  at  least  in  part,  by  'Jniversity  faculty, 
students  or  staff.  The  remainder  vere  conluctel  by  IRRL  or  library  staff 
either  as  practice  as  they  learned  to  use  the  system  or  for  general  reference 
or  '.'■erifi cation  -rork. 

Beyond  providing  general  direction  and  coordination,  IRRL  performs  the 
follo'Ving  specific  tasks  related  to  the  Online  Search  Service:  negotiation  of 
contracts  vith  online  vendors;  training  of  new  searchers;  consultation  on 
system  problems  or  difficult  searches;  keeping  statistics  on  all  aspects  of 
the  Search  Service;  accounting  work  associated  vith  collecting  fees  and 
paying  bills. 
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1  V . 6  PatibasQ  Analys-  i* 

IRRL  is  conducting  ongoing  analyses  of  the  MARC  larabase.  This  vork 
involves  processing  the  annua'  files  distributed  by  the  Library  of  Congress 
and  generating  a  detailed  statistical  analysis  of  the  rate  of  occurrence  and 
sise  distribution  of  each  field  of  the  records.  Che  resulting  tables,  vhich 
have  appeared  in  the  Journal  of  Library  Automation,  shov  for  each  field  ''  ) 
the  percentage  of  records  in  vhich  it  appears,  "’)  the  rate  of  occurrence 
;vithin  records  Yhere  it  appears,  and  ( ^ the  length  of  the  field  in 
characters.  ?or  (l")  and  (d)  va  report  average,  minimum,  and  maximum  values 
and  standard  deviations. 

Che  distribution  of  records  in  Dewey  and  LC  classes  is  also  reported 
along  with  average  record  lengths  for  each  class.  Cummary  statistics  are  also 
listed  separately.  Finally,  ve  report  the  changes  between  reporting  periods 
(showing,  for  instance,  which  fields  are  used  more  now  than  at  the  last 
report)  . 

We  have  also  beg’un  a  second  project  using  the  same  computer  programs  to 
analyme  the  J0N3ER  database  for  the  Council  on  Library  Resources. 

These  data  allow  users  of  the  databases  to  optimise  their  applications 
for  the  actual  characteristics  of  the  databases. 

1 i . 7  Automatic  Detection  of  Database  Overlap** 

This  research  program  addresses  the  problem  of  overlap  coverage  of 
journal  articles  among  multiple  databases.  Overlap  leads  to  duplicative 
processing  of  material  by  database  producers;  duplicative  processing  and 
storage  of  material  by  online  and  batch  vendors  of  database  ser'/ices;  and 
retrieval  of  duplicative  material  by  users  of  databases.  The  net  result  is 
reflected  in  costs  associated  with  the  databases  and  an  unnecessary 
expenditure  of  time  on  the  part  of  users. 


*  This  work  was  supported  by  the  council  on  Library  Resources. 

**  This  work  was  supported  by  the  National  Ocience  Foundation  under  O-ant 
W37  loT  70-c’  :)l  3. 
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The  true  extant  of  tha  problem  ia  uaknovm.  Re'/’eral  atulies  have 
ailreaaei  tha  ^lueation  of  journal  covarage  overlaps,  some  have  aiireaaad  tha 
ijueation  of  article  overlaps  in  narrow  subject  areas,  ani  one  has  attacked  the 
problam  of  identifying  maximum  possible  overlaps  among  specific  databases, 
^one  has  studied  the  problem  of  determining  the  actual  overlaps  among  database 
services . 

■rfe  are  developing  a  methodology  for  automatically  determining  overlap 
among  databases  and  tasting  and  implementing  this  methodology  on  seven 
databases:  AOTITOLA,  3I03I3  Previews,  CACon,  Rxcerpta  'dedica,  Medline, 

3cience  Citation  Index  (3CI),  and  Social  Science  Citation  Index  (3371).  The 
methodology  involves  analysis  of  the  databases  to  determine  the  data  elements 
present  in  and  common  to  multiple  databases  in  order  to  determine  effective 
keys  for  matching  records  in  multiple  files.  The  matching  process  will  be 
done  in  at  least  two  passes.  The  first  pass  or  broad  screen  will  use  a  match 
key  with  high  discrimination  power  to  identify  candidate  duplicate  records,  i 
set  of  finer  screens  will  be  used  for  verifying  pairs  of  records  as 
duplicates.  The  second  pass  will  employ  Marrison  keys  '(-which  are  bit  string 
representations  of  titles)  and  Hamming  distance  techniques  (for  overcoming 
slight  variations  in  representation  of  titles)  together  with  matches  on  other 
elements  as  needed. 

Results  will  be  tested:  (l  )  by  using  manual  verification  of  a 
statistically  representative  sample  in  order  to  determine  mismatches,  -and  (3'' 
by  using  a  kno'wn  set  of  duplicate  records  to  see  if  they  are  identified  by  the 
screens.  The  eval'oation  in  terms  of  mismatches  (wrong  matches)  and  missed 
matches  will  provide  a  measure  of  the  effectiveness  of  the  techniques. 
Further  eval-uation  in  terms  of  gross  costs  will  be  carried  out  in  order  to 
determine  the  efficiency  of  the  techniques. 

The  study  will  produce  overlap  statistics  among  seven  databases,  data 
element  frequency  statistics,  and  a  methodology  for  duplicate  detection  that 
might  be  used  by  database  producers  or  database  vendors.  The  resultant 
methodology  could  be  used  as  a  oreprocessing  step  to  reduce  storage  and 
processing  costs  or  as  a  post  retrieval  step  to  eliminate  the  duplicative 
outout  delivered  to  users. 
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1  X. 

Duriri'?  thia  first  yesr  of  tha  project  tapes  hava  been  acquire!  frrni  the 
iatabaae  proiucera  covering  the  '173  -  i  I'^q  perioi.  lanple  aubaeta  have  been, 
real  ani  prelminary  planning  for  the  file  atruoturea  ini  processing  flow  has 
been  iona.  The  contents  of  the  samples  have  been  stuiie-i  statistically  to 
ilentify  iata  elements  useful  in  the  matching  process  ani  to  allow  estimation 
of  storage  requirements. 


15.  COMMUNICATIONS 


Faculty  and  Senior  Staff 


T.  Basar 

B.  S.  Hajek 

M.  B.  Pursley 

T.  U.  Basar 

J.  G.  Kieffer 

D.  V.  Sarwate 

A.  H.  Haddad 

R.  J.  McSliece 

S.  Tantaratana 

H.  V.  Poor 

Graduate  Students 

B.  Aa zhang 

F.  D.  Garber 

G,  M.  Stadler 

C.  J.  Alston 

E.  A.  Geraniotis 

W.  E.  Stark 

A.  C.  Blumer 

J.  3.  Lehnert 

T.  G.  Van  Loon 

C .  I .  Chang 

R.  G.  Ogier 

!<.  3.  Vastola 

R.  Cruz 

Y.  Rivani 

3.  Verdu 

?.  Enge 

M.  3.  Wallace 

15.1  Multiple-Terminal  Digital  Communications* 

This  is  a  major  research  area  vfhich  includes  various  problems  in 
spread-spectrum  communications  and  random-access  communications  networks.  The 
topics  that  were  investigated  during  the  past  year  include  the  performance  of 
spreai-spectr’ura  multiple-access  communications  systems,  the  effects  of  fading 
on  spread-spectrum  communications,  communications  in  the  presence  of  jamming, 
conflict  resolution  in  random-access  communications  systems,  packet 
communications  systems,  and  signal  detection  in  multiple-access  channels. 
Progress  in  each  of  these  areas  is  described  in  the  subsections  that  follow. 


*  This  research  was  supoorted  by  the  Army  Research  Office  under  grant 
DAAG29- ■^5-1-01  U  and  contracts  DAAt}29-31 -'<-0062  and  DAAD2 9-31 -<-0064.  by  the 
Joint  oer/ices  Electronics  Program  under  contract  N0001 1-79-0-0421,  by  the 
Naval  Research  Laboratory  under  contract  N0001 4-30-C-0302,  and  by  the  National 
Science  Foundation  under  grant  ENG  73-06650. 
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1  5  •  1  . 1  Spread- 3  pec  tram  '*!ul‘:iple- Access  lo'cnmnications  Systems 

’Saring  the  past  year  m  have  obtained  many  ne>r  results  on  the  problem  of 
evaluating  the  average  probability  of  error  for  asynchronous  direct -sequence 
spread-spectr'jm  multiple-access  (SS'l.A)  comm'unications  systems.  Because  of  the 
nature  of  the  multiple-access  Interference,  standard  numerical  integration 
techniques  require  an  excessive  amount  of  computation  in  order  to  provide 
approximations  that  meet  a  prespecified  degree  of  accuracy.  We  have  developed 
upper  and  lower  bounds  which  can  be  eval'uated  with  a  moderate  amo'int  of 
c^'^putation  for  direct-sequence  33M.\  systems  with  up  to  four  or  five 
simultaneous  transmitters.  The  bc'unds  are  presented  in  '.Bo]  along  with 
n'umerical  results  which  illustrate  their  application. 

Ve  have  also  levelopei  acc'urate  approximations  to  the  average  probability 
of  error  for  direct-sequence  S3‘!A  systems.  These  approximations  are  based  on 
the  integration  of  the  characteristi  c  f'anction  of  the  multiple-access 
interf erence .  7ne  primary  aivantage  of  these  approximations  is  that  the 
amoimt  of  computation  required  for  their  evaluation  grows  only  linearly  as  the 
number  of  simultaneous  transmitters  increases,  whereas  the  amount  of 
computation  required  by  standard  numerical  integration  'and  most  other  'ciown 
approximation  techniques)  increases  exponentially.  In  [I’S.Sh"'  we  present 
results  on  the  approximation  of  the  average  error  probability  for  binary  ?S'< 
direct-sequence  S3'-IA  systems.  Extensions  to  quaternary  systems  with  possibly 
offset  in-phase  and  quadrature  compouents  are  given  in  ’".Bl].  In  ’.B' ]  we  also 
allow  chip  waveforms  of  arbicrary  shapes  so  that  the  model  includes  as  special 
cases  the  forms  of  direct-sequence  modulation  Scno'wn  as  quadriphase-shift 
!<eying  (??3X),  offset  QP3X,  and  minimum-shift  Xeying  (BSX).  The  approximation 
method  that  we  present  gives  sufficient  accuracy  for  nearly  all  appl ications . 
However,  if  necessary,  even  greater  accuracy  can  be  achieved  by  using  thia 
approximation  to  obtain  an  expansion  point  for  a  Taylor  series  representation 
of  the  actual  probability  of  error.  By  employing  this  combination  of  the 
characteristic  function  method  followed  by  a  series  -  expansion  method,  we  can 
obtain  any  prespecified  degree  of  accuricy. 


Turther  results  have  also  been  obtained  on  the  performance 
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Jo ubl 5- frequency  ternis  ere  oonailered.  In  ']l5.2]  the  signal- to-noise  ratio  is 
analyze'!  for  general  quaternary  D3/33!tA  systems.  Numerical  results  are  given 
for  3P5'<,  offset  QPSK,  and  MSX  systems.  In  flo.io]  optimal  phases  (:vith 

respect  to  signal- to-noise  ratio)  are  Jeterminei  for  m-sequences  of  periods  51 
and  55.  "Jsing  the  phases  Jeterminei  in  [l5.10]  and  the  system  analysts  of 
M5.2I  we  have  shown  that  the  signal- to-noise  ratio  can  be  significantly 

increased  by  careful  selection  of  the  phases  of  a  set  of  signature  sequences. 
Our  subsequent  work  on  the  probability  of  error  (e.g.,  ^.15.59],  [pll.  sind 
[95])  has  shown  that  similar  improvements  in  the  bit  error  rate  are  also 

obtained.  The  effects  of  the  shape  of  the  chip  waveform  on  the  performance  of 

direct-sequence  S3'dA  systems  has  been  investigated.  Results  on  the  sine  pulse 
vs.  the  rectangular  pulse  are  given  in  '[15.2]  and  [R1],  and  a  more  extensive 
study  covering  various  pulse  shapes  is  presented  in  ']l5.54l. 

Nearly  all  of  the  analytical  results  on  the  performance  of 
lirec t-sequance  33'dA  communications  systems  are  restricted  to  systems  for 
which  the  number  of  chips  per  data  bit  is  an  integer-multiple  of  tha  period  p 
of  each  signature  sequence.  That  is,  there  are  an  integral  number  of  periods 
of  tha  signature  sequence  in  each  data  bit.  7or  certain  applications, 
however,  tha  necessity  to  discriminate  against  multipath  and  to  reiuca  the 
probability  of  intarception  dictates  the  use  of  very  long  signature  sequences. 
Typically  this  will  in  turn  necessitate  the  use  of  direct-sequance  modulation 
with  multiple  data  bits  per  period  of  the  signature  sequence,  in  which  case  'I 
may  ba  much  smaller  than  p.  In  a  DS/SSMA  system  -with  multiple  data  bits  per 
parioi,  the  performance  depends  on  the  partial  correlation  properties  of  the 
signature  sequences  ' rather  than  tha  aperiodic  correlation  properties  as  in 
the  special  case  p  =  'I.  In  '[1 5.52]  we  analyze  the  performance  of  a  T3'55NA 
system  for  which  p  >  'I  and  N  and  p  are  relatively  prime.  This  analysis  takes 
into  acco'unt  the  effects  of  partial  correlation.  In  particular  we  compare  the 
multiple-access  capability  of  two  classes  of  73/33^5  systems  with  N  chips  per 
bit  when  the  performance  measure  is  the  mean-square  value  of  the 
crosscorrelation.  Ne  find  that  the  multiple-access  capability  is  less  for  the 
class  in  which  p  and  'I  are  relatively  prime  than  for  the  class  with  one  bit 
per  period.  This  is  primarily  because  of  improvements  that  result  from  the 
optimisation  of  the  phases  of  the  signature  sequences  in  the  case  p  *  N. 
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In  aliition  to  tha  ■^o^ic  iescribei  above  on  iirect-aeqaance  3T14 
ooxxunications  aystajts,  aevaral  investigations  of  freqency-hoppe'l 
ooxmunioations  systeas  nave  also  been  coapletai.  Both  f ast-frequenoy-’nopping 
'_!o.?3]  and  slow- frequency-hopping  have  been  consiierel.  ?or  slow-frequency- 
hopping  we  have  investigate!  DP3K  [l5.3o]  and  7SK  Jo-io]  lata  aoiulation,  an! 
results  have  been  obtained  for  hybrid  direct-sequence  frequency-hopped 
spread-spectr'um  systens  '_15.3o]. 

15.1.2  Effects  of  Fading  on  3  p  read-3  pec  tm-in:  lonimunications 

Several  of  our  current  research  projects  are  concerned  -with  the  effects 
of  fading  on  the  perfomance  of  spreai-spectr’uni  coTim’unicat ions  systenis.  The 
fading  nay  be  either  selective  or  nonsalective  and  either  Rician  or  Rayleigh. 
Roth  frequency-hopping  and  direct-sequence  spread- spec tr'jia  systens  have  been 
considered . 

In  :J5.dO]  bounds  and  approicinations  are  derived  for  the  average 
probability  of  error  in  an  asynchronous  73X  slow-frequency-hopped  33'U 
conn’unications  system.  Three  types  of  channels  are  considered:  the  nonfading 
additive  white  diussian  noise  channel,  the  nonselective  Rician  fading  channel, 
and  the  doubly-selective  wide-sense-stationary  uncorrelated-scattering  ^V331f3' 
fading  channel.  Tlie  two  types  of  hopping  patterns  considered  are 
deterministic  frequency-hopping  patterns  based  on  Reed-3olomon  codes  and 
random  frequency-hopping  patterns.  Results  for  nonselective  Rayleigh  fading 
are  obcained  as  a  special  case  of  the  results  for  nonselective  Rician  fading. 

In  we  investigate  the  performance  of  slow-frequency-hopped  lo'dl 

using  Ii?3K  data  modulation.  The  channels  considered  are  Rayleigh  '.ifRoIJS  fading 
channels.  Frequency-selective  as  wall  as  time-selective  fading  channels  are 
considered  for  systems  with  various  data  pulse  waveforms  and  channel 
correlation  functions.  Analytical  results  for  the  average  probability  of 
error  are  presented,  and  numerical  examples  of  the  most  important  cases  are 
given. 

The  performance  of  direct-sequence  spread-spectrum  comm'xn Lcations  In  a 
specular-multipath  fading  environment  is  investigated  in  ''j5.5u  The  bounding 
methods  of  ,R5!.  which  were  developed  for  mul  tipi  e-access  in‘:?rf?’'?nce,  are 
applied  to  specular-multipath  interference  in  'J5.5'''.  It  is  sho-wn  that  our 
approximation  to  the  probability  of  error  that  is  based  on  the  signal- to  -  no ise 


84 


I").  lOMMUNIOATIOl^S 

ratio  [15.2]  also  ^ives  very  good  results  when  applied  to  the 

speoular-aiultipath  probleni. 

-he  improvements  that  oan  be  obtained  by  employing  error-correcting  codes 
on  a  slow-frequency-hopped  3S'tA  system  have  been  investigated.  One  conclusion 
of  our  results  (e.g.,  [l5.40])  is  that  in  a  Rayleigh  fading  channel,  bit  error 
rates  on  the  order  of  10”’  to  10“’’  cannot  be  attained  in  most 
slow- frequency-hopped  3S'dA  communications  systems  without  the  use  of 
error-correcting  codes.  However,  the  application  of  certain  error-correcting 
codes  -will  bring  the  bit  error  rate  down  to  this  range. 

13.1.3  Oommunications  in  the  resenca  of  J amming 

Last  year  we  began  an  investigation  of  the  problem  of  comm’unicating  in 
the  presence  of  hostile  .jamming.  Ve  have  formulated  the  problem  in 
game-theoretic  terms,  using  Shannon's  mutual  information  as  the  payoff 
function.  Ye  showed  that  under  certain  restrictions  the  players  in  this  game 
(the  coder  and  the  jammer)  have  simultaneously  optimal  ( saldlepoint) 
strategies.  Tnese  strategies  are  generally  memoryless,  and  so  in  an 
information-theoretic  sense  the  commonly-used  interleaving  strategies  are 
superfluous.  He  have  also  determined  that  in  the  presence  of  binary  coding, 
uniform  laussian  jamming  is  in  general  inferior  to  pulsed  jamming.  However, 
we  found  that  provided  a  modest  amount  of  bandwidth  expansion  is  tolerable, 
suitable  coding  can  neutralise  a  pulsed  jammer.  These  results  are  presented 
in  J3.1.91  3,y3tems  which  utilise  phase- shift-xeying  modulation. 

For  the  case  of  frequency  shift  xeying  (?3K)  modulation  we  investigated 
the  capacity  in  the  presence  of  partial-band  jamming  and  nonselective  fading. 
He  found  that  there  is  an  optimal  positive  code  rate  that  minimizes  the  energy 
per  bit  required  to  achieve  reliable  communication.  Also,  we  showed  that  for 
low  code  rates  ['  .30  suffices)  'uniform  jamming  is  optimal  for  the  jammer. 
Finally  we  analyzed  the  performance  of  a  specific  coding  scheme  for 
nartial-band  jamming.  These  results  are  presented  in  '13.33). 

13.1.1  Tonflict  Resolution  in  Random- Access  Tommunication  Hetworixs 

Tonflict  resolution  consists  of  strategies  which  allow  remote  stations  to 
share  a  comm'unication  channel  of  li.mited  capacity  in  a  decentralized  fashion, 
■fni ch  strategy  is  appropriate  depends  on  the  demands  of  the  stations  and  on 
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the  feeibac'.^  tafomation  available  to  the  stations.  Ve  have  ierived  upper 
bo'inis  on  ranioTi- access  throughput,  anl  ve  have  ievelopei  ranioni- access 
strategies  for  situations  vhen  feeiback  is  severely  liiitei  and  'under  several 
types  of  station  deniands. 

1  bounding  technique  inspired  by  the  state-space  approach  of  stochastic 
control  is  introduced  in  [l5.5j.  It  is  used  to  establish  an  upper  bound  of 
.5126  packets/ slot  for  the  throughput  of  a  time-slotted  multi-access  broadcast 
channel  subject  to  an  infinite  population  of  user  stations  (xodeled  by  a 
Poisson  process)  -using  feedback  information  consisting  of  a  3,  1  ,  or  e  to 
denote  a  slot  -.vith  sero ,  one,  or  at  least  tvfo  packets,  respectively. 

In  '.15.15]  and  :_l5.d2]  a  class  of  simple  recursive  retransmission  control 
strategies  is  given  for  the  random- access  broadcast  channel ,  again  subject  to 
an  infinite  population  of  users.  They  provide  stable  throughput  up  to  e"' 
packets/slot  under  a  variety  of  feeiback  models,  and  they  are  insensitive  to 
fluctuations  in  traffic  intensity.  Analytical  and  numerical  performance 
evaluations  are  given. 

A  class  of  random- access  strategies,  termed  acknowledgement  based 
retransxiss Ion  control  (A3HC)  strategies,  are  analysed  in  for  both  finite 

and  infinite  populations  of  users.  The  only  feedback  needed  to  implement  an 
A3!lj  strategy  is  that  ac'soiowledgement  of  a  successfully  transmitted  packet 
must  be  given  to  the  transmitter.  3n  the  basis  of  an  equilibrium  analysis  and 
the  local  Poisson  approximation  introlued  in  ]'5.i5l  it  is  s'no'vn  that  the  A3P2 
strategies  provide  satisfactory  stable  throughput  if  there  is  an  appropriately 
chosen  upper  bo'und  on  the  n-umber  of  retransmission  attempts  allowed  to  any  one 
user.  Expressions  for  delay,  throughput,  and  probability  of  rejection  are 
found  and  suggestions  for  the  choice  of  parameters  are  given. 

Another  type  of  -user  demand  is  addressed  in  [l  5.4-1].  '-{ere  the  users  each 
have  a  steady  stream  of  packets  to  transmit  and  thus  a  TT'dl  access  mode  is 
appropriate.  However,  if  there  is  limited  comm-unication  among  the  users,  and 
moreover,  if  the  population  of  active  'users  is  time-varying  (as,  for  example, 
in  a  mobile  communication  network)  then  some  strategy  for  lo eking  into  Tl'dA 
with  limited  information  must  be  -used.  Preliminary  versions  -ef  such 
strategies  and  performance  evaluations  are  given  in  _15.44l. 
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15.1.5  jeneral  '1  umerical  Methods  and  Bounding  Techniques 
for  ? acket  Tooiunication  Bystems 

One  numerical  technique  ve  have  investigated  is  the  method  of  phase- type 
Markov  processes.  Phase- type  processes  are  Markov  processes  t  >  0 

vhere  the  component  (^-t)  is  itself  a  Markov  process,  called  the  phase 
process.  For  example,  may  represent  the  number  of  packets  queued  for 

transmission  at  a  communication  station,  and  may  renresent  the  phase  of  a 
TOMA  protocol,  or  may  model  fluctuations  in  incoming  traffic  intensity,  or 
both,  numerical  results  of  this  applicat’on  are  reported  in  ^56]. 

Nev  basic  results  in  the  computational  and  conceptual  aspects  of 
phase- type  Mar'rov  processes  are  presented  in  [l5.1i].  In  particular,  a 
systematic  method  for  treating  general  boundary  states  is  found. 

The  method  of  phase-tyqje  Markov  processes  is  also  applied  in  '".M6]  to  find 
the  average  queueing  delay  in  a  link  of  a  communication  network  when  the 
incoming  traffic  is  obtained  from  a  Poisson  stream  by  a  generalized  round 
robin  routing  rule.  'Jnder  IMP  routing,  successive  packets  coming  into  a 

node  are  routed  out  of  the  node  in  a  pre- determined  sequence  according  to  the 
order  of  arrival.  The  delay  for  3RP  strategies  is  often  signif icantly  smaller 
than  for  the  standard  independent  Bernoulli  splitting  routing  strategy. 

A  general  bo'onding  technique  based  on  drift  analysis  is  given  in 
Simple  conditions  are  given  to  ensure  that  drift  implies  convergence.  One 
application  in  .P2j  is  a  proof  of  channel  stability  under  the  retransmission 
policies  introduced  in  [l5.15].  A  stability  result  and  bounds  on  waiting 
times  in  OI/l/l  queues  are  also  presented  in  '"pp]. 

In  Pd.]  it  is  shown  that  among  all  arrival  processes  ]not  necessarily 
stationary  or  renewal  typel  for  an  exponential  server  queue  with  specifiel 
arrival  and  ser'/ice  rates,  that  the  arrival  process  which  minimizes  the 
average  delay  and  related  quantities  is  the  process  with  constant  interarrival 
times.  The  proof  is  based  on  a  newly  discovered  convexitv  property  of 
exponential  server  queues  which  is  of  independent  interest.  Bince  the  traffic 
need  not  be  stationary  or  renewal  type,  the  theorem  provides  lower  bo'unds 
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(which  are  readily  computable  by  exiating  methods')  to  the  average  delay  in  a 
network  of  queues  under  rather  general  routing  disciplines.  Some  related 
bounds  and  comparison  with  some  numerical  results  are  presented  in 

15.1.6  Signal  Detection  in  Multiple- Access  Noise 

A  project  in  multi-user  communications  that  has  been  initiated  recently 
is  the  study  of  optimim  signal  detection  procedures  for  application  in  the 
multi-user  environment.  Such  procedures  are  inherently  nonlinear  because  of 
the  non-Caussian  nature  of  the  multiple-access  interference  in  the  channel. 
Thus  the  initial  goal  of  this  study  is  to  determine  the  degree  of  performance 
improvement  exhibited  by  optimum  receivers  over  the  more  conventional  linear 
(matched  filter)  receivers,  and  to  identify  those  situations  in  which 
performance  is  sufficiently  enhanced  to  justify  the  use  of  the  more  complex 
nonlinear  systems.  Darly  progress  in  this  study  includes  the  derivation  of  an 
optimum  (minimum-error- probability)  receiver  structure  for  a  two-user  channel 
and  the  development  of  simplifyi-ng  approximations  to  this  optimum  receiver 
that  will  facilitate  the  desired  performance  analysis. 

15.2  Signal  Deteotion,  Estimation,  and  Filtering* 

Progress  in  this  general  area  has  involved  two  primary  continuing  areas 
of  study  as  well  as  the  initiation  of  a  new  area.  The  two  continuing  topics 
of  study  are:  (l  )  robustness  in  detection,  estimation,  and  filtering,  and 
(2)  the  design  and  analysis  of  digital  detection  systems.  Progress  in  chess 
two  areas  is  described  in  the  folio-wing  subsections.  The  new  area  of  sfudy  is 
signal  detection  in  multiple-access  channels,  and  early  progress  in  chis  study 
is  reported  briefly  in  Section  15.1.6. 

15.2.1  Robustness  _i^  Detection,  Bstimation,  and  Filtering 

In  general,  the  area  of  robustness  deals  -with  the  design  of  systems  and 
procedures  which  are  relatively  insensitive  (in  terms  of  performance)  to  small 
deviations  from  an  assumed  model.  Robust  techniques  are  of  interest  in  a  wide 

*  This  research  was  supnorted  by  the  Army  Research  Office  under  contract 
DAAD29-31  -<-0062,  the  foint  Ser'/ices  Electronics  Program  under  contract 
'100011-79-0-0124,  the  Office  of  Naval  Research  under  contract 
N0001 1-31 -0-001 1,  and  the  National  Science  Foundation  'under  grant 
ECS--^9-1  6155. 
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variety  of  applioations  including  radar,  sonar,  and  seisxology,  inasmuch  as 
Lnaccurats  models  are  the  norm  for  such  situations;  and  robustness  has  been 
studied  videly  in  the  contexts  of  communications,  control,  and  inferential 
statistics.  ITne  basic  assumption  for  the  analytical  study  of  robustness  in 
detection  and  filtering  is  that  the  statistics  of  signals  and/or  noise  are  not 
'down  exactly,  but  rather  are  known  to  be  within  some  (usually  nonparametric) 
classes  representing  uncertainty  in  the  'underlying  model.  Tne  primary 
technique  for  designing  a  robust  system  is  to  seek  a  system  achieving  the  best 
worst-case  performance  over  the  relevant  uncertainty  classes;  i.e.,  the 
primary  design  philosophy  is  minimax.  Within  this  context  several  detection 
and  estimation  problems  have  been  considered  'under  this  project,  and  these  are 
described  in  the  following  paragraphs. 

The  problem  of  robust  linear  smoothing  of  a  stationary  random  signal  with 
uncertain  speotr'um  observed  in  additive  noise  -with  'uncertain  spectr'um  is 

considered  in  ^5.21  1.  Here,  a  general  solution  to  this  problem  is  given  for 
spectral  'uncertainty  classes  of  a  general  type  based  on  Tnoquet  capacities. 

This  type  of  model  includes  standard  -uncertainty  models  such  as  contaminated 
mixt'ures  as  well  as  several  topological  models  of  'uncertainty.  The  'usefulness 
of  these  results  is  extended  in  ^I5.5‘5]  by  generalizing  the  class  of  models 
that  can  be  treated  in  the  framework  of  [15.21].  "doreover,  an  extensive 
n-umerical  analysis  [lo.oo]  indicates  that  the  worst-case  performance  of  the 
proposed  techniques  is  generally  much  better  than  that  resulting  from  designs 
which  simply  ignore  the  presence  of  'uncertainty.  More  general  problems  of 
rob'ust  estimation  of  stationary  signals  (incl'uiing  smoothing,  filtering,  ani 
prediction)  also  have  been  considered  for  the  case  of  discrete  time  Jo-o"^]. 
Here,  a  general  minimax  res'ult  is  given  from  -which  robust  solutions  to  a 

variety  of  problems  follow  straightforwardly. 

Also  considered  in  this  general  context  is  the  problem  of  state 
estimation  in  linear  stochastic  systems  with  -uncertainties.  T-;wo  approaches 
have  been  considered  for  this  problem.  One  of  them  is  to  consiier  -cinimax 
designs  for  systems  -with  uncertain  state  and  observation  noise  statistics 
[15.23,15.31].  Here,  under  the  ass'umption  that  the  noise  is  -white  -with 

uncertain  componentwise  correlation,  it  is  sho-wn  that  the  minimax- 
mean- sq’uare-error  state  estimator  is  the  Malman  filter  for  a  least-favorable 
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xolel,  ^nd  several  examples  are  considered.  delated  results  for  xinimax 
linear-quadratic  control  are  found  in  ■_15.4^]  and  are  liscussed  in  the 
Decision  and  Control  'unit  of  this  report.  Another  approach  to  state 
estimation  with  'uncertainties  concerns  systems  'vith  random  jumps  '/rhich  are 
modeled  by  considering  the  standard  linear  system  model  -vith  adiitional  marked 
Poisson  input  noise.  1  suboptimal  sequential  smoothing  algorithm  for  this 
model  is  developed  and  analyzed  in  [15.5^.  This  scheme  involves  detecting  the 
incidence  times  and  estimating  the  marks  of  the  Poisson  input  process  and  then 
reconstr’ucting  the  state  -vith  a  (conditional)  Kalman  smoother.  The 
performance  of  this  scheme  is  seen  to  be  s’uperior  to  that  of  optimum  causal 
estimators  and  linear  nonca-usal  estimators,  and  the  implementation  of  the 
scheme  is  superior  to  the  optim'um  (nonlinear)  smoother,  -vhich  is  described  by 
an  infinite  sequence  of  filtering  equations.  .1  related  st'udy  is  found  in 
Jo. 5'’]  '-riiich  considers  the  problems  of  minixax  (or  robust)  filtering  and 
state  estimation  for  obser/ed  nonhomogeneo'us  Poisson  processes  vith  -uncertain 
rate  f'unctions.  In  this  situation,  analogies  '.vith  the  case  of  contin'uous 
observations  can  be  exploited  to  derive  minimax  designs. 

Ceveral  topics  in  rob'ust  signal  detection  have  also  been  considered. 

Three  basic  problems  have  been  treated  in  this  context.  Pirst,  the  role  of  a 
generalised  measure  of  signal- to-noise  ratio  ('xno’rfn  as  the  deflection)  in 
designing  robust  detection  systems  is  considered  in  (15.25].  It  is 

demonstrated  that  this  criterion  is  often  much  simpler  than  earlier  criteria 
for  use  in  robust  design  because  a  constructive  method  for  finding  robust 
solutions  is  available  for  the  new  criterion  -whereas  no  such  general  method  is 
available  for  earlier  formulations,  A  second  problem  involves  the  design  of 
rob’ust  systems  for  detecting  signals  in  the  presence  of  dependent  noise 

J5.25].  The  proposed  technique  uses  a  moving- average  modal  for  noise 
dependence  -and  considers  designs  that  are  minimax  robust  to  first  order  in  the 
degree  of  dependence.  The  third  effort  in  this  area  involves  finding 

solutions  to  the  problem  of  robust  discrete- time  matched  filtering  J5.55]. 
The  problem  of  robust  matched  filtering  has  been  considered  previously  -within 
a  general  Hilbert  space  formulation  ^15.2'^],  and  J5.53i  derives  specifi: 
sol'itions  for  the  particular  problem  of  discrete- time  matched  filtering  -rfithin 
a  variety  of  -unoertainty  models  for  signals  and  noise. 
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15.2.2  Digital  Detection  Syetema 

The  second  continuing  area  of  atudy  in  ietection,  estimation,  and 
filtering  involves  the  analysis  and  optimum  design  of  digital  signal-detection 
systems.  This  problem  is  of  interest  in  a  -^file  variety  of  applications  in 
v/hich  receiver  digitization  is  increasingly  important  for  efficient 
implementation.  Recent  progress  in  this  general  area  can  be  divided  roughly 
into  two  main  topics:  optimum  data  quantization  for  detection  systems  and 
optimum  nonlinear  designs  for  signal  ietection  in  dependent  noise. 

Progress  in  the  first  of  these  areas  is  reported  in  M5.4T,  J5.29j,  and 
']l5.50].  In  '.15.4I,  results  of  an  earlier  study  on  optimum  data  quantization 
for  coherent  detection  systems  are  extended  to  systems  for  the  ietection  of 
purely  stochastic  signals.  Optimum  and  suboptimum  schemes  are  compared,  and 
it  is  seen  that  significant  detection  efficiency  can  be  gained  by  implementing 
in  optim’um  quantization  scheme  rather  than,  for  example,  by  using  standard 
minimum- disto rtion  quantization  in  a  quadratic  receiver.  These  results  are  of 
Lntsrest  in  sonar,  radio  astronomy,  and  seismological  applications  in  which 
signals  are  often  modeled  as  nondeterministic  stochastic  processes.  A  further 
study  of  quantization,  ;15.50],  considers  the  design  and  analysis  of 
quantizers  for  use  in  memoryless  (recursive)  signal  detection  systems  for 
detection  in  m-iepenient  noise.  It  is  demonstrated  here  that,  by  utilizing 
the  dependence  structure  among  noise  samples,  one  can  improve  memoryless 
ietection  performance  considerably  without  increasing  complexity.  Moreover  it 
is  demonstrated  that,  in  the  presence  of  strong  dependence,  quantizer- 
detectors  utilizing  even  small  numbers  of  quantization  levels  can  be  nearly  as 
efficient  as  the  optimam  (Neyman-'^earson)  detector  for  the  case  of  Daussian 
noise.  A  third  problem  considered  is  that  of  designing  optim’um  input 
amplitude  compressors  for  coherent  detection  systems  utilizing  uniform 
quantization.  It  is  shown  in  ^15.29]  that  the  optimum  such  input  compressor 
is  a  scaled  version  of  the  locally  optimum  detection  nonlinearity,  and  that 
the  performsince  degradation  (over  locally  optim’um  detection)  resulting  from 
the  proposed  quantization  scheme  is  negligible  for  only  a  moderate  number  of 
quantization  bits. 
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Procaiarea  for  lasi^iiiag  effloiant  systama  for  ietaoting  signals  in 
iepenlsnt  noise  bacicgrounda  have  been  liscussed  in  two  of  the  above  oontexta, 
notably  [lo.2o]  -rfhich  ieala  with  robust  detection  in  weakly  dependent  noise 
and  [l5.50]  'which  deals  with  optixim  quantization  for  detection  in  m-dependent 
noise.  Optini’jai  detection  for  the  -wealc-dependence  x.odel  of  J5.26]  is 
considered  in  [15.25].  It  is  shown  here  that,  for  a  aiovin®- average  dependence 
lodel ,  efficient  detection  can  be  achieved  by  reshaping  the  independent-noise 
optim-oni  detection  nonlinearity  -with  an  additive  linear  correction  tern.  in 
investigation  of  several  commonly  used  noise  models  indicates  that  the 
resulting  performance  improvement  is  most  significant  for  impulsive  types  of 
noise.  In  [l5.46]  similar  modified  structures  are  considered  for  detecting 
signals  in  4>-mixing  noise.  It  is  demonstrated  that  the  il>-mixing  model  with 
only  the  mixing  coefficients  given  is  not  sufficiently  descriptive  of 
dependence  structure  to  admit  a  design  which  is  uniformly  better  than  the 
independent-noise  design  (ihis  is  in  contrast  to  the  moving-average  model  of 
[i5.25l).  Ho  wsver,  optimum  memoryless  designs  for  p-mixing  and  other 
dependence  models  can  be  derived  using  general  results  of  'J5.22]  provided 
that  the  second-order  distributions  of  the  noise-process  are  ’xno'wn.  The 
optim'om  system  in  this  case  is  the  solution  to  a  Fredholm  operator  equation  of 
the  second  kind  and,  as  such,  can  be  characterized  in  terms  of  the 

3arrett-Lampard  expansion  of  the  noise  process.  A  further  result  in  this 
general  area  involves  the  comparison  of  memoryless  detection  systems  to 

systems  with  linear  memory  (such  as  the  'leyman-’l’earson  detector  for  the  case 
of  laussian  noise).  In  [15.34]  it  is  3ho:wn  that,  for  a  general  situation 
involving  stationary  noise,  memoryless  systems  are  asymptotically  as  efficient 
as  systems  "with  linear  memory.  This  indicates  that  implementation  complexity 
can  be  reduced  '.without  sacrificing  performance  for  such  situations  involving 
large  sample  sizes. 

15.5  Tata  Compression  Theory  and  Techniques 

lur  research  on  data  compression  has  been  extended  into  several  new 
classes  of  problems  including  the  construction  of  universal  source  codes, 
bounds  on  the  redundancy  of  universal  data  compression  techniques,  and  an 

investigation  of  the  redundancy  of  Huffman  codes.  Existence  results, 

red'undancy  bo'unds,  and  construction  technique.e  for  xniversal  source  codes  for 
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ala3383  of  meitoryleaa  sources  are  given  in  [)5.9].  The  codes  considered  are 
fixed-length- to-variable-length  (?L-VL)  source  codes.  New  results  have  been 
obtained  on  other  types  of  source  coles  including  variable- length- to- 
fixei-length  (7L-FL)  codes.  For  short  block  lengths  we  have  fo’unl  that 
properly  designed  7L-FL  ^universal  source  codes  provide  lower  redundancy  than 
our  best  FL-VL  source  codes  previously  reported  in  [l5.9l. 

Our  investigation  of  FL-7L  universal  source  codes  for  'unifilar  Narkov 
sources  has  provided  several  new  results.  A  computational  formula  was  found 
for  evaluating  upper  and  lower  bounds  on  the  redundancy  of  weighted  ’oniversal 
coding  for  a  binary  two-state  Narkov  source,  where  the  weighting  distribution 
of  the  transition  probabilities  is  'oniform.  This  formula  was  also  analyzed 
for  large  block  lengths  and  found  to  agree  asymptotically  with  the  upper  bound 
on  the  red'jndancy  for  minimax  coding.  Farlier  results  on  upper  bounds  on  the 
red'ondancy  for  minimax  coding  for  unifilar  sources  were  extended  to  cases 
where  the  structure  of  the  unifilar  source  was  'unkno-wn. 

Another  problem  •.'ihich  has  been  investigated  is  that  of  the  redundancy  of 
Huffman  coding.  Huffman's  algorithm  is  a  well-known  procedure  for 
constructing  a  variable-length  code  of  minim’jm  possible  average  length  for  a 
source  -with  precisely  'ioiown  statistics.  The  red’undancy  of  a  Huffman  code  is 
defined  to  be  the  difference  between  its  average  length  and  the  source 
entropy.  Ne  have  developed  a  powerful  general  technique  for  estimating  this 
red'undancy.  It  implies  for  example  that  most  sources  have  redundancy  near  the 
value  logjlog^  e"'  *  .025766575,  and  that  apart  from  a  few  exceptional 

cases,  no  code  has  redundancy  exceeding  1  -  log,  e  ^  log,  e  log,  e  » 

.085071552.  This  work  is  reported  in  ^15. 48]. 
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16.1  Computer  Algorithms  for  Integrated  Circuit  Design* 

In  this  research  a  study  was  carried  out  on  the  use  of  partitioning  and 

tearing  algorithms  for  the  analysis  of  large-scale  circuits  and  systems.  In 

this  approach  the  circuit  is  partitioned  into  subcircuits,  each  subcircuit  is 

solved  separately,  and  then  the  subcircuit  solutions  are'  combined  to  obtain 

the  solution  of  the  entire  circuit.  Tearing  provides  many  computational 

advantages,  especially  in  the  simulation  of  large-scale  integrated  circuits. 

First,  computer  storage  can  be  reduced  if  identical  subcircuits  are  created. 

Secondly,  computer  time  can  be  saved  by  exploiting  subcircuit  latency;  that 

is,  subcircuits  that  are  inactive  during  a  given  time  intar/al  can  easily  be 

by-passed  in  the  solution  process.  Thirdly,  the  use  of  tearing  methods  allows 

*  This  worlc  -iraa  supported  by  the  Joint  Services  Electronics  Program 
('J.3.  Army,  U.3.  Navy,  ‘  and  U.3.  Air  Force)  under  Contract  N'DOOl  A-'^9-'S-'D12l, 
and  by  the  National  Science  Foundation  under  Grant  2CS  SO-D'^O'^S. 
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for  the  parallel  processing  of  subcircuit  equations.  Finally,  tearing 
facilitates  the  use  of  reduced-order  or  macromodels  in  the  analysis  and 
permits  the  application  of  mixed-mode  simulation. 

In  our  research,  equation  solution  techniques  for  the  analysis  of 
circuits  and  systems  were  investigated.  In  particular,  when  tearing  methods 

are  employed,  there  are  a  number  of  solution  techniques  that  can  be  used  to 

exploit  latency.  These  techniques  have  been  studied  and  for  electronic 

circuits  the  most  efficient  solution  technique  has  been  identified  [l6.'5S]. 
In  conjunction  with  this  work  a  fast  and  reliable  ordering  technique  for  the 
modified  nodal  aquation  formulation  was  found  [16.6].  Our  ordering  procedure 
eliminates  the  time-consuming  pivot  selection  process  that  is  usually  required 
for  the  solution  of  the  modified  nodal  equations.  In  addition  to  the  above 
work  on  the  solution  methods,  various  latency  criteria  for  the  transient 
analysis  of  digital  circuits  'were  studied.  In  order  to  test  our  algorithms 

the  3PICS2  circuit  analysis  program  -ms  modified  to  include  tearing  methods. 
This  new  circuit  analysis  program  is  called  SLATS  [16.31];  in  general,  this 
program  runs  faster  than  the  3PICB2  program  for  the  transient  analysis  of 
digital  logic  circuits  and  yields  the  same  accuracy. 

In  related  wo'rk,  algorithms  for  efficient  solution  updating  when  large 
changes  in  the  values  of  some  system  parameters  occur,  have  been  derived 
[15.4.].  This  problem  has  also  been  investigated  when  changes  occur  in  large, 
partitioned  systems  [l6.22]. 

In  addition,  the  problem  of  finding  the  shortest  path  in  a  network,  which 
occurs  in  network  routing  among  other  applications,  has  been  studied  and  an 
efficient  updating  algorithm  when  changes  occur  in  the  network  parameters  has 
been  derived  [16.211. 

In  many  cases  electronic  device  models  are  represented  in  tabular  form. 
This  tabular  form  can  be  used  to  approximate  the  nonlinearities  by 
piecewise-linear  functions.  In  our  research  the  properties  of  the  solutions 
of  piecewise-linear  circuits  have  been  investigated  in  detail  [ 1 6. 5, 1 5. 20] . 

Research  has  also  been  done  on  the  use  of  piecewise  nonlinear  models. 
The  use  of  piecewise  nonlinear  models  for  the  exponential  nonlinearity  of  a  pn 
junction  offers  great  flexibility  in  modeling  [  1  6. 29, 1 6.  However  these 

models  do  not  satisfy  the  requirements  of  the  standard  iterative  solution 
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techniques.  In  this  work  we  have  explored  new  solution  algorithms  tailored  to 
the  piecewise  nonlinear  models  and  compared,  computationally,  the  new 
algorithms  and  models  to  many  of  the  traditional  methods  used  in  programs  for 
the  calculation  of  the  DC  operating  point  of  bipolar  circuits.  Our  results 
indicate  that  the  piecewise  nonlinear  models  and  algorithms  compare  very  well 
to  the  best  methods  that  are  available  today  6.30, 1 5.36]. 

Future  work  in  this  area  will  include  theoretical  studies  of  the 
algorithms  for  piecewise  nonlinear  models.  Also  the  applicability  of  this 
t3rpe  of  modeling  for  use  in  MOS  circuits  will  be  studied. 

1 5.2  ^  New  Structured  Logic  A rray* 

The  use  of  highly  structured  architectures  is  becoming  necessary  to  take 
advantage  of  the  advent  of  VLSI  circuits.  One  area  that  lends  itself  to 
structuring  is  that  of  logic  arrays.  These  arrays,  such  as  ROM's,  PLA's,  gate 
arrays,  and  custom  logic  cells,  are'  widely  used  throughout  the  IC  industry. 
The  success  of  these  structured  logic  arrays  is  in  large  part  due  to  the 
development  of  successful  GAD  tools  for  their  design,  layout,  and  routing. 

There  are  two  broad  classes  of  structured  arrays:  customer 

reprogrammable  and  customer  nonreprogrammable .  In  this  work  [15.34]  we  are 
investigating  a  new  architecture  for  a  customer  reprogrammable  logic  array. 
This  array  structure,  known  as  REPLICA  (Random,  Electronically  Programmable, 
Logic  Integrated  Custom  Arrayl  is  a  regular  structure  of  NOR  gates  with 
restricted  fan-in  and  fan-out.  We  have  developed  CAD  tools  for  the 
minimization  of  logic  realizations  in  REPLICA  similar  to  the  folding 
algorithms  used  for  PLA's.  We  have  also  developed  guidelines  for  the  optimal 
form  for  the  logic  equations  to  be  Implemented  in  REPLICA.  Several  examples 
of  logic  functions  implemented  in  REPLICA  indicate  that  it  competes  well  with 
PLA's  and  has  some  features  that  may  make  it  more  iesireable  than  PLA's  in 
situations  such  as  realizations  of  sequential  logic.  Work  is  progressing  on 
comparisons  of  REPLICA  with  existing  structured  arrays  and  on  improving  the 
CAD  tools  necessary  for  realizing  larger  teat  examples  in  the  REPLICA 
architecture . 

♦  This  work  was  supported  by  the  Joint  Ser'/ices  Electronics  Program 
'■J.3.  Army,  'J.3.  Navy,  ’  and  U.3.  Air  Forced  under  Contract  N0001 4-"’?-C-0d24. 
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1 0 . 3  Hierarchical  Design  Techniques  for  VLSI  Systems* 

With  the  complexity  of  the  systems  that  can  be  realized  on  a  single  VLSI 
chip  the  design  time  is,  using  conventional  techniques,  on  the  order  of  100 
man-years.  In  order  to  reduce  this  unacceptable  design  time,  structured 
design  techniques  are  being  investigated.  In  this  work  [l6.37l  the  impact  of 
a  structured  design  methodology  on  the  design  of  a  VLSI  system  was  studied. 
The  design  and  testing  of  the  chip  were  both  considered  from  the  beginning  of 
the  design  process.  As  a  test  example  a  CORDIC  function  (used  for  realizing 
trigonometric  functions)  was  designed.  The  function  to  be  implemented  was 
studied  by  simulation  in  order  to  finalize  the  architecture  of  the  design. 
Then  followed  a  structural  design,  a  physical  design  including  layout,  and  a 
verification  using  different  simulation  tools. 

The  investigation  clearly  showed  that,  with  some  sacrifice  in  chip  area, 
a  structured  design  approach  was  feasible  for  VLSI.  Various  other  systems  are 
now  being  studied  for  implementation  in  VLSI  using  structured  design 
techniques.  The  goal  of  these  studies  is  the  refinement  of  the  tools  and 
techniques  necessary  in  the  design  hierarcy  for  VLSI  systems. 

16.4  The  Use  of  Multiple  Criterion  Optimization  for  the  Design  of 
Digital  Filters** 

The  design  of  digital  filters,  both  FIR  [l  6. 1 2 , 1 6. 32 1  and  HR,  using 
competing  criteria  [16.3]  is  being  studied.  In  particular,  various 
formulations  of  the  design  problem  which  trade-off  between  magnitude  and  phase 
performance  are  being  considered.  These  formulations  involve  the 
specification  of  a  desired  phase  behavior  as  well  as  the  standard  desired 
magnitude  behavior. 

Performance  constraints,  such  as  monotonic  magnitude  response  in  the 
transition  regions,  are  also  Imposed  on  the  behavior  of  the  filter.  In  order 
to  solve  these  rather  difficult  optimization  problems  we  are  modifying  the 
sequential  quadratic  method  for  constrained  optimization  proposed  by  Han  and 

*  This  work  was  supported  by  the  Joint  Services  Flectronics  Program 
(U.S.  Army,  U.3.  Navy,  and  U.3.  Air  Force)  under  Contract  N0001 4-''9-C-0424. 

**  This  work  was  supported  by  the  Joint  Ser/ices  Electronics  Program 
('J.3.  Army,  U.S.  Navy,  and  U.S.  Air  Force)  under  Contract  N0001 4-’9-C-0421, 
and  by  the  National  ocience  Foundation  under  Grant  ECS  30-070"'5. 
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Powell.  ¥e  are  also  exploring  the  use  of  a  dsmamic  frequency  grid  which  will 
allow,  theoretically,  convergence  to  the  true  solution  of  the  design  problem 
and  not  the  approximation  that  is  found  in  all  existing  digital  filter  design 
programs  today. 

15.5  Analytical  and  Computer-Aided  Design  of  New  Circuit  Structures  for 

Switched  Capaoitor  Filters* 

Research  is  being  conducted  in  optimum  switched  capacitor  (SC)  filter 
structures  and  in  computer-aided  analysis  and  design  (CAD)  of  SC  circuits. 
Recent  studies  of  optimum  SC  filter  structures  resulted  in  two  different 
designs  using  the  bilinear  z- transformation.  One  design  uses  a  state  space 
technique  [l6. 1 1 , to. 14, 16. 1?]  while  the  other  uses  a  follow- the-leader 
feedback  (FL?)  topology  [ib.lO].  Both  techniques  allow  SC  filters  to  be 
implemented  'with  stray  insensitive  first  order  modules  (such  as  forward  Buler, 
backward  Suler,  or  bilinear  descrete  integrators)  or  second  order  biquadratic 
modules.  Optimization  procedures  were  applied  to  maximize  dynamic  range  while 
minimizing  sensitivity,  capacitance  ratios,  and  total  capacitance.  Impressive 
improvements  in  sensitivity,  as  analyzed  by  Honte  Carlo  analyses,  have  been 
obtained  for  an  example  of  an  optimized  fifth-order  Chebychev  LP  filter 
designed  with  the  state  space  and  FL?  approaches. 

In  the  area  of  CAD,  a  user-oriented  program,  SCAP  II,  based  on  modified 
nodal  analysis,  was  completed.  Frequency  and  sensitivity  analyses  were 
demonstrated  on  a  wide  variety  of  useful  SC  circuits  taken  from  the  literature 
[16.25,16.26,16.55].  SCAP  II  is  capable  of  performing  frequency  analysis  and 
frequency  domain  sensitivity  analysis  on  SC  circuits  that  contain  M03 
switches,  capacitors,  and  voltage  controlled  voltage  sources,  and  which 
operate  with  a  l-phase  50'^  duty  cycle  switching  sequence.  The  analysis  is 
efficiently  implemented  with  sparse  matrix  numerical  methods  and  a  specially 
designed  aquation  ordering  strategy  which  minimizes  fill-in  during  computation 
of  the  LU  factors.  Sensitivity  analysis,  -iriiich  is  an  implementation  of  first 
order  relations  obtained  by  differentiating  the  original  modified  nodal 
equations,  is  efficiently  implemented  for  capacitors,  capacitor  ratios,  and  op 

*  This  work  is  supported  by  the  National  Science  Foundation  'jnler  Grant 

3^0-73-1^315. 
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amp  gains  by  using  the  same  LU  factors  as  already  computed  for  the  frequency 
analysis.  3CAP  II  requires  approximately  twice  the  machine  time  as  the 
previously  reported  block-partitioning  technique  of  oCAP  I,  and  hence  appears 
to  be  a  useful  aid  in  SC  circuit  design  [  1  6. 1 5. 23] . 

The  most  recent  activity  in  SC  CAD  has  been  the  development  of  a  very 
general  N-phase  SC  analysis  routine,  referred  to  as  SCAPN  [l6.13j.  In 
addition  to  using  sparse  matrix  methods  and  proper  reordering  techniques  SCAPN 
uses  a  specially  designed  "pre-LU  factorization  algorithm"  which  minimizes  the 
amount  of  computation  that  must  be  done  to  update  the  LU  factors  at  each  new 
frequency.  SCAPN  also  features  a  very  efficient  sensitivity  analysis  based  on 
the  adjoint  network  concept  [lo.ld],  and  a  nonlinear  distortion  analysis  based 
on  the  Volterra  series  technique  for  analyzing  the  effects  of  driving  MOS 
capacitors  and  amplifiers  into  nonlinear  operation  ['5. 15]. 

16.5  Automatic  Tuning  Algorithms  for  Analog  Filters* 

To  meet  filter  response  specifications  analog  filters  usually  must  be 
tuned  or  adjusted,  preferably  by  computer  automation  if  the  production  level 
is  high.  In  this  research  three  generalized  tuning  algorithms  have  been 
studied  on  the  basis  of  their  architecture,  computational  complexity,  and 
effectiveness,  further,  a  method  has  been  found  for  the  tuning  element  and 
frequency  selection  problem,  a  problem  relevant  to  all  three  methods.  Monte 
Carlo  simulations  have  been  done  on  several  active  filter  circuits  in  order  to 
enhance  the  comparison  and  provide  a  demonstration  [l 5.231. 

15.7  Fault  Isolation  in  Analog  Circuits** 

The  objective  of  this  research  is  to  be  able  to  isolate  faulty  components 
in  analog  circuits  from  test  point  measurements.  Presently  the  only  approach 
that  has  had  limited  success  in  practice  is  the  fault  dictionary  method. 

*  This  work  was  supported  by  the  National  Science  Foundation  'under  Irant  3NG 
78-1 1788. 

**Thi3  work  was  supported  by  the  U.3.  Naval  Air  Fngineering  Center  'under 
Contract  N68338-30-C'-'C506  and  by  the  Joint  Services  Slectronics  Program 
(U.S.  Army,  U.3.  Navy,  and  U.S.  Air  Force''  under  Contract  N0001  d-’b— C-012i. 
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Unfortunately  this  method  requires  that  one  anticipate  all  possible  fault 
conditions  in  advance.  This  is  usually  an  impossible  task,  and  reliable  field 
data  is  required  in  order  to  update  the  dictionary.  Typically  it  is  also 
necessary  to  confine  the  dictionary  to  single  catastrophic  faults  in  order  to 
avoid  excessive  stora«;e  requirements  and  Ion?  search  times  for  the  fault 
conditions.  Finally,  in  the  fault  dictionary  it  is  not  easy  to  determine  the 
testability  of  the  circuit  from  a  given  set  of  test  points. 

In  our  work  we  have  shown  that  a  sensitivity  matrix  can  be  used  to 
isolate  both  hard  and  soft  faults.  The  idea  behind  the  method  is  that  if  the 
columns  of  the  sensitivity  matrix  are  linearly  independent,  then  any  fault  in 
a  component  will  cause  a  unique  shift  in  the  test  point  voltages.  The  method 
can  locate  single  as  well  as  multiple  faults  in  analog  circuits.  Furthermore, 
the  independence  of  the  colimns  of  the  sensitivity  matrix  are  a  good  measure 
of  the  testability  of  the  circuit  -with  respect  to  the  given  test  points.  The 
sensitivity  matrix  is  computed  for  the  linearized  circuit  model  at  the  nominal 
parameter  values  in  the  pre-test  phase.  Thus,  only  a  few  simple  mathematical 
operations  are  required  in  the  test  equipment.  A  sensitivity  matrix  must  be 
computed  for  each  test  signal,  and  assuming  that  the  rank  of  the  sensitivity 
matrix  is  equal  to  the  number  of  test  points  for  any  column  ordering,  then  the 
number  of  faults  in  the  circuit  must  be  less  than  the  number  of  test  points  in 
order  for  the  method  to  give  reliable  results. 
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17.1  Control  and  Decision  Strategies  for  Systems  Under  Imperfect  Information* 

Uncertainties  arise  because  of  unknown  system  parameters,  unknown  signal 
environments,  and  hardware  tolerances.  Related  eomplexities  arise  in 
situations  involving  multiple  system  performance  criteria  and  multiple 
decision  makers.  The  objective  of  this  project  is  to  gain  a  basic 
understanding  of  the  behavior  and  control  of  complex  systems.  During  the  past 
year,  attention  has  been  given  to  several  topics.  Highlights  are  summarised 
in  the  following. 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
(U.3.  Army,  U.3.  Navy,  and  U.3.  Air  Force)  ’under  Contract  N0001 1-79-0-0421. 
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Several  results  ooncerning  the  sensitivity  ani  robustness  of  control 
systems  to  'uncertainties  in  the  plant  model  have  been  obtained.  We  have 
ieveloped  bounds  on  the  magnitude  of  allowed  perturbations  which  guarantee 
that  the  perturbed  system  performs  better  than  an  equivalent  open  loop  system. 
These  bounds  are  expressed  in  terms  of  the  singular  values  of  the  return 
difference  matrix.  We  have  also  developed  the  singular  value  sensitivity 
function  as  a  robustness  analysis  tool.  Singular  value  sensitivities  can  be 
used  to  analyze  perturbations  'such  as  multiplying  incident  perturbations  or 
strucrured  perturbations)  for  -which  the  guaranteed  bounds  do  not  apply. 

In  the  area  of  f-uniamental  system  theory  basic  to  control  strategy 
development,  three  general  topics  were  the  subject  of  investigation  this  past 
year:  (l  )  the  input-output  behavior  of  some  nonlinear  feedback  loops,  (2)  the 

invertibility  of  input-output  systems,  and  (3)  the  existence  and  properties  of 
"strange  attractors"  in  system  theory.  We  give  below  a  brief  account  of  these 
studies,  concentrating  on  the  last  one,  since  it  is  new. 

1 .  We  have  treated  the  "locking"  of  a  class  of  nonlinear  feedback  loops 
to  a  certain  set  of  input  signals;  sufficient  conditions  for  the 
locking  have  been  obtained.  Riurther,  sufficient  conditions  for  the 
"tracking"  of  the  input  signal  by  the  output  have  been  found.  These 
results  are  being  extended;  a  manuscript  is  in  preparation. 

2.  liven  an  output,  what  was  the  input  to  the  system?  It  is  this  -which 
constit-utes  the  question  of  invertibility.  Some  partial  results  have 
been  achieved,  in  the  sense  that  sufficient  conditions  for  the 
invertibility  of  a  class  of  input-output  systems  have  been  found. 
These  conditions  are  fairly  straight-forward  extensions  of  previous 
-work.  Currently,  efforts  are  being  mads  to  generalize  the  conditions 
so  that  they  are  applicable  to  a  wide  class  of  systems. 

3.  A  variety  of  systems  are  s-uspected  of  exhibiting  "strange"  behavior. 
3y  this  is  meant  something  of  the  follo'wing  sort  (the  definitions 
varying  according  to  the  investigator);  There  exist  trajectories  of 
the  system  which,  starting  close  to  each  other,  diverge  exponentially 
initially;  these  trajectories  remain  bounded,  however,  for  all  time. 
?inally,  the  trajectories  approach  a  set  of  points  ^  the  "strange 
attractor"  1  -which  are  not  periodic  orbits  or  a  finite  set  of  singular 
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points . 

Tha  force!  van  ier  Pol  equation  with  very  large  damping  and 
proportionately  large  sinusoidal  forcing  and  similar  such  equations 
are  the  only  differential  equations  for  which  it  has  been  proven  that 
a  strange  attractor  exists.  Numerical  solutions  of  the  Lorenz  system 
and  of  a  version  of  the  forced  Duffing  equation  exhibit  behavior  much 
like  that  of  strange  attractors,  but  proofs  of  this  are  absent. 

It  has  been  suggested  that  some  of  the  properties  of  turbulent 
mixing  and  of  a  number  of  other  phenomena  could  be  explained  if 
strange  attractors  existed  for  the  governing  equations. 

Ne  have  started  to  study  strange  attractors,  by  repeating  some 
of  the  numerical  work  reported  in  the  literature,  and  by  a  review  of 
the  theory.  Ve  were  unable  to  duplicate  the  results  on  the  forced 
Duffing-like  aquation  for  a  long  time  -until  we  -used  an 
Adams-Sashforth  predictor-corrector  integration  routine  for  the 
numerical  solutions  of  the  differential  equation;  no'w  our  results 
agree  with  those  published.  Ve  are  in  the  initial  stage  of  what  may 
be  a  long  term  project. 


Wsll-poselness  of  leader- follower  strategies  obtained  from  singularly 
perturbed  reduced  order  models  has  been  investigated  in  a  manner  similar  to 
our  approach  previously  developed  for  Nash  strategies  for  reduced  order 
models.  Ye  have  derived  order  reduction  proced’ures  -which  lead  to  well-posed 
formulations  of  the  leader-follower  control  of  singularly  pert-urbed  dynamic 
systems.  Ye  have  investigated  the  strong  connection  between  information 
structure  and  well-posedness  of  singularly  perturbed  Nash  games  and  we  have 
identified  -why  the  natural  order  red-action  is  ill- posed  as  -we  reported 
earlier.  Ither  structures  which  preseir/ed  information  were  demonstrated  to 
lead  to  -well-poaed  solutions.  Details  are  described  in  the  journal  and 
meeting  papers. 

The  problem  of  controlling  a  linear  stochastic  system  so  as  to  minimize  a 
quadratic  steady  state  cost  functional  with  uncertain  process  and  obser/ation 
noise  statistics  was  investigated.  The  problem  was  formulated  as  a  minimax 
optimization.  Ye  showed  that  the  minimax  solution  is  obtained  by  the  usual 


104 


17.  DECISION  AND  CONTROL 


linear  quadratic  deterministic  ?ain  !aultipl3ring  the  state  estimate  provided  by 
a  linear  filter.  The  linear  filter  is  a  Kalman-Busy  filter  for  a  particular 
worst  case  pair  of  spectral  density  matrices.  These  results  parallel  those  of 
the  minimax  state  estimation  problem  discussed  in  Area  15. 

17.2  Implementation  Constrained  Decomposition  and  Hierarchical  Control* 

In  large-scale  systems,  control  taslcs  may  be  decomposed  and  hierarchical 
levels  may  be  imposed  on  the  basis  of  analytically  derived  strategies.  This 
project  is  devoted  to  the  fundamental  problem  of  investigating  an  analytical 
framework  for  incorporating  realistic  constraints  of  a  computer  network. 
Questions  concerning  information  flow,  time  scales,  and  loss  of  feedback  loops 
are  examined  in  coordinating  the  multiple  controllers  in  a  network. 

We  have  investigated  and  developed  several  classes  of  hierarchical 
optimization  algorithms.  Glasses  of  algorithms  which  exploit  a  weak 
structural  coupling  between  variables  of  an  optimization  problem  were  analyzed 
using  local  convergence  theory.  Natural  decompositions  were  defined  in  terms 
of  splittings  which  directly  exploit  the  structure  of  the  problem  in  the 
decomposition.  Vs  demonstrated  that  iterations  defined  by  natural 
decompositions  would  converge  whenever  the  coupling  is  sufficiently  weak.  Ve 
also  have  developed  a  hierarchical  optimization  algorithm  which  can  be  used  to 
solve  the  linear  quadratic  robust  control  problem.  The  algorithm  exploits  the 
structure  of  the  dual  problem  to  define  a  simple  hierarchical  algorithm  which 
does  not  require  a  stabilizing  initial  point  and  which  has  guaranteed 
convergence  properties. 

A  common  feature  in  large  scale  system  practice  is  that  different 

decision  makers  assume  different  simplified  low-order  models  of  the  same  large 

scale  system.  This  may  be  due  to  the  necessity  of  easing  the  computational 

burden  in  evaluating  strategies  or  the  lack  of  adequately  modeled  dynamics  of 

some  parts  of  the  large  scale  system.  Based  on  these  simplified  models,  the 

decision  makers  would  then  attempt  to  arrive  at  low-order  decentralized 

strategies  which  would  be  economically  feasible  to  implement,  at  the  same  time 

preser/ing  near-optimality.  The  problem  of  designing  control  strategies  for 

decision  makers  under  a  situation  where  they  have  different  models,  different 

*  This  work  was  supported  by  the  Joint  oervices  Electronics  Program 
(H.3.  Army,  U.3.  Navy,  and  U.3.  Air  Force)  'onder  Contract  N0001 4.-79-C-0424'. 
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infomation  seta,  and  diffarent  objective  functionals,  has  been  addressed 
previously  within  the  framework  of  multiparameter  singular  perturbations.  In 
this  framework,  a  large  scale  system  is  viewed  as  consisting  of  a  "slow"  core 
coupled  to  a  number  of  "fast"  subsystems.  A  multimodel  situation  results  when 
each  decision  maker  retains  the  dynamics  of  one  fast  subsystem  and  the  slow 
core,  while  neglecting  the  dynamics  of  the  remaining  fast  subsystems.  This 
year  the  general  problem,  wherein  the  fast  subsystems  need  not  be  weakly 
coupled  and  each  fast  subsystem  might  be  affected  by  more  than  one  decision 
maker,  has  been  addressed.  A  procedure  has  been  formulated  to  obtain  near 
optimal  decentralised  Nash  strategies  in  the  two  player  case  under  a 
multimodel  situation  'Within  the  framework  of  multiparameter  singular 
perturbations. 

1  .  3  Control  Strategies  for  Complex  Systems  for  'Jse  in  Aerospace  Avionics* 

’.Vhenever  model  -uncertainty  is  present  or  a  range  of  operating  conditions 
is  anticipated,  engineering  analysis  and  design  must  deal  -with  questions  of 
sensitivity.  Performances  of  manufact'ured  components  -are  necessarily 
specified  -with  nonzero  tolerances.  The  parameter  values  'which  characterize 
these  components  are  therefore  uncertain.  These  parameter  val-ues  may  change 
during  the  operating  lifetimes  of  the  components  due  to  aging  and  due  to 
changes  in  environmental  conditions,  such  as  pressure  and  temperature.  In 
addition,  mathematical  models  used  for  analysis  and  design  of  actual  systems 
cannot  possibly  lead  to  predicted  performance  -which  exactly  matches  the 
performance  of  the  actual  systems.  It  is  useful  to  regard  some  of  the 
parameters  of  the  models  as  uncertain  in  order  to  make  these  models  mors 
realistic.  The  combined  effects  of  parameter  uncertainties  on  overall  system 
behavior  are  of  principal  consideration  in  any  system  design. 

Ve  investigate  uncertainty  using  three  techniques.  Dne  point  of  view 
regards  the  parameters  as  unknown  but  deterministic,  and  studies  the 
uncertainty  situation  -using  parameter  sensitivity  methods.  New  results  here 
include  establisment  of  some  relationships  bet'ween  singular  val-ue  robustness 
approaches  and  the  now-standard  comparison  sensitivity  methods.  Design- 
oriented  techniq-ues  are  being  developed  for  large  parameter  variations  in 
linear  and  nonlinear  systems. 

•  This  -work  was  supported  by  the  7.3.  Air  Force  -mder  Grant  AFD3R-t'3-^6'^3 . 
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A  second  approach  models  the  parameters  as  random  variables.  The  tools 
for  this  approach  are  stochastic  control  and  dual  control  theory,  includini? 
stochastic  multi-controller  situations.  The  interplay  between  information 
structure  and  control  has  been  emphasized  in  our  recent  work. 

We  have  expanded  our  basic  research  to  encompass  an  area  of  increasing 
importance  in  control  research  and  aerospace  applications  -  adaptive  control. 
Our  approach  focuses  on  the  role  of  singular  perturbation  theory  in  the 
development  of  model  reference  adaptive  systems  when  the  plant  and  reference 
model  have  different  orders.  An  important  possible  interaction  has  been 
iiscovered  between  frequencies  present  in  a  "suff iciently  rich"  adaptive  input 
signal  and  neglected  high  frequency  terms  in  a  reduced -order  model  in  a  model 
reference  adaptive  system.  Details  of  these  results  may  be  found  in  the 
journal  articles  and  conference  papers  listed  in  the  publications  section. 

17.4  Large  Scale  Systems* 

Several  new  results  have  been  obtained  in  multimodeling,  that  is,  in 
modeling  and  control  situations  in  which  different  decision  makers  employ 
different  models  of  the  same  system.  Multiparameter  singular  perturbation 
techniques  are  employed  to  capture  the  multimodel  nature  of  fast  dynamic 
subsystems  interconnected  through  slow  variables.  In  one  problem  studied,  a 
decentralized  filtering  and  control  scheme  was  developed  for  obtaining  low 
order  Nash  equilibrium  strategies  in  a  multimodel  setting.  The  decision 
makers  have  decentralized  information  structures  and  are  constrained  to  use 
finite  dimensional  compensators.  In  another  study,  physical  examples  (e.g., 
networks)  are  used  to  reexamine  the  role  of  modeling  in  large  scale  system 
analysis  and  design.  A  general  property  of  the  systems  considered  is  that 
they  are  strongly  coupled  (physically)  in  the  slow  time  scale,  and  weakly 
coupled  in  the  fast  time  scale.  This  leads  to  a  multimodel  situation  in  which 
every  subsystem  controller  can  neglect  all  other  fast  subsystems  except  his 
own. 

In  the  area  of  stochastic  multiplayer  strategies,  an  indirect  method  to 

obtain  the  Stackelberg  solution  of  two-person  nonzero-sum  differential  games 

with  general  dynamics  and  cost  functionals  has  been  developed  with  the  leader 

*  This  work  was  supported  by  the  National  Science  Foundation  under  3rant  ECS 
79-19-596. 
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having  access  to  samplsd  values  of  the  state.  This  new  Tiethod  involves  a 
three-step  optiaization,  and  converts  the  original  dynamic  problem  into  two 
open-loop  optimal  control  problems  and  one  static  finite-dimensional 
optimization  problem.  Solutions  of  these  optimization  problems  are  obtained 
as  solutions  of  a  set  of  differential  equations  with  mixed  bo’undary 
conditions,  which  in  turn  leads  to  the  Stackelberg  cost  value  of  the  leader, 
as  well  as  to  a  class  of  strategies  that  yield  that  cost  level. 

Chained  aggregation,  an  iterative  approach  to  large  scale  system  model 
reduction,  has  been  developed  further.  The  resulting  system  representation, 
the  Ceneralized  Hessenberg  Representation  (CHR),  has  been  shown  to  exhibit 
explicitly  the  system  observability  structure.  Ceometric  (in  the  sense  of 
Wonham)  properties  of  the  CHR  representation  have  been  exposed  and  exploited 
in  the  development  of  a  three -component  decentralized  hierarchical  feedback 
design  technique. 

Finally,  coherency,  a  concept  widely  used  in  power  system  studies,  has 
been  examined  from  an  aggregation  point  of  view.  Specifically,  the 
relationship  among  coherency,  aggregation,  and  observability  has  been 
clarified,  and  a  fundamental  relationship  between  "alow  coherency"  and  weak 
coupling  has  been  established.  Details  of  these  and  other  results  are  given 
in  the  list  of  .journal  articles. 

17.5  Power  System  Normal  and  Security  State  Assessment* 

Normal  and  security  state  assessment  of  large  scale  power  systems 
requires  the  analysis  of  power  system  performance  ijnder  all  levels  of  likely 
load  and  contingencies.  Performance  is  measured  in  terms  of  bus  voltage 
levels,  line  power  flows,  stability  margins,  and  their  sensitivities  to  large 
changes  in  network  conditions.  Explicit  load  flow  solutions  have  been 
leveloped  to  provide  near  closed  form  analysis  of  network  respon.ae  to  total 
load.  These  explicit  forms  may  be  used  to  produce  time  varying  solution 
tra.jectories  and  exact  stochastic  load  flow  solutions. 


*  This  work  was  suoported  by  the  National  science  Foundation  under  Irant  “TH 
75-0559i,  and  the  'Jniversity  of  Illinois  Power  Affiliates  Program. 
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Many  optimal  load  flow  algorithms  require  stability  margin  constraints. 
These  constraints  have  typically  hinged  largely  on  bus  voltage  angle  limits, 
'ifhen  finite  VAR  control  is  included,  these  angle  limits  have  been  aho-<m  to  be 
erroneous  and  in  some  cases  unfeasible.  Extension  of  these  realistic 

stability  margins  to  an  arbitrary  network  is  currently  being  considered. 

17.5  Control  of  Stochastic  Systems  Containing  Parameter  Uncertainty* 

The  main  result  of  this  research  is  the  determination  of  algorithms  for 
Sensitivity  Adaptive  Feedback  with  Estimation  Redistribution  (SAFER)  control. 
This  type  of  control,  which  we  proposed  and  developed,  has  dual  properties  in 
the  sense  that  the  control  takes  into  account  future  parameter  estimation 

accuracy.  The  parameter  sensitivities  are  altered  by  changing  the  sensitivity 
weighting  coefficients  in  a  composite  performance  index.  Our  early  work  on 
this  problem  dealt  with  a  general  multivariable  state  space  formulation  which 
led  to  quite  complicated  algorithms.  We  continually  simplified  the  control 
algorithms  by  first  considering  input-output  representations,  and  finally 
single-input  single-output  representations.  Much  of  the  application  of 
adaptive  control  is  for  this  latter  class  of  problems.  Very  substantial 
reductions  in  computation  are  achieved  not  only  because  of  the  minimal  number 
of  'unknown  parameters  in  the  single- input  single-output  case,  but  also  because 
of  our  use  of  only  two  sensitivity  functions  to  capture  the  impact  of  all  the 

parameters  in  a  SAFER  control  setting.  Our  simpler  algorithm  for  SAFER 

control  has  been  demonstrated  to  yield  performances  which  are  superior  to  the 
certainty  equivalent  control. 

This  project  is  a  cooperative  effort  between  the  Coordinated  Science 
Laboratory  of  the  University  of  Illinois  and  the  Systems  Engineering 
Laboratory  of  the  Venezuelan  Institute  for  Scientific  Research. 


*  This  work  was  supoorted  by  the  National  Science  Foundation  under  Urant  INT 
77-20969. 
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17.7  Aggregated  Hierarohical  Control* 

The  objective  of  this  research  is  to  investigate  hierarchical 
lecomposition  algorithms  for  which  the  problem  solved  at  each  level  of  the 
hierarchy  is  approximated.  Standard  hierarchical  decomposition  theory 
requires  that  an  optimization  problem  be  decomposed  exactly  and  that  the 
iteration  defined  by  the  decomposition  converge  to  the  exact  solution.  Ve 
have  relaxed  the  exactness  requirement  by  allowing  the  exact  problems  at  each 
hierarchical  level  to  be  replaced  by  an  operator  theoretic  aggregate.  The 
hierarchical  structure  is  determined  by  the  composition  of  an  exact 
decomposition  and  an  aggregation.  Conditions  'under  which  the  approximate 
decomposition  iterates  converge  to  an  approximation  of  the  exact  solution  have 
been  determined  for  linear  problems.  The  extension  of  these  results  to 
nonlinear  problems  through  the  use  of  nonlinear  aggregation  and  nonlinear 
splitting  functions  is  under  current  investigation. 


*  This  worh  was  supported  by  the  National  Science 
irNG-71-0377e. 
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13.1  ^  Invest  illation  of  Residue  Nu-aber  Architectures  for  Digital  Filters 

with  Fault  Tolerance* 

This  project  involves  a  study  of  how  the  theory  of  redundant  residue 
number  systems  (RNS)  can  be  developed  and  applied  to  the  isolation  and 
detection  of  hardware  failures  in  digital  processors,  -vi  th  subsequent  "soft 
failure",  that  allows  a  faulty  processor  to  continue  operating  with  reduced 
capabilities.  The  class  of  systematic  redundant  RNS  codes  was  analyzed  and 
simulated,  thereby  establishing  preferred  design  procedures  and  demonstrating 
failure  resistant  behavior  in  recursive  digital  filters  [iS-ipL  At  this 
time,  the  theory  of  failure  resistant  design  for  the  systematic  RNS  codes  is 
relatively  complete,  including  bounds  on  the  necessary  red’undancy  and 
procedures  for  identifying  faulty  modules  and  for  reallocating  resources 
folio-ring  soft  failure.  The  major  issues  now  concern  reliability  of  the  error 
checlcing  hardware  itself,  as  well  as  new  techniques  for  minimizing  the 

*  This  worlc  was  supported  by  the  National  Science  Foundation  under  Grant 
BNG-79-01686. 
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complexity  and  cost  of  the  overhead  hardware. 

Progress  has  been  reported  in  the  design  of  complex  residue  arithmetic 
for  use  in  spectral  analysis,  processing  of  bandpass  signals  represented  in 
analytic  form,  and  in  complex  cepstral  analysis  for  speech 

[ 1 3. 1 3, 1 3. SO, 1 3. 32 ].  Also  a  design  has  been  presented  which  permits  the  block 
elements  of  a  general  mixed  radix  converter  to  be  controlled  in  such  a  way 
that  the  general  converter  can  be  electronically  adapted  to  sequentially 

generate  all  projections  that  are  required  for  error  location  [ 13.261.  Since 
this  is  a  significant  step  toward  being  able  to  implement  an  entire  error 
checker  on  a  single  VLSI  integrated  circuit,  a  great  leal  of  effort  is  being 
expended  to  extend  this  concept  in  order  to  bring  failure  resistant  concepts 
into  practical  usage. 

A  new  effort  has  been  started  to  evaluate  the  advantages  and  to  extend 
the  theory  of  redundant  RNS  product  codes  for  digital  signal  processing. 
Rather  than  simply  adding  redundant  residue  digits,  data  is  encoded  in  a 
product  code  by  multiplication  with  a  cods  generator,  kn  RNS  product  code  is 
an  arithmetic  code  which  is  preserved  'under  error  free  arithmetic  operations. 
Tnile  not  too  much  is  known  about  the  theory  of  RNS  product  codes  in  signal 
processing,  initial  indications  are  that  they  have  useful  properties  worthy  of 
further  study. 

13.2  Roundoff  Noise  and  Limit  Cycles  in  Digital  Filtering* 

It  is  well  kno'wn  that  canonical  implementations  of  narrowband  recursive 
digital  filters  can  exhibit  large  error  due  to  arithmetic  roundoff.  However, 
this  error  can  be  reduced  at  the  expense  of  additional  hardware  complexity. 
For  example,  Mullis  and  Roberts  have  developed  a  procedure  for  finding  the 
minimam  ro'undoff  noise  linear  state-space  'LSSi  filter  structure  ^Rl"!. 
Another  realization  for  reducing  errors  in  digital  filters  is  the  error 
spectrum  shaping  (ESS)  structure  rR2,R3l.  The  ESS  structure  utilizes 
quantizer  feedback  and  generally  has  a  simpler  Implementation  than  the  LoS 
structure. 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
(■J.3.  Army,  H.S.  Navy,  J.3.  Air  Force)  'under  Contract  N0Q01 4-79-C-0424. 
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We  have  compared  the  noise  reduction  performance  of  the  ESS  structure 
with  that  of  the  optimal  L3S  structure  for  second  order  digital  filter 
sections.  It  was  found  that  optimal  direct  form  1  and  direct  form  2  ESS 
realizations  always  have  a  higher  signal- to-noise  ratio  than  the  optimal  L3S 
structure.  Several  suboptimal  ESS  structures  with  simple  hardware 
implementations  were  also  considered  and  optimal  values  for  the  ESS 
coefficients  in  these  structures  were  derived.  For  filters  with  zeroes  at  z  = 
-1  ,  it  was  found  that  several  of  the  simple  ESS  structures  can  outperform  the 
optimal  L3S  structure.  For  elliptic  filters  it  was  found  that  the  suboptimal 
direct  form  2  ESS  structures  perform  poorly,  but  that  the  direct  form  1  ESS 
structures  perform  well.  We  are  currently  investigating  high-order  ESS 
structures  and  developing  ESS  implementations  with  low  coefficient 
sensitivity. 

Our  worh  concerning  limit  cycles  covers  two  topics.  First,  we  have 
studied  the  concept  of  accessibility,  introduced  by  Claasen  _et  [R4-1.  A 

zero- input  limit  cycle  is  accessible  if  it  can  be  reached  from  a  state  which 
is  not  a  state  of  the  limit  cycle.  It  is  important  to  study  accessible  limit 
cycles,  since  they  are  most  frequently  observed.  We  have  shown  that 
accessibility  places  restrictions  on  the  error  sequences  introduced  at  the 
filter  quantizers.  In  the  case  of  two  rounding  quantizers,  these  restrictions 
can  be  used  to  derive  a  lower  bound  on  the  minimum  amplitude  accessible  limit 
cycle,  A  related,  but  different,  result  holds  for  two  tTrancating  quantizers. 
We  have  also  shown  that  only  certain  states  can  directly  access  a  limit  cycle 
and  that  large  limit  cycles  tend  to  be  accessible. 

The  second  limit  cycle  topic  is  that  of  determining  exact  maximum 
amplitude  limit  cycles  in  digital  filters.  We  have  recently  devised  an 
algorithm  for  finding  the  maximum  amplitude  limit  cycle  of  a  particular 
period.  This  is  accomplished  by  considering  all  possible  initial  states, 
which  can  be  maximum  amplitude  states  of  a  limit  cycle.  Results  have  been 
obtained  which  significantly  reduce  the  number  of  states  requiring 
consideration,  thereby  resulting  in  a  computationally  feasible  algorithm. 
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1 3. 3  New  Techniques  in  Di^itsl  3i^nnl  Processing  for  Synthetic  Aperture  Padnr* 

Research  in  synthetic  aperture  radar  (3AR')  has  been  conducted  in  four 
distinct,  but  related  projects:  (l )  analysis  of  inherent  phase  distortion  in 
rectangular  format  FPT  algorithms,  (2")  a  SAR  computer  simulation  vith  polar 
format  ??T  processing,  (3)  a  study  of  2D  interpolators  for  polar- to- 
rectangular  transformation,  and  (i)  an  investigation  of  number  theoretic 
concepts  for  high  speed  failure  resistant  digital  processors  required  in 
real-time  SAR.  Recently,  special  emphasis  was  placed  on  analysis  of  spotlight 
mode  operation,  including  requirements  for  2D  spatially-varying  interpolators 
for  polar- to-rectangular  coordinate  transformation  required  in  ??T  processing. 
Five  2D  interpolators  have  been  analyzed  and  simulated:  ( 1  ")  1  st  order  inverse 
distance  (nonseparabla)  ,  (2''  Ist  order  separable  bilinear.  ^3)  3rd  order 
separable  Lagrange,  (i)  separable  optimal  FIR  filtering,  and  (o'*  zero  order 
interpolation  with  over-sampling.  Attention  was  concentrated  on  low  order 
interpolators  because  they  are  simple  enough  to  implement  in  practice. 
Although  many  questions  remain,  it  appears  that  the  use  of  optimal 
space- invariant  FIR  interpolation  coupled  'with  nearest-neighbor  selection 
(zero  order  interpolation)  may  be  the  most  efficient  computational  means  to 
achieve  high  resolution  performance. 

It  was  discovered  that  in  polar  format  FFT  SAR  processing,  each  line  of 
recorded  data  is  the  Fourier  transform  of  a  density  integral  of  the  target 
area,  collected  at  a  specified  loolc  angle.  3y  interpreting  the  data  recording 
in  this  "way,  it  was  possible  to  establish  that  the  SAR  principle  does  not 
depend  on  instantaneous  motion  between  the  target  and  the  radar,  but  rather  on 
the  total  view  angle  through  which  projections  are  taken.  This  principle 
seems  to  be  obscured  in  the  present  literature,  where  SAR  is  discussed  in 
terms  of  Doppler  frequency  analysis,  thereby  implying  a  dependence  cn 
instantaneous  Doppler  frequency.  Recent  work  has  concentrated  on  learning  how 
reconstruction  can  best  be  accomplished  from  incomplete  frequency  domain  data, 
which  results  in  SAR  due  to  limited  view  angles. 


*  This  work  was  supnorted  by  the  Air  Force  Office  of  Scientific  Research  under 
Orant  AFOSR-'^O-OOT'OA. 
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1 3.4  Ney  Directions  in  Multidimensional  Signal  Processing* 

We  have  recently  begun  a  program  of  basic  research  in  multidimensional 
signal  processing  spanning  the  areas  of  digital  array  processing,  linear  space 
varying  processing,  and  nonlinear  processing  using  order  statistics.  In  array 
processing,  efforts  are  being  concentrated  on  the  development  of  a  common 
mathematical  framework  for  computer  ailed  tomography  (CAT)  and  spotlight-mode 
synthetic  aperture  radar  (3AR).  The  underlying  principle  in  CAT  is  the 
projection-slice  theorem.  We  have  shown  that  this  same  theorem  also  forms  the 
basis  for  spotlight-mode  3AR;  hence  similar  processing  is  required  in  both 
systems.  Future  work  will  explore  the  possibility  of  applying  efficient  CAT 
algoritms  to  3AR  and  vice  versa.  In  the  area  of  linear  space  varying 
processing  we  have  extended  recent  results  concerning  required  sampling  rates 
for  time  varying  digital  systems.  Unlike  the  time  invariant  case,  the 
required  sampling  rate  depends  on  both  the  input  bandwidth  and  the  desired 
analog  impulse  response.  We  are  presently  investigating  the  use  of  simple 
time  varying  digital  filters  to  perform  time  invariant  filtering.  In  the  area 
of  nonlinear  filtering  we  have  studied  a  generalized  median  filter  whose 
output  is  a  weighted  sum  of  order  statistics.  Assuming  a  constant  signal  in 
white  background  noise,  a  procedure  has  been  developed  for  finding  optimal 
weight  values.  Deterministic  properties  of  the  generalized  median  filter  and 
various  possible  recursive  order  statistical  filters  are  also  being  studied. 

1 3. 5  Image  Sequence  Processing  and  Dynamic  Scene  Analysis** 

The  processing  of  images  involving  motion  has  become  increasingly 
important.  Applications  include  target  tracking,  dynamic  robot  vision,  image 
bandwidth  compression,  dynamic  medical  imaging,  and  highway  traffic 
monitoring.  Image  sequence  processing  involves  a  large  amount  of  lata. 
However,  because  of  the  rapid  progress  in  computer,  I3I,  and  VLSI 
technologies,  many  useful  processing  tasks  for  image  sequences  can  now  be 
performed  in  a  reasonable  amount  of  time. 

*  This  -irork  was  sunnorted  by  the  Joint  Services  Electronics  Program 
(7.3.  Army,  U.S.  Navy,  IT.S.  Air  Force)  under  Contract  N0001 4-79-0-0421. 

**  This  work  was  supported  by  the  National  Science  Foundation  -under  Crant 
MCS-73-13906. 
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A  central  issue  in  im.a?e  sequence  processing  is  notion  estimation. 
During  the  past  year  xe  have  developed  a  number  of  efficient  algorithms  for 
estimating  the  three-dimensional  motion  parameters  of  rigid  bodies  'both 
translation  and  rotation)  from  two  time-sequential  image  frames.  The 
following  uniqueness  theorems  have  also  been  proven:  (i)  Por  a  rigid  planar 
patch,  two  solutions  for  the  motion  parameters  'vill  exist  in  general, 
(ii)  Por  a  rigid  curved  surface,  the  solution  for  the  motion  parameters  is 
unique . 


1  3. 5  Dnder:'fater  Image  Processing* 

■rfe  have  been  collaborating  'With  the  Naval  Research  Laboratory  on  a  basic 
research  project  in  underwater  image  processing  using  both  optical  and  sonar 
sensors.  Applications  include  searching  for  man-made  objects  on  the  ocean 
floor  and  ocean  mining.  During  the  past  year,  our  research  has  concentrated 
on  studying  the  effects  of  ocean  water  on  optical  and  sonar  images  and  on 
finding  shape  and  texture  features  of  images  of  the  ocean  floor  that  can 
distinguish  between  man-made  and  natural  objects.  Successful  algorithms  for 
extracting  shape  features  have  been  developed.  Ne  have  also  studied  the 
characteristics  of  various  acoustical  Imaging  methods  and  concluded  that 
synthetic  aperture  sonar  and  holography  are  the  two  most  promising  approaches 
for  'underwater  imaging. 
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19.1  Summary  of  Progress 

This  group  has  two  primary  areas  of  interest.  The  first  is 
computer-assisted  education  in  population  and  food  and  their  impact  on  other 
social  and  economic  areas  of  society.  The  second  is  concerned  with  climate 
and  food  and  the  impact  of  large  interannual  variations  of  climate  on  the 
availability  of  food  for  the  world's  growing  population. 

In  the  first  area,  population  and  food,  we  have  designed  a 
computer-assisted  education  program  concerned  with  world  population  as  well  as 
specific  countries.  The  program  is  available  in  3A3IT  and  has  been  adapted 
for  use  on  the  Apple  computer.  The  program  is  used  in  over  200  educational 
institutions  each  year.  Many  of  the  programs  are  used  at  the  high  school 
level  where  distribution  of  the  programs  has  just  beg-on  luring  the  past  year. 

In  the  second  area,  food  and  climate,  a  retrospective  study  has  been  made 
of  the  relationship  of  crop  production  to  climate  anomalies.  Various 
interesting  climate  sequences  have  been  observed  where  a  certain  type  of 
climate  anomaly  in  one  region  of  the  world  is  very  frequently  followed  by 
another  anomaly  in  another  part  of  the  world. 

*  This  work  was  supoorted  by  the  National  Science  Foundation  'cnder  Irant  3ED 
3501  and  United  States  Department  of  Agriculture  Contract  SEA 
5S-5t9B-0-33d. 


118 


n.  ’OPILATIOV/POOD/CLIMAT^  oT'TOIB^ 

These  sequences  of  climate  anomalies  can  be  usei  as  a  means  of  crop 
forecasting  in  the  region  of  the  later  anomaly.  Various  techniques  are  being 
ievelopei  to  attempt  to  recognize  these  sequences  in  the  TOrll  climate  pattern 
of  the  past  twenty  years,  '.fnen  special  patterns  are  obseir/si  they  are  teste! 
against  longer  lata  bases.  Some  of  the  iata  are  contains!  in  time  series 
which  go  back  to  1  365.  A  paper  iescribing  the  relationship  of  crops  in 
Australia  an!  Inlonesia  to  the  climate  anomalies  of  the  tropical  Pacific  Ocean 
has  been  accepte!  for  publication  by  the  journal  Olimate  Ohange 
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